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Abstract 


This  report  summarizes  the  work  completed  in  the  year  June  1998-May  1999  for  the  prin¬ 
cipal  investigator’s  Young  Investigator  Program  grant,  the  last  scheduled  year  for  this  three 
year  grant.  The  report  is  divided  into  five  chapters  covering  new  analytical  solutions  for  dy¬ 
namic  cracks  in  orthotropic  composites.  Each  chapter  stands  alone  and  has  been  submitted, 
accepted  or  prepared  for  publication  in  a  refereed  journal.  First,  closed  form  solutions  for 
uniformly  loaded  semi-infinite  cracks  in  orthotropic  materials  are  presented.  Next,  solutions 
regarding  the  stress  intensity  factor  at  the  tip  of  propagating  cracks  in  orthotropic  materi¬ 
als  are  presented.  In  addition,  closed  form,  dynamic  Green’s  functions  solutions  are  given 
for  semi-infinite  cracks  in  orthotropic  materials.  Previously,  in  last  year’s  report,  numerical 
Green’s  function  solutions  were  presented  for  this  case  the  first  time.  These  newer,  closed 
form  solutions  are  much  easier  to  produce  and  implement  in  more  complex  solutions.  The 
problem  of  having  a  crack  in  transversely  isotropic  material  rotated  with  respect  to  the  crack 
front  is  also  solved  with  some  effort.  The  methodology  used  in  finding  all  these  solutions  may 
be  used  in  other  problems  including  penetration  mechanics.  Such  applications  are  currently 
being  explored.  Finally,  experimental  investigation  of  the  application  of  these  solutions  to 
dynamic  crack  initiation  in  composites  under  impact  conditions  is  being  completed.  It  is 
expected  that  a  report  on  that  work  will  be  submitted  in  six  months  with  the  final  report  for 
this  contract.  (All  funds  for  the  grant  have  been  expended.  A  six  month  no  cost  extension 
has  been  granted.) 
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Chapter  1 


A  New  Method  for  Examining 
Dynamically  Loaded,  Semi-infinite 
Stationary  Cracks  in  Orthotropic 
Materials 

Co-authored  with  C.  Rubio- Gonzalez  and  accepted  for  publication  in  the  Journal  of  Mechanics 
and  Physics  of  Solids 

The  elastodynamic  response  of  an  infinite  orthotropic  material  with  a  semi-infinite  crack 
under  impact  loads  is  examined.  Three  loading  modes  are  considered,  i.e.,  opening,  in-plane 
shear  and  antiplane  shear.  Solution  for  the  stress  intensity  factor  history  around  the  crack 
tip  is  found.  Laplace  and  Fourier  transforms  along  with  the  Wiener-Hopf  technique  are 
employed  to  solve  the  displacement  formulation  of  the  equations  of  motion.  The  asymptotic 
expression  for  the  stress  near  the  crack  tip  is  analyzed  which  lead  to  a  closed-form  solution 
of  the  dynamic  stress  intensity  factor  for  each  loading  mode.  It  is  found  that  the  stress 
intensity  factors  are  proportional  to  the  square  root  of  time  as  in  the  isotropic  case.  Results 
for  orthotropic  materials  are  shown  to  converge  to  known  solutions  for  isotropic  materials 
derived  independently. 
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1.1  Introduction 


Several  problems  have  been  examined  in  dynamic  fracture  mechanics  for  isotropic  materi¬ 
als.  In  fact,  there  are  many  closed  form  solutions  for  stationary  and  propagating  cracks  in 
isotropic  materials  under  dynamic  loading  (Freund  1990,  Chen  and  Sih  1977).  The  stationary 
semi-infinite  crack  under  uniform  step  loading  in  the  crack  faces  was  considered  first  by  Maue 
(1954).  The  equivalent  propagating  case  was  studied  by  Baker  (1962).  For  orthotropic  mate¬ 
rials  the  available  solutions  are  fewer  in  number,  finite  cracks  under  impact  loading  have  been 
analyzed  by  Kassir  and  Bandyopadhyay  (1983),  Shindo  and  Nozaki  (1991),  Rubio-Gonzalez 
and  Mason  (1998)  using  integral  transform  methods.  These  authors  reduce  their  problems  to 
a  Fredholm  integral  equation  in  the  Laplace  transform  domain  which  is  solved  numerically. 
The  stress  intensity  factor  is  recovered  in  the  time  domain  by  numerical  Laplace  inversion 
of  this  solution.  Although  this  approach  can  be  used  to  solve  problems  with  non-symmetric 
loads  (Rubio-Gonzalez  and  Mason  1998)  it  is  restricted  to  finite  cracks  and  numerical  solu¬ 
tions.  With  the  growing  use  of  composites  in  many  engineering  applications  it  is  desirable 
to  have  closed  form  solutions  for  fundamental  problems  in  cracked  orthotropic  materials. 
The  semi-infinite  crack  is  a  fundamental  problem  because  cracks  under  dynamic  loads  can 
usually  be  considered  semi-infinite  for  a  short  duration  of  time  immediately  after  loading.  In 
the  present  work  the  semi-infinite  crack  problem  under  impact  loads  in  orthotropic  material 
is  analyzed,  a  closed  form  solution  for  the  dynamic  stress  intensity  factor  is  obtained  for 
each  loading  mode,  i.e.,  opening,  in-plane  shear  and  antiplane  shear.  Following  the  example 
of  Baker  (1962),  Laplace  and  Fourier  transforms  are  applied  to  the  displacement  equations 
of  motion  and  combined  with  the  Wiener-Hopf  technique  to  find  expressions  for  the  stress 
ahead  of  the  crack  tip  and  the  displacements  on  the  crack  faces.  Asymptotic  analysis  of  the 
stress  near  the  crack  tip  leads  to  expressions  for  the  stress  intensity  factors  Kj{t),  Ku(t), 
and  Km(t). 
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1.2  Governing  Equations,  in-plane  problems 


Consider  the  plane  problem  of  an  infinite  orthotropic  medium  containing  a  semi-infinite 
crack,  Figure  1.1.  Let  Ei,  ^  and  ( i,j  =  1,2,3)  be  the  engineering  elastic  constants 
of  the  material  where  the  indices  1,  2,  and  3  correspond  to  Cartesian  coordinates  (x,y,z) 
chosen  to  coincide  with  the  axes  of  material  orthotropy.  The  crack  faces  are  along  the 
negative  x-axis  and  the  origin  of  the  xy  axes  is  the  crack  tip.  Uniform  tractions  are  applied 
to  the  crack  faces  in  directions  depending  upon  the  problem  considered. 


°o  H(t) 


(a)  (b) 

Figure  1.1:  Schematic  of  the  semi-infinite  crack  geometry,  (a)  Normal  loading,  (b)  In-plane 
shear  loading. 


Restricting  the  problem  to  two  dimensions  with  wave  propagation  limited  to  the  x  —  y 
plane  by  setting  all  the  derivatives  with  respect  to  2  to  be  zero,  it  is  readily  shown  that  the 
displacement  equations  of  motion  (Nayfeh  1995)  reduce  to 


d2u  d2u  /1  ^  d2v  __  1  d2u 

Clls?  +  V  + '  +  c,2,a^  ~  ^a?’ 

d2v  d2v  ,  x  d2u  _  1  d2v 

a?  +  C22  V  +  (  12,aSa^  “  ^aF’ 


(1.1) 

(1.2) 


where  u  and  v  are  the  x  and  y  components  of  the  displacement  vector  and  Cn,  C12  and  c22 
are  non-dimensional  parameters  related  to  the  elastic  constants  by  the  relations: 


11  /x12[l  -  (E2/E1)u22Y 

C22  =  ( E2/E\)c\\ , 


C12  —  ^12C22  —  ^21  Cn, 
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for  generalized  plane  stress,  and  by 

Cl1  =  ;£k(1  - 

022  =  j^(1  ~  I'13‘/s,)’ 

E\  !  E-i  N 

Cl2  =  +  e[v'^2)’ 

A  =  1  —  ^12^21  —  ^23^32  —  ^31^13  “  ^12^23^31  —  ^13  ^21  ^32  5 


(1.4) 


for  plane  strain.  In  the  orthotropic  solid,  cs  =  \jpvil  P  represents  the  velocity  of  the  in-plane 
shear  wave  propagating  along  the  the  principal  material  axes  and  p  is  the  mass  density.  The 


stresses  are  related  to  the  displacements  by  the  equations: 


0"x 

Pl2 

° y 

P'12 

Txy 

Pl2 


du 


dv 


cn  o — b  Ci2-q-, 
ox  ay 

du  dv 

Cl2x - b  C22VT~) 

dx  dy 


du  dv 
dy  +  dx‘ 


(1.5) 


1.3  Normal  Impact 


In  the  case  of  normal  tractions  to  the  crack  faces,  figure  1.1(a),  a  spatially  uniform  pressure 
of  magnitude  cr0  is  assumed  to  be  suddenly  applied  at  t  =  0.  Exploiting  symmetry  and 
taking  only  the  upper  half  plane  y  >  0,  the  corresponding  boundary  conditions  are 

ay(x,  0,  t)  =  —a0H(t)  for  —  00  <  x  <  0, 

rxy(x,0,t)  =  0  for  -  00  <  x  <  00,  (1.6) 

v(x,  0,  t)  =  0  for  x>0, 

where  H(t)  is  the  Heaviside  step  function.  In  addition,  the  conditions  of  zero  displacements 
at  infinity  and  zero  initial  conditions  are  assumed. 

The  method  of  solution  of  the  governing  equations  presented  here  follows  that  described 
in  Freund  (1990)  for  the  isotropic  case  with  some  significant  differences;  displacement  poten¬ 
tials  are  not  used  and  no  assumptions  are  made  about  the  form  of  the  unknown  functions. 
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In  equations  (1.1)  and  (1.2),  the  time  variable  may  be  removed  by  application  of  the  Laplace 
transform 

r°o  1  r 

f*(p )  =  JQ  /(*)  dt,  f(t )  =  —  ^  f*(p)  ept  dt,  (1.7) 

where  Br  denotes  the  Bromwich  path  of  integration  which  is  a  line  parallel  to  the  imaginary 
axis  in  the  p-plane.  Applying  relations  (1.7)  to  equations  (1.1)  and  (1.2)  and  assuming 
zero  initial  conditions  for  the.  displacements  and  velocities,  the  transformed  field  equations 


become 


d2u*  d2u*  .  .  d2v 

+  +  ^  +  Cl2^ 

d2v* 


dy 2  dxdy 

d2v*  .  d2u* 

+  c22  o  o  +  (1  +  C12)- 


■2 

=  0, 


P 2 

=  °’ 
US 


(1.8) 

(1.9) 


ax2  axay 

where  the  transformed  displacement  components,  u*  and  v*,  are  now  functions  of  the  vari¬ 
ables  x,  y,  and  p.  The  application  of  the  Laplace  transform  to  the  boundary  conditions  (1.6) 


gives 


al(x,0,p) 
i~xy(x,  0  ,p) 

v*{x,0  ,p) 


-an-  for  —  00  <  x  <  0, 
P 

0  for  —  00  <  x  <  00, 

0  for  a:>  0. 


(1.10) 


To  obtain  a  solution  of  the  differential  equations  (1.8)  and  (1.9)  subject  to  conditions 
(1.10),  the  Fourier  transform  is  applied  in  this  work, 

/oo  1  roo 

f(x)  e'sx  dx,  f{x)  =  —  /  F(s)  e~isx  ds ,  (1.11) 

-00  Z7T  J—  00 

rather  than  the  double  sided  Laplace  transform.  It  is  assumed  that  the  displacements  in  the 
Laplace  transform  domain  have  the  form 

1  r°° 

u*(x,y,p)  =  —  /  A(s,y,p)e  tsx  ds,  (1.12) 

Z7T  y»oo 
1  f°° 

v*{x,y,p )  =  —  J^B{s,y,p)e  tsx  ds,  (1.13) 

where  A  and  B  are  the  Fourier  transforms  of  the  Laplace  transform  of  the  displacements, 
u*  and  v*,  respectively,  and  are  yet  to  be  determined.  Substituting  these  transforms  into 
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equations  (1.8)  and  (1.9),  the  functions  A  and  B  are  found  to  satisfy  the  simultaneous 
ordinary  differential  equations 


(cnS2  +  p2 / c2)A  -  +  (1  +  ci2)is— 


(s2  +  p2/c2s)B-c22^  +  {l+cl2)is 


fB 

dy 2 


dA 

dy 


0, 

0. 


The  solution  of  these  equations  which  vanishes  for  y  — )■  oo  is 


(1.14) 

(1.15) 


^(s,y,p)  =  Ai(s,p)e  ™  +  A2(s,p)e  l2V , 

B(s,y,p )  =  — —Ai(s,p)e~liy  -  —  A2(s,p)e~j2y, 
s  s 


(1.16) 


where  Ax  and  A2  are  arbitrary  functions  and  onj(s,p)  stands  for  the  functions 

a.(  )_c^+p2/c;-7i 

A’P>  (1  +  c12)7j  ’  J  ’ 

with  7j  and  'y2  being  two  distinct  roots  of  the  quadratic  equation 


(1.17) 


C2274  +  [(c12  +  2c12  -  CnC22)s2  -  (1  +  c22)p2/c2] 72  +  (cns2  +  p2/c2s)(s 2  +  p2/c2)  =  0.  (1.18) 


It  can  be  shown  that  for  many  materials  the  roots  71  and  72  are  real  and  positive  and  the 
expressions  for  the  displacements  in  the  Laplace  transform  domain  become: 

1  r°° 

u*  =  —  /  (Aie~™  +  A2e~™)e~tsx  ds,  (1.19) 

Z7T  J  — 00 

_ 2  /*oo  p—isx 

v*  =  —  (aiAie-™  +  a2A2e-™) - ds.  (1.20) 

Z7T  J—  00  s 

And,  using  (1.5)  the  corresponding  expression  for  r*y  is  given  by 

K,  =  ~  /j(a  1  +  Ti Mie-™  +  (Q2  +  72>2l2e-M]e-“*  ds.  (1.21) 

Applying  the  second  condition  of  (1.10)  to  equation  (1.21)  yields 


M{s,p)  =  ~PiAi(s,p), 

A  =  ^±21. 

«2  +  72 


(1.22) 
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Therefore,  the  expressions  for  the  transformed  components  of  displacement  become 

i  r°° 


u 


1  f°° 

(x,y,p)  =  —  - /31e-^)A1(s,p)e~isx  ds, 

M*>p)  u 


-oo 
—  i  roo 

v*(x,y,p )  =  —  J  -  /W"™) 

and  the  associated  stress  components  are  given  by 


ds, 


(1.23) 

(1.24) 


_* 

= 


a*y  = 


-if*  12 

2?r  ./  — oo 

-*>12  /,°° 


,4i(*»P)  -i 


/OO 

^[(ciis2  -  a;i7ici2)e-7iy  -  (cus2  -  a272c12)Ae_722/]^y^e“isx<i,.25) 
2n  L[{C”S2  ~  £*1Tic22)e_7iy  ~  (ci2S2  -  0:272^22) Ae~72y]  — 1 -  — 1  —  e~^xc(d,.26) 


T*  =  ”^12 

xy 


-00 

roo 


27T 


/oo 

(«i  +  7i)[e-712'  -  e-72l']A1(s,p)e-isids. 

-OO 


Introducing  the  functions 

E(s,p)  =  hai-Pia2)Ai(a,p)t 
s 

a! ^2  r£p' j  (p1 

F(s,p)  =  _  ^iQ,2^[C12s2  -  «l7lc22  -  A(Cl2S2  ~  <*272C22)], 


(1.27) 

(1.28) 
(1.29) 


£  = 


r{(C12  +  C12  -  CiiC22)(ci2771A^2  —  Cn) 


Cn(l  +  c12)(iV1  +  iV2) 

-C22  [cuNlNl  +  cn  (iV2  +  NiN2  +  iV2)]},  (1.30) 

N 1,2  =  (CllC22  -  cf2  -  2ci2  ±  [(cuc22  -  c\2  -  2 c12)2  -  4cnC22]1/2},  (1.31) 


where  the  velocity  Cd  =  ^/ciics  represents  the  dilatational  wave  speed  along  the  x— axis  and 
in  view  of  the  first  and  third  boundary  conditions  in  (1.10),  equations  (1.24)  and  (1.26)  yield 
the  following  pair  of  dual  integral  equations  for  the  determination  of  the  function  E(s,p) 
ip  12C  ro° 


-  F^,r!>  E(s,p)e~i,x  ds  =  -25 

p 


-  oo  <  x  <  0,(1.32) 


— %  r°° 

v*(x,  0,p)  =  —  J  E(s,p)e~isx  ds  =  0  0  <  x  <  oo.  (1.33) 

Let  v*_  (x,p)  be  the  unknown  Laplace  transform  of  the  vertical  displacement  on  the 
negative  rr— axis,  and  cr*  (x ,  p)  be  the  unknown  Laplace  transform  of  the  normal  stress  on 
the  positive  x— axis,  so  that 

»'(*,0,p)  =  l  °  f°r  x>0  „•(*  0  P)  =  /  °%(*,P)  for  *  >  0  (134) 

1  \  vl(x,p)  for  x  <  0  »l  ,  'W  \  -o-o/p  for  x  <  0.  (LMi 
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Then  the  Laplace  transform  of  the  normal  stress  and  the  vertical  displacement  on  the  whole 
boundary  y  =  0  is  given  by, 

ipnC 


I  E^v)e-"x  ds  =  r)+<4tr,p), 

2tt  7-oo  y/s2+p2/c2d  P 


/: 


r°° 

—  J^E(s,p)e~lsx  ds  =  v*_(x,p), 
and,  by  Fourier  transform  inversion,  these  equations  give 


(1.35) 

(1.36) 


7==.  g(«,p) 


V^+WI 


-^  +  E+(«), 
pis 


where 


-<B(«,P)  =  V-(>), 


roo 

E+(s)  =  /  o*+(x,p)elsxdx , 

«/  0 

,p)elsxdx, 


(1.37) 

(1.38) 


f  fca 

7-oo  \  P 


V-(s)  = 

f  v: 

7  —  00 

ei,xdx  = 

_£o_ 

ips 

(1.39) 

(1.40) 


From  the  physics  of  the  problem  it  is  reasonable  to  assume  that  the  function  <7+(a;,p) 
and  v*_{x,p)  are  exponentially  bounded  at  infinity  and  this  ensures  the  existence  of  their 
Fourier  transform  (1.39)  and  (1.40).  In  particular  it  is  shown  by  Noble  (1958)  that 

if  |<7+(a:,p)|  <  Miex~x  as  x  -»  +oo  then  E+(s)  is  analytic  in  Im(s)  =  A  >  A_  =  Ai, 

and  if  \v*_(x,p)\  <  M2ex+X  as  x  — >  — oo  then  V- (s)  is  analytic  in  Im(s)  =  A  <  A+  =  A2. 

1.3.1  Wiener-Hopf  technique 

Eliminating  E(s,p )  from  (1.37)  and  (1.38)  yields  a  Wiener-Hopf  equation 


°o  v  ,  x  c  F(s,p) 

- - E+(s)  =  Pl2C"  /  K-(s) 

Jrp2jc2 


(1.41) 


which  contains  only  the  two  unknown  functions  E+(s)  and  VL(s),  and  now  the  Wiener-Hopf 
technique  can  be  applied  as  follows.  Suppose  that  the  function  L(s)  is  defined  and  factored 
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as 


L(s)  = 


L-{s)=»ut  F{S'P) 


L+is)  \ls2+p2l 


r2 

cd 


then  equation  (1.41)  becomes 


(1.42) 


—L+(s)  -  E  +(»)£+(»)  =  M»)K.(«). 
ips 


(1.43) 


Assume  that  the  function  D(s )  can  be  decomposed  as 


%-L+(s)  =  D(s)  =  D+(s)  +  DM 
ips 


(1.44) 


then  equation  (1.43)  becomes 


D+(s)  -  E+(s)L+(S)  =  L_(s)Vl(S)  -  D.{s)  =  W(s). 


(1.45) 


The  first  member  of  this  equation  is  analytic  in  the  upper  half  of  the  s— plane  Im(s)  = 
A  >  max(Ai,A2)  <  0  and  the  second  member  is  analytic  in  the  lower  half  of  the  s— plane 
Im(s)  =  A  <  0.  Therefore,  the  regions  of  analyticity  overlap.  After  determining  W(s), 
solutions  for  E+(s)  and  VI  (s)  can  be  found. 

After  algebraic  manipulation  it  can  be  shown  that  the  function  F(s,p )  is  given  by 


F{s,p)  = 


{c22Ci27i72  —  7l72(l  +  Ci2)Ci2S2 


+ 


+p2/c2(7l  +  72)  (1  +  C12) 

(CuS2  +  p2 / c2)  [022(71  +  72)  +  7172(1  +  £-12)022  +  S2(ci2  +  cf2  ~  C11C22)  “  C-22P2 / C2]  } 


The  only  zeros  of  F(s,p)  are  of  the  form  ±ip/cR  where  cr  is  the  Rayleigh  wave  speed.  This 
can  be  seen  by  substituting  s  =  ip/v  in  F(s,p),  letting  F(s,p)  =  0  and  dividing  by  the 
non-zero  factors,  then  F(ip/v,p)  =  0  reduces  to 


2  /  C11C22  ~  c22  _  v^\ 
cn  V  c22  c2  c2\ 


1  - 


Cll^ 


=  0 


which  is  the  Rayleigh  function  for  orthotropic  materials  (Ting,  1996).  The  roots  of  this 
function  are  v  =  ±cr. 
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Im(s) 


Figure  1.2:  Branch  cuts  for  F(s)  in  the  s— plane. 


Consequently,  the  first  step  in  factoring  L(s)  is  to  define 

F(s)  =  (1.46) 

s2+p2/4 

thus  F(s)  1  as  s  -*  oo,  (the  constant  £  in  (1.29)  was  chosen  to  make  this  possible).  The 
function  F(s)  is  regular  and  F(s)  ^  0  in  the  s-plane  cut  as  shown  in  figure  1.2.  The  only 
singularities  are  the  branch  points  shared  with  71  and  72.  It  is  readily  shown  that  the  branch 
points  of  71  and  72  are 


for  71; 
for  72; 


*_P 

Cs 


s  =  ± 


ip  _  ±lp 
y/ Oil  Cd 


It  is  well-known  that  factorization  is  accomplished  most  directly  for  functions  that  ap¬ 
proach  unity  as  |s|  — »■  00  and  that  have  neither  zeros  nor  poles  in  the  finite  plane;  F(s)  is 
an  example  of  such  a  function.  Therefore,  using  Cauchy’s  integral  formula  it  can  be  shown 
that  (Freund  1990) 


/  r±  z  —  s 

where  T_  (r+)  is  the  contour  enclosing  the  branch  points  +ip/cd,  +ip/cs,  (- ip/cd ,  -ip/cs). 


Using  the  fact  that  F(s)  =  F(s)  then 


F±(s)  =  exp 


— 1  i‘1/c>  _!  ^Im[F(*pu;)]^  dw 


{  Ti 


1/Cd  W  T  ; 


is  £ 

P  J 
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Returning  to  the  factorization  of  L(s)  it  is  found  that 


therefore 


L(s)  = 


Mf)  „  cHsM+£l 4) 

^+(s)  y.S2  +  P2/Cj 

, F+(s)F-(s)(s  +  »p/ca)(ti  -  ip/cK) 
Js  +  ip/criJs  -  vpjd 


(1.47) 

(1.48) 


(1.49) 

(1.50) 

(1.51) 


with 

L+(s)  -  MO)' 
is 

Each  side  of  equation  (1.45)  is  analytic  in  one  of  the  overlapping  half  planes,  and  the  sides 
coincide  on  the  strip  of  overlap.  Consequently  each  side  of  (1.45)  is  the  analytic  continuation 
of  the  other  into  its  complementary  half  plane;  so  that  the  two  sides  together  represent  one 
and  the  same  entire  function  W  (s).  The  entire  function  will  be  determined  by  its  behavior 
at  |s|  — >  oo  which  is  related  with  the  behavior  of  physical  quantities  near  x  =  0.  First  note 
that  L+(s )  ~  s-1/2  and  L_(s)  ~  s1/2  as  \s\  — >  oo,  therefore,  D+(s)  and  £L(s)  in  (1.52) 
are  bounded  in  their  respective  planes  of  analyticity  and  vanish  at  infinity.  Furthermore, 
a*+{x, p)  is  expected  to  be  square  root  singular  as  x  — >  0+  and  v*L(x,p)  is  expected  to  vanish  as 
x  — >•  0“  to  ensure  continuity  of  displacement.  Consequently  from  the  Abel  theorem  (Noble, 
1958)  relating  asymptotic  properties  of  transforms 


D.(s)  = 


<70  ri+(o) 


(1.52) 


lim  x1l2al(x,p)  "  lim  s'^2E+(s) 

„  vn+  ‘  v  '  .s-4no  ~ v  ' 


for  some  q  >  0.  Therefore,  it  is  expected  that  E+(s)  ~  s-1/2  and  V_(s)  ~  s~1-9  as  |s|  — )■  co, 
thus  the  products  E+(s)L+(s)  and  L_(s)V_(s)  vanish  at  infinity.  Therefore,  each  side  of 
(1.45)  vanishes  as  |s|  — >  oo  in  the  corresponding  half  planes.  According  the  the  Lliouville’s 
theorem,  a  bounded  entire  function  is  constant.  In  this  case,  W(s)  is  bounded  in  the  finite 
plane  and  W(s)  0  as  |s|  -»  oo  so  that  the  constant  must  be  zero;  thus,  W(s)  =  0.  By 
using  then  (1.45)  and  (1.52),  the  functions  of  interest  are  given  by 


E+(s)  = 

V-{»)  = 


D+{') 

L+(s) 

D-(s) 


L-(s) 

1.3.2  Stress  Intensity  factor 


<7p  1  [MO) 

p  is  L+(s) 
<70  1  MO) 
p  is  L_(s)' 


-1 


(1.53) 

(1.54) 


To  find  the  stress  intensity  factor,  an  asymptotic  expression  for  the  normal  stress  near  the 
crack  tip  is  sought.  The  Abel’s  theorem  relating  asymptotic  expressions  between  a  function 
and  its  Fourier  transform  is  the  following  (Noble  1958) 

lim  \fx  al(x,p)  =  lim  E+(s).  (1.55) 

*->0+ v  +v  '  s-h-oq  V  7T 

Clearly,  the  behavior  of  E+(s)  as  s  — >  oo  is  needed.  Note  that 

L+(s)  =  as  s^°° 

and  since  F+( 0)  =  F_(0)  =  y/F( 0)  =  f^F(0,p)  then 


L+(  0) 


P 


y/wi\[F{ <Vp) 


therefore 


E+(s)  = 


<70 


P 


as  s  — ¥  oo. 


(*p)3/2  sl/2 

Using  this  relation  and  the  definition  of  the  stress  intensity  factor  in  (1.55)  gives 
Kj(p)  =  lim  V2ttx  o*Ax,p)  =  lim  e~tn^y/2s  E+(s) 

x— >0+  s— H-oo 


<7q\/2 


P 


p"’P  F(0,p) 


(1.56) 

(1.57) 

(1.58) 
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It  is  readily  shown  that 


so  that 


V'no  ,p)  = 


p\ 


I  \/&l  2 

CsCd£ 


*7(p)  = 


2cs«e  i 


(1.59) 


(1.60) 


^  ^  P3/2 

and  by  Laplace  inversion  the  dynamic  stress  intensity  factor  in  the  time  domain  for  this 
loading  mode  is 

(1.61) 

For  the  isotropic  case  and  with  plane  strain  conditions,  substituting  Ey  —  E2  =  E, 
v12  =  1^13  =  v2%  —  v  and  n,i2  =  E/ 2(1  +  v)  in  equations  (1.4)  and  the  results  in  (1.30)  it  is 
found  that  c22  =  2(1  -  v)/(l  -  2v)  and  £  =  1/(1  -  u)  and  equation  (1.61)  reduces  to  the  well 
known  result  for  isotropic  materials  (Freund  1990) 


K 


yQl  -  2^)/tt 

i(t)  =  2go  (i  —) -  vi- 


1.4  In-Plane  Shear  Loading 


Consider  the  crack  geometry  illustrated  in  figure  1.1(b).  The  crack  faces  are  subjected  to 
suddenly  applied,  spatially  uniform  shear  traction  of  magnitude  To-  Exploiting  asymmetry 
and  examining  the  half  space  y  >  0,  the  corresponding  boundary  conditions  are 


rxy(x,0,t) 
cry(x,0,t) 
u(x ,  0,  t) 


—ToH(t)  for  —  oo  <  x  <  0, 
0  for  —  oo  <  x  <  oo, 

0  for  x  >  0, 


(1.62) 


where  H (t)  is  the  Heaviside  step  function.  In  addition  the  condition  of  zero  displacements 
at  infinity  and  zero  initial  conditions  are  assumed. 
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The  method  of  solution  is  similar  to  that  used  for  normal  impact.  Applying  Laplace 
transform  to  boundary  conditions  (1.62)  gives 

r*y(x,0,p)  =  -r0/p  for  —  oo  <  x  <  0, 

<7*(x,0,p)  =  0  for  —  oo  <  x  <  oo,  (1.63) 

u*(x,0,p)  —  0  for  x  >  0. 


Assuming  the  displacement  field  (1.12)  and  (1.13)  yields  the  same  system  of  ordinary  differ¬ 
ential  equations  (1.14)  and  (1.15)  whose  solution  is  given  by  equations  (1.16)  where  ctj  are 
defined  by  (1.17)  and  7 y  are  obtained  from  the  solution  of  (1.18).  Using  (1.19),  (1.20)  and 
(1.5)  indicates  that  a*  is  given  by 

%  j[/i  o  poo 

<rfa,y,p)  =  -  —  /  [(s2ci2-ai7iC22)A1e-7lV  +  (s2Ci2-a272C22)A2e-72!'] - ds.  (1.64) 

Z7T  J  —  OO  S 

Applying  the  second  condition  of  (1.63)  to  equation  (1.64)  yields 


A2(s,p) 

A 


-/Mi  (a,p), 
s2c12  -  o;i7iC22 
S2Ci2  —  «272c22 


(1.65) 


Therefore,  the  expressions  for  the  transformed  components  of  displacement  become 


1  r°° 

u*{x,y,p )  =  —  j^{e-'yiy  - /32e~^y)A1(s,p)e~tsx  ds,  (1.66) 

v*(x,y,p)  =  —^  [  (a1e~7iy  -  p2a2e~'r2y)—  ^,p\~Ux  ds,  (1.67) 

27 r  J- oo  s 


and  the  associated  stress  components  are  given  by 

<  =  r[(cnS2-cv17iC12)e-7^-(CllS2-a272Ci2)/32e-72y]^^e-i-<468) 

Z7T  J  —oo  S 

v*y  =  Cjc^s2  ~  «i7ic22)(e-7lj/  -  e-™)Al(-Ss’Ph-isxds,  (1.69) 

<y  =  ^/_J(ai+7i)e-7iy-(a2  +  72)/?2e-72j/]A1(S,p)e-is3;dS.  (1.70) 

Introducing  the  functions 


C(s,p )  =  (1  —  P2)Ai(s,p), 


(1.71) 
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(1.72) 


G(s,p )  =  ^  _j£y"  [(Qii  +  7i)  ~  ft  (<*2  +  72)1, 

77  =  Cl1^  (1.73) 

where  £  and  iV1>2  are  given  in  (1.30)  and  (1.31)  and  in  view  of  the  first  and  third  boundary 
conditions  in  (1.63),  equation  (1.66)  and  (1.70)  yield  the  following  pair  of  dual  integral 
equations  for  the  determination  of  the  function  C(s,p), 

-P12V  ro° 


P) 


7: 


C(s,p)e~isx  ds  =  --  -  00  <  a;  <  0(1.74) 


V 


u 


2  7T  J-00  ^s2  +  p2/ 

1  r°o 

*(x,0,p)  =  w-  C(s,p)e~lsx  ds  =  0  0  <  a:  <  00,  (1.75) 

27T  7—oo 

Let  u*_(x,p )  be  the  unknown  Laplace  transform  of  the  displacement  on  the  negative 
x— axis,  and  r£(x,p)  be  the  unknown  Laplace  transform  of  the  shear  stress  on  the  positive 
x— axis,  so  that 

*/«  \  fO  for  x  >  0  */'n  \  (  T+(xtP )  f°r  ^  >  0  (-1  na\ 

u  (x,  0,p)  -  |  for  x<o  Txy(xAp)~  |  _To/p  for  £  <  0.  (  *76^ 

Then,  the  Laplace  transform  of  the  shear  stress  and  the  displacement  on  the  whole  boundary 
y  —  0  is 

-P12V  f°° 

J — 00  .  /o2  _L  <r>2  /^»2  P 

-isa:ds  =  7il(x) 


2n 


1  r°° 

j  C(S,v)e~' 


(1.77) 

(1.78) 


and  by  Fourier  transform  inversion,  these  equations  give 

,G(S’?).  :C{s,p)  =  -5-  +  7VW 


\/s2  +p2A 


C(«,p)  =  £7-(«) 


where 


r;(x)eisidx 


r+w  =  r 

./0 

C/_(s)  =  f°  u*_(x)eisxdx 
J  —  OO 


(1.79) 

(1.80) 

(1.81) 

(1.82) 


/ 


-To 


eisxdx  = 


-00  V  p 


—  OO 

To 

ips 
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Eliminating  C(s,p)  from  (1.79)  and  (1.80)  results  in  a  Wiener-Hopf  equation 


ro  m  _  ..  _  G(S’P)  tt 

T+(s)  —  fj,l2rj—7====U-(s) 


(1.83) 


ips  \fs2  +  p2/Cs 

which  contains  only  the  two  unknown  functions,  T+(s)  and  C/_(s),  and  can  be  solved  using 
the  Wiener-Hopf  technique  as  in  the  normal  impact  analysis.  The  result  is 


where 


D+(s)  -  T+(s)L+(s )  =  L-(s)U-(s)  -  £>_(s)  =  0 


r  (  \  AS  -*PM/>,  t  ^ 

L-{s)  =  -put  ,  —  G-(s), 


L+(s) 


s  ip/cs 

+  ip/cs  1 
(s  +  ip/cR)  G+(s)  ’ 


and 


with 


where 


D{s)  =  £l+(s)  =  ^ 
ips  p 


L+(s)  —  L+(0)  ^  L+( 0) 


is 


is 


D+(s)  =  Tj 


L+(s)  —  L+(  0) 


is 


and  D _  (s) 


=  D+(s)  +  T>_(s), 

L+  (0) 


n 

p 


IS 


G(s)  = 


G(s,p) 


s2+p2/c2R 

It  can  be  shown  that  G(s)  -»  1  as  s  — >  oo.  Thus,  the  solution  of  (1.84)  is 


T+(s ) 
C/_(S) 


g±(f) 

L+(s) 


TO 

p  is 


L+(0) 


L+(s) 


D-(s)  T0  1  L+(0) 


(1.84) 


(1.85) 

(1.86) 

(1.87) 


L-(s)  p  is  L-(s)' 

with  L±(s)  defined  in  (1.85)  and  (1.86)  and  G±(s)  given  by 

G±(s)  =  exp  (=1  /1/C*  tan"1  -  , . 

(  7T  J\/cd  \Re[G(ipw;)]/ 

To  calculate  the  stress  intensity  factor  T+(s )  as  s  — >■  oo  is  needed.  First,  note  that 

L+(s)  =  — ^  as  s  y  oo 


(1.88) 

(1.89) 

(1.90) 
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and  since  G+(0)  =  G_(0)  =  ^<5(0)  =  ^ijG(0,p)  then 


L+(0)  = 


Vwc»yjG{  0,p) 


Therefore, 


A00  =  77 


-T0  1 


P 


as  s  — ¥  oo. 


(zp)3/2  s1/2  ^G(0,p) 

Using  (1.55)  and  the  definition  of  the  stress  intensity  factor  gives 


Kjr(p)  =  lim  \f2nx  rl(x,p )  =  lim  e  l7r/4-\/2s  T+(s) 

iiv  '  I-J-0+  +  s-H-oo  ' 

roV2  p 


It  can  be  shown  that 


so  that 


and  by  Laplace  inversion 


s/g(  o.p)  = 


p 


*vVp)  =  ToV/2c~P  -J7J, 


P 


(1.91) 

(1.92) 

(1.93) 

(1.94) 

(1.95) 


(1.96) 


For  the  isotropic  case  and  with  plane  strain  conditions  77  =  1/(1  —  and  equation  (1.96) 
reduces  to  the  well  known  result  for  isotropic  materials  (Freund  1990) 


Kn(t)  =  27oi 


2c, 


7r(l  —  v) 


Vi. 


1.5  Antiplane  Shear 


Consider  a  semi-infinite  crack  in  orthotropic  material  shown  in  figure  1.3,  subjected  to  an¬ 
tiplane  shear  loading  of  magnitude  r0.  The  governing  differential  equation  for  two  dimen¬ 
sional  antiplane  motions  of  homogeneous,  orthotropic,  linearly  elastic  solids  is  (Nayfeh  1995) 


2d2 w  d2w  _  1  d2w 

dx2  +  dy2  c\  dt2  ’ 


(1.97) 
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where 


32  C55 


P  -  75-  and  ch 

O44 


—  \JCu/p- 


Again,  the  principal  material  axes  of  orthotropy  are  chosen  to  coincide  with  the  ( x ,  y,  z)  axes, 
such  that  the  in-plane  and  antiplane  motions  are  not  coupled.  In  terms  of  the  engineering 
constants,  C44  and  C55  are  given  by  C44  =  //23  and  C55  =  p 31.  The  non-zero  stresses  are 
related  to  the  displacements  by  the  equations 


dw 

OX 


Tyz  =  C44^. 


dw 

dy 


(1.98) 


Exploiting  symmetry  and  examining  the  upper  half-plane  y  >  0,  the  boundary  conditions 
for  this  loading  mode  are  given  as 


ryz(x,0,t)  =  —T0H(t)  for  x  <  0 
w(x,  0,  i)  =  0  for  x  >  0. 


(1.99) 

(1.100) 


t0H(t) 

Figure  1.3:  Schematic  of  the  semi-infinite  crack  geometry.  Antiplane  shear  loading. 


The  method  of  solution  is  very  similar  to  that  used  for  the  in-plane  problems.  Applying 
Laplace  transform  to  the  governing  equation  gives 


2  d2w*  d2w* 
dx2  '  dy2 

0 

II 

* 

1 

(1.101) 

and  to  the  boundary  conditions  gives 

Tyz(x^P)  =  ~To/P 

for  x  <  0 

(1.102) 

w*(x,0,p)  =  0 

for  x  >  0. 

(1.103) 
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Assuming  a  displacement  field  of  the  form 

1  r°° 


w 


i  r°° 

(x,y,p)  =  —  /  A(s,y,p)e  tsx  ds, 

Z7T  J— oo 


(1.104) 


and  inserting  it  into  (1.101)  results  in  the  ordinary  differential  equation 

cPA 


dy 2 


(s2/32+p2/c2)A  =  0. 


The  solution  of  this  equation,  which  vanishes  for  y  — >  oo  is 


A{s,y,p)  =  A1(s,p)e  iy, 


with 


7=  \//?2s2  +P2/C1- 


Therefore,  the  non-zero  displacements  and  stresses  are 

1  r00 


w 


(x,y,p)  =  -  ds, 

—  iClcc  fOO 

J  s  A\{s,p)e~lve~tsx  ds, 


tL(x’V,p)  = 

Tyz(X^y,P)  = 


27T  J — oo 

-<?44  '°° 
27T 


/oo 

7  A1(s,p)e-7S/e-“x  ds. 

-00 


(1.105) 

(1.106) 

(1.107) 

(1.108) 


Applying  the  boundary  conditions  (1.102)  and  (1.103)  yields  a  system  of  dual  integral 
equations  for  the  determination  of  the  function  Ai(s,p) 

y°°  A1(s,p)e~'yye-isx  ds  =  0  for  x  >  0, 

/OO 

'y  A\{s,p)e~'iye~lsx  ds  =  —  r0/p  for  x  <  0, 

-OO 


-cu  r°° 

27 r 


or  on  the  whole  boundary  y  =  0 

1  f°° 


—  T  A1(s,p)e^ye~isx  ds  =  -w*_ (x,p), 

Z7T  J—  00 

f°°  7  Ai(s,p)e~'yye~isx  ds  =  --H{-x)  +  r;(x,p). 

d — OO  P 


-C44 

27T  J  —oo 


By  Fourier  transform  inversion 


Ai(s,p)  =  W_(s), 


-CulfAi(s,p)  =  -^-+T+(s), 
ips 


(1.109) 

(1.110) 
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where 


W _ 


'-(s)  =  [  w*_(x,p)etsxdx, 

J —00 

roo 

T+(s )  =  /  r?  (x,p)etsxdx. 

Jo 

Eliminating  Ax  from  equations  (1.109)  and  (1.110)  gives  the  Wiener-Hopf  equation 

^-T+(s)  =  CuW-(s). 

ips 

Using  the  Wiener-Hopf  technique  as  in  the  previous  problems  it  is  found  that 


r+(»)  =  — ~7 


Mg)  ; 

***  /  \ 


where 


and 


M»)  = 


7+(«) 


p  is  [£-($) 

M»)  =  c«7-(»). 


7+(s)  =  \[iJs  +  ip Jcii,  7-(s)  =  -  ip/ch. 

To  find  the  stress  intensity  factor  T+(s)  as  s  M  oo  is  needed.  This  is  given  by 


rp  I  \  T0  \/ [lCh 

i+(s)  =  Tin  as  s-¥  oo. 


(*p)3/2  Sl/2 

Therefore,  the  dynamic  stress  intensity  factor  is 

W 2 


V 


and  in  the  time  domain 


(l.in) 


(1.112) 


(1.113) 


(1.114) 


For  transversely  isotropic  materials  with  fibers  along  the  x— axis,  C55  =  pu  and  (3ch  — 
yfpn/p  —  cs,  this  equation  becomes 

K m(t)  =  2T0^2cs/i{\/t. 

For  isotopic  materials,  the  same  expression  for  Kuj(t ),  with  the  shear  wave  velocity  cs  = 
\Jp/p  is  found. 
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Stress  Intensity  Factor,  Semi-infinite  crack 


Figure  1.4:  Stress  intensity  factor  history  for  a  semi-infinite  crack  in  different  materials,  real 
time. 


Graphite 

Epoxy 

E-Glass 

Epoxy 

Boron 

Epoxy 

Epoxy 

Steel 

Cll 

20.77 

8.38 

32.67 

3.91 

C22 

2.29 

3.91 

3.5 

Cl2 

1.91 

1.5 

7.48 

5.5 

6.4 

1.96 

76.92 

Began 

wmmaiM 

0.73 

nm 

KB 

Table  1.1:  Mechanical  properties  used  for  the  analysis. 


1.6  Results  and  Conclusions 

Closed  form  expressions  of  the  dynamic  stress  intensity  factors  Ki(t),  Ku(t)  and  Km(t) 
have  been  determined  for  semi-infinite  cracks  in  orthotropic  materials  under  the  three  loading 
modes,  i.e.,  opening,  in-plane  shear  and  antiplane  shear  loadings,  respectively.  The  method  of 
solution  differs  from  that  typically  used  in  the  isotropic  case.  The  displacement  formulation 
of  the  equations  of  motion  is  solved  without  the  use  of  Helmholtz  potentials.  It  has  been 
shown  that  the  orthotropic  formulation  includes  the  isotropic  results  as  special  cases.  That  is, 
for  each  loading  mode  the  dynamic  stress  intensity  factor  for  isotropic  materials  is  recovered 
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Semi-infinite  Crack,  Normal  and  Shear  Loading 


Figure  1.5:  Stress  intensity  factor  history  Kj(t )  and  Kn(t)  for  a  semi-infinite  crack  in  differ¬ 
ent  materials.  The  isotropic  case  corresponds  to  epoxy  and  the  orthotropic  to  graphite-epoxy. 

from  the  orthotropic  expressions  with  the  proper  substitution  of  the  elastic  constants. 

Figure  1.4  shows  the  stress  intensity  factor  history  for  the  opening  mode  for  different 
materials.  The  material  properties  used  are  shown  in  table  1.1.  As  in  the  isotropic  case, 
the  stress  intensity  factors  are  proportional  to  the  square  root  of  time;  hence  there  is  no 
equivalent  quasi-static  problem  for  the  semi-infinite  crack  under  uniform  loadings  and  there 
is  no  long-time  equilibrium  value.  It  can  be  seen  that  the  fiber  reinforcement  of  epoxy  leads 
to  an  increase  of  the  stress  intensity  factor  for  a  given  time,  i.e.,  Kj(t)  for  graphite-epoxy 
composite  is  greater  than  the  corresponding  Kj{t)  for  epoxy  at  a  given  time  t.  Figure  1.5 
shows  the  stress  intensity  factors  Ki(t )  and  Ku(t )  for  isotropic  (epoxy)  and  orthotropic 
(graphite-epoxy)  materials,  note  that  with  same  load  amplitude  the  stress  intensity  factor 
Kn(t)  is  greater  than  Kj{t)  at  a  given  time.  In  both  cases,  the  introduction  of  fibers  results 
in  an  increase  in  the  stress  intensity  factor. 

It  is  interesting  to  note  that  Ki(t)  for  a  composite  with  fibers  parallel  to  the  x— axis 
and  for  a  composite  with  fibers  parallel  to  y-axis  is  the  same.  This  result  is  contained  in 
equation  (1.61)  since  the  ratio  £/ yfon  is  invariant  to  the  replacement  of  C22  by  cu  and  Cn 
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Finite  Crack,  Orthotropic  Material,  Uniform  Load 


Figure  1.6:  Comparison  of  the  stress  intensity  factor  Kj(t)  for  a  semi-infinite  crack  and  a 
finite  crack  of  length  2 a  in  orthotropic  material,  the  latter  is  a  numerical  solution  obtained  us¬ 
ing  integral  transform  methods  by  Kassir  and  Bandyopadhyay  (1983)  and  by  Rubio-Gonzalez 
and  Mason  (1998)  using  FFT  for  Laplace  transform  inversion.  K0  =  o^s/va. 

by  C22.  This  result  does  not  hold  for  the  in-plane  shear  loading  however,  since  Kji(t)  for  the 
fibers  oriented  along  the  y— axis  is  less  than  that  for  fibers  oriented  along  the  re— axis  by  a 
factor  of  (C22/C11)1/4.  This  factor  is  given  in  Table  1.1  for  the  materials  examined  here. 

Integral  transform  methods  have  been  applied  to  solve  dynamic  problems  dealing  with 
finite  cracks  in  orthotropic  materials.  Figure  1.6  shows  the  stress  intensity  factor  for  a  finite 
crack  of  length  2a  in  orthotropic  material  under  uniform  impact  loading,  this  is  an  approxi¬ 
mated  solution  obtained  by  Kassir  and  Bandyopadhyay  (1983)  and  by  Rubio-Gonzalez  and 
Mason  (1998)  using  a  more  suitable  technique  for  numerical  inversion  of  the  Laplace  trans¬ 
form.  It  is  known  that  for  short  times  the  finite  crack  behaves  like  a  semi-infinite  crack. 
As  illustrated  in  the  figure,  the  closed  form  solution  derived  here  agrees  with  the  numerical 
solution  for  short  times  ( cst/a  <  2)  as  expected. 
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Chapter  2 


Dynamic  Stress  Intensity  Factor  for  a 
Propagating  Semi-infinite  Crack  in 
Orthotropic  Materials. 


co-authored  with  C.  Rubio- Gonzalez 

The  elastodynamic  response  of  an  infinite  orthotropic  material  with  a  semi-infinite  crack 
propagating  at  constant  speed  is  examined.  Solution  for  the  stress  intensity  factor  history 
around  the  crack  tip  is  found  for  the  loading  modes  I  and  II.  Laplace  and  Fourier  transforms 
along  with  the  Wiener-Hopf  technique  are  employed  to  solve  the  equations  of  motion.  The 
asymptotic  expression  for  the  stress  near  the  crack  tip  is  analyzed  which  lead  to  a  closed-form 
solution  of  the  dynamic  stress  intensity  factor.  It  is  found  that  the  stress  intensity  factor  for 
the  propagating  crack  is  proportional  to  the  stress  intensity  factor  for  a  stationary  crack  by 
a  factor  similar  to  the  universal  function  k(v)  from  the  isotropic  case.  Results  are  presented 
for  orthotropic  materials  as  well  as  for  the  isotropic  case. 
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2.1  Introduction 

Problems  of  crack  propagation  at  constant  speed  can  be  classified  into  three  classes  depending 
on  the  boundary  conditions  [5].  The  first  class  is  the  steady  state  crack  growth.  Here,  the 
crack  tip  moves  at  constant  speed  for  all  the  time  and  the  mechanical  fields  are  invariant 
with  respect  to  an  observer  moving  with  the  crack  tip.  The  second  class  of  problems  is  the 
self-similar  crack  growth  subjected  to  time-independent  loading.  In  this  case,  the  crack  tip 
moves  at  constant  speed  from  some  initial  instant,  and  certain  mechanical  fields  are  invariant 
with  respect  to  an  observer  moving  steadily  away  from  the  process  being  observed.  The  third 
category  of  problems  corresponds  crack  growth  due  to  time-dependent  loading  usually  in  the 
form  of  crack  face  pressure  or  normally  incident  stress  pulse. 

The  prototype  problem  of  the  first  category  is  the  two  dimensional  Yoffe  problem  [15]  of 
a  crack  of  fixed  length  propagating  in  an  isotropic  body  subjected  to  uniform  remote  tensile 
loading.  Even  though  this  mathematical  problem  is  not  a  realistic  model  of  a  physical 
situation  because  of  the  feature  that  the  crack  closes  at  one  end  at  the  same  rate  at  which 
it  opens  at  the  other  end;  it  is  important  because  some  field  quantities  are  independent 
of  the  fictitious  crack  length  2a,  in  particular  the  angular  variation  of  the  near  tip  stress 
field.  Propagating  finite  cracks  under  steady  state  conditions  have  also  been  considered 
for  orthotropic  materials  by  several  authors  using  different  techniques.  Kassir  and  Tse  [6] 
through  an  integral  transform  technique  reduced  the  related  boundary  value  problem  to  a 
system  of  dual  integral  equations.  Piva  and  Viola  [11]  solved  the  same  problem  through  a 
complex  variable  approach  reducing  the  equations  of  motion  to  a  first  order  elliptic  system  of 
the  Cauchy-Riemann  type.  More  recently,  Lee,  et  al  [8]  obtained  higher  order  expressions  for 
the  stress  and  displacements  around  the  tip  of  a  crack  propagating  in  orthotropic  materials. 
In  that  work  a  power  series  representation  of  complex  potentials  was  used.  Atkinson  [1] 
studied  the  steady-state  propagation  of  a  semi-infinite  crack  in  aelotropic  materials  by  means 
of  the  Cauchy  integral  formula.  In  all  of  these  works  the  emphasis  was  in  the  angular 
variation  of  the  stress  field  near  the  crack  tip.  Kousiounelos  and  Williams  [7]  analyzed  the 
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problem  of  a  dynamically  propagating  crack  in  a  highly  orthotropic  fiber  composite  infinite 
strip  subjected  to  constant  displacement  mode  I  loading.  The  basic  assumption  was  to 
neglect  the  strain  along  the  fibers  and  simplify  substantially  the  mathematical  treatment. 
Using  Fourier  transforms  and  the  Wiener-Hopf  technique  they  obtained  expressions  for  the 
dynamic  stress  intensity  factor  and  dynamic  energy  release  rate  under  this  simplified  model. 

The  Broberg  problem  [3]  is  the  prototype  of  the  second  category.  In  self-similar  problems 
attention  is  focused  on  the  symmetric  expansion  of  a  crack  at  constant  rate  from  zero  initial 
length.  A  thorough  discussion  of  the  analysis  of  self-similar  mixed  boundary  value  problems 
in  elastodynamics  was  presented  by  Willis  [14].  He  considered  problems  in  both  two  and 
three  space  dimensions  as  well  as  isotropic  and  anisotropic  materials. 

The  prototype  of  the  third  category  is  the  Baker  problem  [2]  for  isotropic  materials.  He 
considered  a  transient  problem  in  which  a  semi-infinite  crack  extends  at  constant  velocity 
after  a  step  stress  loading  on  the  crack  faces  is  applied  at  time  t  =  0.  This  type  of  loading 
induces  a  transient  traction  distribution  on  the  crack  plane  ahead  of  the  crack  tip.  If  the 
material  is  of  limited  strength,  the  crack  will  begin  to  extend  at  some  later  time,  t  =  to,  after 
the  application  of  the  load,  where  to  Is  the  delay  time  for  the  process.  Freund  [4]  showed 
that  the  stress  intensity  factor  is  independent  of  the  time  delay  and  that  it  has  the  form  of 
the  stress  intensity  factor  for  a  stationary  crack  subjected  to  the  same  stress  wave  loading 
multiplied  by  the  stress  universal  function  ki(v),  i.e. 

KI(t,v)  =  kI(v)KI(t,0).  (2.1) 

Furthermore,  Freund  [4]  also  demonstrated  that  the  relation  (2.1)  holds  for  crack  propagation 
at  non-uniform  speed  being  v  the  instantaneous  crack  tip  velocity  and  Kj(t,  0)  the  stress 
intensity  factor  for  the  crack  as  it  had  been  at  the  instantaneous  position  for  all  time. 
Note  the  importance  of  the  universal  function  of  the  crack  tip  velocity  ki(v)  in  problems  of 
arbitrary  loading  and  non-uniform  crack  speed,  even  though  ki(v)  could  have  been  obtained 
in  the  solution  of  a  less  general  problem  like  that  of  a  semi-infinite  crack  propagating  at 
constant  velocity  and  under  a  step  loading  (Baker  problem). 
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The  extension  of  the  Baker  problem  to  orthotropic  materials  under  mode  I  and  mode  II 
loading  is  the  aim  of  this  work.  That  is,  the  problem  to  be  examined  consists  of  an  unloaded 
crack  propagating  at  constant  velocity  v  in  an  undeformed  unloaded  orthotropic  material 
that  remains  at  rest.  At  some  time,  t  =  0,  uniform  stress  is  applied  to  the  crack  faces.  The 
resulting  history  of  the  stress  intensity  factor  at  the  crack  tip  is  a  function  of  the  velocity 
as  well  as  the  material  properties  and  time,  Ki  =  Ki(t,Cij,v)  where  i  =1,  II.  As  will  be 
shown  here,  the  dynamic  stress  intensity  factor  for  the  stationary  semi-infinite  crack  due  to 
Rubio-Gonzalez  and  Mason  [12]  is  recovered  letting  v  =  0.  In  [12]  it  was  shown  that 


A/M) 

=  2cr0^ 

2cs£(0)  n 

, —  v  *-> 

7T\/C22 

(2.2) 

A//M) 

=  2t0i/2c577(0)/7t  Vi , 

(2.3) 

which  represent  one  of  the  few  available  closed  form  solutions  for  the  dynamic  stress  intensity 
factors  for  cracks  in  orthotropic  materials.  Making  use  of  these  results,  it  will  be  possible 
to  write  Ki(t,v )  and  Kn(t,v)  in  the  form  of  (2.1)  and  consequently  derive  expressions  for 
the  universal  functions  of  the  crack  speed  ki(v )  and  kn(v)  in  orthotropic  materials.  The 
motivation  to  achieve  this  task  is  that  kj(v)  and  kn(v)  may  be  useful  in  more  general 
problems  as  in  the  case  of  isotropic  materials. 

2.2  Governing  Equations,  in-plane  problems 

Consider  the  plane  problem  of  an  infinite  orthotropic  medium  containing  a  semi-infinite 
crack,  Figure  2.1,  propagating  at  constant  velocity,  v  along  the  re'— axis.  Let  Ei,  Hij  and 
Uij  ( i,j  —  1,2,3)  be  the  engineering  elastic  constants  of  the  material  where  the  indices  1, 
2,  and  3  correspond  the  Cartesian  coordinates  ( x y',  z ')  chosen  to  coincide  with  the  axes  of 
material  orthotropy.  The  crack  faces  are  parallel  to  x'— axis. 

Restricting  the  problem  to  two  dimensions  with  wave  propagation  limited  to  the  x'  —  y' 
plane  by  setting  all  the  derivatives  with  respect  to  z  to  be  zero,  it  is  readily  shown  that  the 
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Figure  2.1:  Schematic  of  the  semi-infinite  crack  propagating  at  constant  speed. 


displacement  equations  of  motion  [9]  reduce  to 

d2ux  d2ux  d2uy 

ciiTrT2-  +  TnT  +  (l  +  ci2) - 


_  1  d2ux 

dx'2  '  dy'2  +  ^  °12^  dx'dy'  c2  dt 2 

|  _  .  / -i  i  .  \  dV  _ 

<9z'2  22  dy'2  12  dx'dy'  c2  dt 2 


(2.4) 

(2.5) 


where  ttx  and  uy  are  the  a:  and  y  components  of  the  displacement  vector  and  cn,  ci2  and  c22 
are  non-dimensional  parameters  related  to  the  elastic  constants  by  the  relations: 

Ei 


Cll  = 


/r12[l  -  {E2/Ei)u22y 
c22  =  {E2/Ei)cu, 


(2.6) 


C 12 

for  generalized  plane  stress,  and  by 

Ei 


Vl2C22  =  l^2i  On, 


Cu  = 


(1  —  V2zV$2), 


p 12^ 

E2 

C22  =  -^-(1 -^13^31), 


C12  = 


//12A 

Ei  (  ,  e2 

T^A^21  +  E-W^’ 


(2.7) 


A  —  1  —  ^12^21  —  ^23^32  —  ~  ^12^23^31  —  v13l'2\vZ2^ 


for  plane  strain.  In  the  orthotropic  solid,  cs  =  \Jp\2/ p  represents  the  velocity  of  the  in-plane 
shear  wave  propagating  along  the  the  principal  material  axes  and  p  is  the  mass  density. 

A  new  coordinate  system  (x,  y)  is  attached  to  the  propagating  crack  tip,  the  relation 
between  the  fixed  and  moving  coordinates  are 


x  =  x'  —  vt  y  =  y' . 
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Steady-state  conditions  are  not  assumed,  consequently  the  motion  equations  in  the  moving 
coordinates  are 


f  v2\  d2ux 


d2Uy 

dx 2 


+  c22 


d^Ux 
dy 2 

d2Uy 

dy 2 


+  (1  +  ci2) 


d2Uy 

dxdy 


+  (1  +  cia) 


d2ux 

dxdy 


—  1 

(  d2ux 

c2 

^  dt 2 

. 

dxdt  J  ’ 

1 

(  d2Uy 

2v8'2uA, 

T2 

\  dt 2 

dxdt)’ 

The  stresses  are  related  to  the  displacements  by  the  equations: 


g X 

H12 

Ml2 

Txy 

Vl2 


dux 

dx 

dux 


cn-z~  +  ci2-a  > 
dx  dy 


dUy 

'  dy 
du„ 


Cl  2T; - 1"  c22' 


5a: 


dy ’ 


5ux 

5y 


+ 


5% 

dx 


(2-8) 

(2.9) 


(2.10) 


2.3  Normal  Impact 


A  spatially  uniform  pressure  of  magnitude  a0  is  applied  suddenly  on  the  crack  faces  at 
t  =  0.  Exploiting  symmetry  and  taking  only  the  upper  half  plane  y  >  0,  the  corresponding 
boundary  conditions  are 


<7y(x,  0,t)  —  -cr0H(t)  for  -  00  <  x  <  0, 

rxy(x,0,t)  =  0  for  -  00  <  x  <  00,  (2-11) 

uy(x,  0,  t)  =  0  for  x  >  0, 

where  H(t)  is  the  Heaviside  step  function.  In  addition  the  condition  of  zero  displacements 
at  infinity  and  zero  initial  conditions  are  assumed. 

The  method  of  solution  of  the  governing  equations  presented  here  follows  that  described 
by  Baker  [2]  for  the  isotropic  case  with  some  significant  differences;  displacement  potentials 
are  not  used.  In  equations  (2.8)  and  (2.9),  the  time  variable  may  be  removed  by  application 
of  the  Laplace  transform 

r(p)  =  r  m  e-p‘  m,  m  =  Z-.[  rw)  <*,  (2.12) 

Jo  Zm  J  Br 
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where  Br  denotes  the  Bromwich  path  of  integration  which  is  a  line  parallel  to  the  imaginary 
axis  in  the  p-plane.  Applying  relations  (2.12)  to  equations  (2.8)  and  (2.9)  and  assuming 


zero  initial  conditions  for  the  displacements  and  velocities,  the  transformed  field  equations 
become 


Cn - o  I  Cu 


d2u*x  d2u*x  „  .  d2u*y 

— r  +  ^f  +  (l+ci2) - y~ 


dx 2 


1  - 


d2K 

dx2 


dy2 

d2u 


+  °22~dy2  +  ^  +  °12^ 


dxdy 

d2ul 


4(?2<-2,p^ 


1 


dxdy 


\  [p2<  ~  2vpih 


=  0,  (2.13) 

=  0,  (2.14) 


where  the  transformed  displacement  components,  ux  and  u*,  are  now  functions  of  the  vari¬ 
ables  x,  y,  and  p.  The  application  of  the  Laplace  transform  to  the  boundary  conditions 


(2.11)  gives 


a*y(x,0,p) 

Txy(XT°’P) 

U*y(X,0,p) 


=  — cr0-  for  —  oo  <  x  <  0, 

P 

=  0  for  —  oo  <  x  <  oo, 
=  0  for  x  >  0. 


(2-15) 


To  obtain  a  solution  of  the  differential  equations  (2.13)  and  (2.14)  subject  to  conditions 
(2.15),  the  Fourier  transform  is  applied, 

/oo  1  poo 

f(x)eisxdx,  f(x)  =  —  F(s)  e~isx  ds.  (2.16) 

-oo  L'K  J — oo 

It  is  assumed  that  the  displacements  in  the  Laplace  transform  domain  have  the  form 


u. 


u. 


1  r°° 

■.(x,y,P)  =  tsx  ds, 

1  r°° 

,(x,V,v)  -  ^  J  B(s,y,p)e  ds, 


(2.17) 

(2.18) 


where  A  and  B  are  the  Fourier  transforms  of  the  Laplace  transform  of  the  displacements, 
u*x  and  u*,  respectively,  and  are  yet  to  be  determined.  Substituting  these  transforms  into 
equations  (2.13)  and  (2.14),  the  functions  A  and  B  are  found  to  satisfy  the  simultaneous 
ordinary  differential  equations 


cn - o  s  + 


2  vpis  p2 


+  — 
r2  ^  „2 
cs  cs  J 


,  d2A  dB 

A~w  +  (l  +  Ca)tsn 


0, 


v“ 


2 vpis  i  p2 

n  l  o 


S  J 


„  d2B  ..  dA 

B  ~  C22lS  +  {  n)'Sdy  “ 


(2.19) 

(2.20) 
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The  solution  of  these  equations  which  vanishes  for  y  — >  oo  is 


A{s,y,p)  -  A\{s,p)e  7l3/  +  A2(s,p)e 

B(s,y,p )  =  — —A1(s,p)e~'lliy  -  —  A2(s,p)e~'r2y, 


(2.21) 


where  A\  and  A2  are  arbitrary  functions  and  otj(s,p)  stands  for  the  functions 


aj(s,p)  = 


(1  +  c12)jj 


,  j  =  1,2 


(2.22) 


with  71  and  7I  being  two  distinct  roots  of  the  quadratic  equation 


C2274  +|  ^2  _  C22  ^CH  —  ^2^  +  c12  +  2ci2 


s2  —  (1  +  C22)-^  —  (1  +  C22) 


2upis ) 


v2  \  2  ,  2vPis  ,  F 


Cu - ,  I  s  + 


r2  r2 


2  2 

C5  CS 

v2\  2  2vpis  p: 


1 


72  + 


1  _  _  + 


+  “7 


c2  c2 

us  J 


=  0.(2. 23) 


It  can  be  shown  that  for  many  materials  the  roots  71  and  72  are  real  and  positive  and  the 
expressions  for  the  displacements  in  the  Laplace  transform  domain  become: 


J-  [°°  (. Aie +  A2e~™)e~isx  ds, 

271  7 — 00 

(2.24) 

_ q  POO  p  isx 

„  /  {axAxe-™  +  <x2A2e-™)  ds, 

2n  J- 00  5 

(2.25) 

and  using  (2.10)  the  corresponding  expression  for  r*y  is  given  by 

r*xy  =  /_J(ai  +  7iMie-™  +  (o2  +  72 )A2e-™]e~isx  ds.  (2.26) 

Applying  the  second  condition  of  (2.15)  to  equation  (2.26)  yields 


A2(s,p) 

A 


~PiAi(s,p), 
<*1  +7i 
0:2  +  72  ' 


(2.27) 


Therefore  the  expressions  for  the  transformed  components  of  displacement  become 


<(z,y,p)  = 

i  r  _  /jle-^y)A,(s,p)e-i”  ds, 

Z7T  7—oo 

(2.28) 

u*y{x,y,p)  = 

^  f°  (o^e"™  ds, 

27T  7—00  5 

(2.29) 
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and  the  associated  stress  components  are  given  by 

K  =  ^  r[(cllS2-ai7lCl2)e-1*»-(c11s2-a272c12)fte-^]^l^e-‘“<«.30) 

Z7T  J-oc  S 

a !  =  ^  l(cas2  -  aiKCv)e-™  -  (cns2  - 
y  Z7T  J-oo  S 

r*s  =  (2-32) 


Introducing  the  functions 


E{s,p) 
F(s,p ) 


-(oi  -/5ia2)^i(s,p), 

5 

r  2 


(«1  -  PiCX2)^ 


[Ci2S2  -  Q!i7iC22  -  A(Ci2S2  -  C*272C22)], 


(2.33) 

(2.34) 


i  -  £)» » + ^ + K 


(2.35) 


(cn  -  u2/c2)(l  +  c12)(iV1  +  AT2)(1  -  r>2/4) {(C?2  +  °12  CuC22  + 

(cMN2  -  cn  +  v2/c2)  -  c22[cuN2N2  +  (cu  -  v2/c2)(N2  +  N,N2  +  Nflfc 6) 


; — {-v2/c2  +  c22(cn  -  v2/c2)  -  c\2  -  2 c12  ±  [{v2/c2s 
C22 

-c22(cn  -  v2/c2 )  +  c2 2  +  2ci2)2  -  4c22(cii  -  v2/c2)(  1  -  v2/c2)]1/2},  (2.37) 


where  the  velocity  q  =  ^/c77cs  represents  the  dilatational  wave  speed  along  the  x— axis  and 
in  view  of  the  first  and  third  boundary  conditions  in  (2.15),  equation  (2.29)  and  (2.31)  yield 
the  following  pair  of  dual  integral  equations  for  the  determination  of  the  function  E(s,p) 


°l(x,0,p) 


iput  f°°  F(s,p ) 


E(s,p)e~isx  ds  = 


oo  <  x  <  0,  (2.38) 


— i  r°° 

u*y(x,0,p)  =  —  J  E(s,p)e~is 


0  0  <  x  <  oo. 


(2.39) 


Let  u*_(x,p)  be  the  unknown  Laplace  transform  of  the  vertical  displacement  on  the 
negative  x— axis,  and  cr*+(x,p)  be  the  unknown  Laplace  transform  of  the  normal  stress  on 
the  positive  x— axis,  so  that 


u*y{x,0,p)  = 


0  for  x  >  0 

u*_(x,p )  for  x  <  0 


o-y(z,0,p) 


a^_(x,p)  for  x  >  0 
—Oq/p  for  x  <  0. 


(2.40) 
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Then  the  Laplace  transform  of  the  normal  stress  and  the  vertical  displacement  on  the  whole 
boundary  y  =  0  is  given  by 
*^12^  r°° 


2*  J-x  E(s’p)e~‘“  ds  =  “ (-*)  +  <(*. p). 


/: 


- %  ro o 

2tt  j —oo 


/CO 

jB(s,p)e_lsI  ds  =  u*  (x,p), 

-OO 

and  by  Fourier  transform  inversion,  these  equations  give 


(2.41) 

(2.42) 


where 


Vd 

=  -—  +  E+(«), 
pis 

(2.43) 

-iE(s,p) 

=  v.(«), 

(2.44) 

E+M  = 

poo 

/  a*+  (x ,  p)elsx  dx , 

(2.45) 

V-(s)  = 

J  ^u*y_(x,p)elsxdx, 

(2.46) 

L{  ;°W  = 

^0 
ips * 

From  the  physics  of  the  problem  it  is  reasonable  to  assume  that  the  function  o*+(x,p) 
and  u*_(x,p )  are  exponentially  bounded  at  infinity  and  this  ensures  the  existence  of  their 
Fourier  transform  (2.45)  and  (2.46).  In  particular  it  is  shown  in  [10]  that 

if  |cr+(a;,p)|  <  M]_ex~x  as  x  — >  -t-oo  then  £+(s)  is  analytic  in  Im(s)  =  A  >  A_, 

and  if  |u*_(x,p)|  <  M2ex+X  as  x  — >  — oo  then  V_(s)  is  analytic  in  Im(s)  =  A  <  A+. 

2.3.1  Wiener-Hopf  technique 

Eliminating  E(s,p )  from  (2.43)  and  (2.44)  yields  a  Wiener-Hopf  equation 

(2.47) 


^  -  E+w  = 

ips  0d 


which  contains  only  the  two  unknown  functions  £+(s)  and  V_(s),  and  now  the  Wiener-Hopf 
technique  can  be  applied  as  follows.  Suppose  that  the  function  L(s )  is  defined  and  factored 
as 

T.  (c\  T?( o  ■n'l 

(2.48) 


i+W 


9d 
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then  equation  (2.47)  becomes 

-  E+(*)M»)  =  £-(»)*-(»)• 

ips 

Assume  that  the  function  D(s)  can  be  defined  and  decomposed  as 


^-L+{s)  =  D(s)  =  D+(s)  +  D.(s), 
ips 


then  equation  (2.49)  becomes 


(2.49) 


(2.50) 


D+(s)  -  S+(s)L+(s)  =  L.(s)V4s)  -  D-(s)  =  W(s). 


(2.51) 


The  first  member  of  this  equation  is  analytic  in  the  upper  half  plane  Im(s)  =  A  >  A_  and  the 
second  member  in  the  lower  half  plane  Im(s)  =  A  <  A+.  Therefore,  if  A+  >  A_  the  regions 
of  analyticity  overlap.  Using  the  Liouville’s  theorem  to  determine  W (s),  solutions  for  £+(s) 
and  V_(s)  can  be  found. 

The  only  zeros  of  the  function  F(s,p )  are  of  the  form  s  =  ip/{±CR  +  v )  where  cr  is  the 
Rayleigh  wave  speed.  This  can  be  seen  by  substituting  s  =  ip/(c  +  v )  in  F(s,p),  letting 
F(s,p)  =  0  and  dividing  by  the  non-zero  factors,  then  F(ip/(c  +  v),p)  =0  reduces  to 


[c22  ( C11C22  -  C? 2  _  C^A 

Cu  V  °22  C2s)  \ 


C 2  C2 


1--^- 


1  - 


cnci 


=  0 


which  is  the  Rayleigh  function  for  orthotropic  materials  [13].  The  roots  of  this  function  are 
c  =  ±cR. 

Consequently,  the  first  step  in  factoring  L(s )  is  to  define 

F(s,p) 


F«  = 


e 


R 


where 


e 


R 


1 - 2  I  s  +  2 

V  CR/  CR 


2  vpis  pi 


+  ^2- 
CR 


(2.52) 


(2.53) 


It  can  be  shown  that  F(s)  — >  1  as  s  — >  00,  (the  function  £(w)  in  (2.34)  was  chosen  to  make 
this  possible).  The  function  F(s)  is  regular  and  F(s)  7^  0  in  the  s— plane  cut  as  shown  in 
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Im(s) 


\  i  p  /  (cs+v) 

\ 

i  p  /  (cd+v)\ 


Re(s) 


p  /  (Cd-V) 

\ 

-i  p  /  (cs-v)\ 


Figure  2.2:  Branch  cuts  of  F(s)  in  the  s— plane. 

figure  2.2,  the  only  singularities  are  the  branch  points  shared  with  71  and  72.  Where  the 
branch  points  of  71  and  72  are 


for 

7i5 

IP 

—ip 

&  1  J 

Cs  +v 

cs~v 

for 

72; 

ip 

—ip 

Cd  +  v' 

cd  —  v 

It  is  well-known  that  factorization  is  accomplished  most  directly  for  functions  that  ap¬ 
proach  unity  as  |s|  — >  00  and  that  have  neither  zeros  nor  poles  in  the  finite  plane;  F(s)  is 
an  example  of  such  a  function.  Therefore,  using  Cauchy’s  integral  formula  it  can  be  shown 
that  [5] 

logF(^) 


F±(s) 


dz 


lr±  Z  —  S  ) 

where  F(s)  =  F+(s)F-(s)  and  T_  [T+]  is  the  contour  enclosing  the  cut  between  the  branch 
points  +ip/(cs  +  v )  and  +ip/ (c<*  +  v ),  [—ip/(cs  —  v )  and  —ip/(cd  —  «)].  Using  the  fact  that 


F(s)  =  F(s)  one  can  write 

F±(s)  =  exp  ^  J  1 3=F 

V  J\/cd^v 

From  equations  (2.35)  and  (2.53)  note  that 
Qd  = 


1^Cs=FUtan“1  die 

*  Ji/c&v  \Re[F(ipw)}J  wT  \ 


N 


1--2  b2  + 


2  vpis  p 2 


Or  — 


v2  \  9  2 vpis  p2 

1  -  —  s2  +  -f  1  v 


-Rj 


-R 
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Returning  to  the  factorization  of  L(s )  we  have 

M«)  „  ,A»)9ii 

i(s)  " 


~  P 12^ 


(2.54) 

(2.55) 


therefore 


L-(s)  = 

.  *>  ■ 


and 


with 


D(s)  =  g-L+{s)  =  ^ 
ips  p 


Si  F+(») 


M»)  -  i+(0)  +  i+(0) 


25 


25 


M») 


P 


M»)  -  Mo) 


25 


£>-(*)  = 


=  D+(s)  +  D_(s) 

"MO)' 


Oo 

P 


IS 


(2.56) 

(2.57) 

(2.58) 

(2.59) 


Each  side  of  equation  (2.51)  is  analytic  in  one  of  the  overlapping  half  planes,  and  the  sides 
coincide  on  the  strip  of  overlap.  Consequently  each  side  of  (2.51)  is  the  analytic  continuation 
of  the  other  into  its  complementary  half  plane;  so  that  the  two  sides  together  represent  one 
and  the  same  entire  function  W (s).  The  entire  function  will  be  determined  by  its  behavior 
at  |s|  — >  oo  which  is  related  with  the  behavior  of  physical  quantities  near  x  =  0.  First 
note  that  L+(s )  ~  s-1/2  and  L_(s)  ~  s1/2  as  |s|  — >  oo,  and,  D+(s )  and  D_(s)  in  (2.59) 
are  bounded  in  their  respective  planes  of  analyticity  and  vanish  at  infinity.  Furthermore, 
cr+(x,p)  is  expected  to  be  square  root  singular  as  x  — >  0+  and  u*_(x,p)  is  expected  to  vanish 
as  x  — »  0_  to  ensure  continuity  of  displacement.  As  a  result,  from  the  Abel  theorem  [10] 
relating  asymptotic  properties  of  transforms,  we  know 

lim  x1^2al(x,p)  ~  lim  s1|/2E+(s), 

a;-K)+  +  s-+  oo  +v  n 

lim_  |a:|_9u*  (x,p)  ~  lim  |s|1+9Vr_(s), 


£->>0 


s— »— oo 


for  some  q  >  0.  Therefore,  it  is  expected  that  E+(s)  ~  s  x/2  and  V-(s)  ~  s  1  9  as  |s|  — >  oo, 
thus  the  products  E+(s)L+(s)  and  L_(s)V_(s)  vanish  at  infinity.  Therefore,  each  side  of 
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(2.51)  vanishes  as  |s|  — »  oo  in  the  corresponding  half  planes.  According  the  the  Liouville’s 
theorem,  a  bounded  entire  function  is  constant.  In  this  case,  W (s)  is  bounded  in  the  finite 
plane  and  W(s)  — >  0  as  |s|  — >  oo  so  that  the  constant  must  be  zero;  thus,  kF(s)  =  0.  By 
using  (2.51)  and  (2.59),  the  functions  of  interest  are  then  given  by 


E+(s) 

V-(s) 


D+(s)  _  _goJ_  ’MO) 
L+(s)  p  is  L+(s) 
D-(s)  _  go  1  M°) 
L_(s)  p  is  L_(s)' 


(2.60) 

(2.61) 


2.3.2  Stress  Intensity  factor 


To  find  the  stress  intensity  factor,  an  asymptotic  expression  for  the  normal  stress  near  the 
crack  tip  is  sought.  The  Abel’s  theorem  relating  asymptotic  expressions  between  a  function 
and  its  Fourier  transform  is  the  following  [10] 

lim  y/x  al(x,p)  =  lim  E+(s).  (2.62) 

x-»0+  +  s-H-oo  V  7T 


Clearly,  the  behavior  of  E+(s)  as  s  — t  oo  is  needed.  Note  that 

yjl  -  v/cd  i 

L+(S)  =  ml-v/cR7^  aS  S  °° 

and  since  F+( 0)  =  FL(0)  =  \/P( 0)  =  ^^/F(0,p)  then 


MO) 


P 

V*PCd\fF(  0,P) 


therefore 


M«)  =  77 


-<70  1  1  -  V/CR 


P 


as  s  — y  oo. 


(ip)z/2  s1'2  y/l-v/a  y/ciy/F( 0,p) 

Using  this  relation  and  the  definition  of  the  stress  intensity  factor  in  (2.62)  gives 


(2.63) 


Kt(p)  =  lim  \/2nx  al(x,p)  =  lim  e  ™l4\p2s  E+(s) 

x->0+  '  s— t+oo 

__  <7q\/2  1  -  c/cr _ p 

PV 2  y/l-c/dy/ciy/F^p) 


(2.64) 

(2.65) 
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it  is  readily  shown  that 


V^(0,P) - 


so  that 


K*j(p)=a0 


P\ 


1  -  c/cfi 


2cs£  1 


(2.66) 


(2.67) 


Vw^1^p3/2 

and  by  Laplace  inversion  the  dynamic  stress  intensity  factor  in  the  time  domain  for  this 
loading  mode  is 

1  -  V/CR 


Ki(t,v )  =  2(j0- 


2cs^(v) 


Vi. 


-  V/Vd\  *V°£ 

Denoting  £(0)  the  function  f  evaluated  at  v  =  0,  equation  (2.68)  can  be  rearranged 


(2.68) 


Ki(t,v ) 


1  —  v/c 


R 


«WU?SiVi 


(2.69) 


The  expression  in  brackets  corresponds  to  the  dynamic  stress  intensity  factor  for  a  stationary 
semi-infinite  crack  in  orthotropic  materials  derived  in  [12],  and  the  remaining  factor  will  be 
k[(v),  the  universal  function  of  the  crack  tip  speed  that  relates  the  stress  intensity  factor  for 
stationary  and  propagating  crack.  That  is 

1-  v/cR 


kT(v)  = 


W 

£(0 )’ 


Kj&v)  =  kj^Kt&O) 


The  expression  for  kj(v)  can  be  approximated,  for  0  <  v  <  cr,  by 

1  -  v/cr 


(2.70) 


ki(v) 


yfl  -  ~v~/ ~cd 


(2.71) 


as  in  the  case  of  isotropic  materials. 


2.4  In-Plane  Shear  Loading 

Consider  the  crack  geometry  illustrated  in  figure  2.1(b).  The  crack  faces  are  subjected  to 
suddenly  applied,  spatially  uniform  shear  traction  of  magnitude  r0  at  time  t  =  0.  Exploiting 
asymmetry  and  examining  the  half  space  y  >  0,  the  corresponding  boundary  conditions  are 

rxy(x,0,t)  =  —T0H(t)  for  —  oo  <  x  <  0, 
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cry(x, 0,  t)  =  0  for  -  oo  <  x  <  oo, 

ux(x,0,  t)  =  0  for  x  >  0, 


(2.72) 


where  H(t)  is  the  Heaviside  step  function.  In  addition  the  condition  of  zero  displacements 
at  infinity  and  zero  initial  conditions  are  assumed. 

The  method  of  solution  is  similar  to  that  used  for  normal  impact.  Applying  Laplace 
transform  to  boundary  conditions  (2.72)  gives 

T*y(x,0,p)  =  -T0/p  for  -  oo  <  x  <  0, 

a*(x,0,p)  =  0  for  -  oo  <  x  <  oo,  (2.73) 

u*x(x,0,p)  =  0  for  x  >  0. 

Assuming  the  displacement  field  (2.17)  and  (2.18)  yields  the  same  system  of  ordinary  differ¬ 
ential  equations  (2.19)  and  (2.20)  whose  solution  is  given  by  equations  (2.21)  where  a tj  are 
defined  by  (2.22)  and  7 j  are  obtained  from  the  solution  of  (2.23).  Using  (2.24),  (2.25)  and 
(2.10)  indicates  that  a*  is  given  by 

oil  roo  g—isx 

a*(x,y,p)  =  — [  [(s2c12-aiiic22)A1e~'riy  +  (s2c12-a2l2C22)A.2e~'r2y]——ds.  (2.74) 

y  Z'K  J —00  s 

Applying  the  second  condition  of  (2.73)  to  equation  (2.74)  yields 

A2  (s,p)  =  -foAi(s,p), 


S2Ci2  -  G!i7iC22 

P2  9  * 

(2.75) 

"  C12  —  <*272^22 

Therefore,  the  expressions  for  the  transformed  components  of  displacement  become 

K(x,y,p)  =  if 

Z7T  J— 00 

(2.76) 

u‘v(x,y,p)  =  ^  r  (a,e~™  -  ds, 

*  Z'K  00  S 

(2.77) 

and  the  associated  stress  components  are  given  by 

= 

J1  [  [(cnS2  ai7iC12)e  liy  (cus2  72y]  e  d&. 78) 

Z'K  J—  00  S 

<  = 

f 12  I"  (cus2  -  aaiC22)(e~™  -  e-™)^^-e-isxds, 

Z'K  J- oo  S 

(2.79) 

<y  = 

7“2  r  l(“i  +  7i)e-™  -  (“2  +  72)A e-™\At(s,p)e-“*ds. 

Z'K  J —00 

(2.80) 
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Introducing  the  functions 


(2.81) 

(2.82) 


C(s,p) 

G(s,p) 


0s{v) 


(1  -p2)Ai(s,p), 
-0, 


(1  -  (32)r) 


[(«l+7l)  -ft(G!2  +  72)]> 


1-4  W2  + 


2  vpis  p2 


+ 


N 

Cn^/l  -  u2/c^l  -  u2/q 
C22N1N2 


Z(v) 


(2.83) 

(2.84) 


where  £  and  AT1j2  are  given  in  (2.36)  and  (2.37)  and  in  view  of  the  first  and  third  boundary 
conditions  in  (2.73),  equation  (2.76)  and  (2.80)  yield  the  following  pair  of  dual  integral 
equations  for  the  determination  of  the  function  C(s,p), 


«*(®»0,p)  =  ^:^ooC'(s,p)e"isa:  ds 


— -  —  00  <  x  <  0,  (2.85) 

P 

0  0  <  x  <  00.  (2.86) 


Let  u*_(x,p)  be  the  unknown  Laplace  transform  of  the  displacement  on  the  negative 
x— axis,  and  r|(x,p)  be  the  unknown  Laplace  transform  of  the  shear  stress  on  the  positive 
x— axis,  so  that 


<(x,0, 


0 

<-(• x,p) 


for 

for 


x  >  0 
x  <  0 


r*y(x,o,p) 


r+(x,p)  for  x>0 
—tq/p  for  x  <  0. 


(2.87) 


Then,  the  Laplace  transform  of  the  shear  stress  and  the  displacement  on  the  whole  boundary 
y  —  0  is 


^  l-ooC^S,P^e~tSX  dS  =  U*x~^ 


and  by  Fourier  transform  inversion,  these  equations  give 

»nV^f-C(s,P)  =  -^  +  T+(S) 
C(s,p )  =  U-(s) 


(2.88) 

(2.89) 


4'5 


(2.90) 

(2.91) 


where 


roo 

T+{s)  =  /  r*+{x)elsxdx 

J  o 

U-(s)  =  f°  u*x_(x)eisxdx 
J— oo 


(2.92) 

(2.93) 


»'sxdx  = 

ips 


Eliminating  C(s,p )  from-  (2.90)  and  (2.91)  results  in  a  Wiener-Hopf  equation 


To  ^  ..  „G(s>P)rr  ^ 

t—  -  T+(s)  =  -HuV—fi - <MS) 

ZpS  U  § 


(2.94) 


which  contains  only  the  two  unknown  functions,  T+(s )  and  £/_(s),  and  can  be  solved  using 
the  Wiener-Hopf  technique  as  in  the  normal  impact  analysis.  The  result  is 


D+(s)  -  T+(s)L+(s)  =  L-(s)U-(s)  -  D_(s)  =  0 


(2.95) 


where 


L-(s)  =  -/ii  2^G_(»), 


and 


with 


where 


and 


£>(s)  =  ^-M«)  =  ^ 

zps  p 


M<)  -  MO)  +  MO) 


zs 


zs 


=  T>+(s)  +  £>_(s), 


(2.96) 

(2.97) 

(2.98) 


D+(s)  = 


70 

P 


L+(s)  -  T+(0) 

,  and  D-(s )  =  — 

[i+(0)l 

is 

P 

is 

G(s)  = 


G&p) 
Or  ’ 


(2.99) 


e,= 


\ 


1-4  »>+ 


2upis  p2 


+ 


<?s 
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It  can  be  shown  that  G(s )  — >•  1  as  s  — >■  oo.  Thus,  the  solution  of  (2.95)  is 

T  (> \  _  g+(g)  =  ro  1  [Mg)  , 

+V' }  L+{s)  p«[L+(s) 

tj  (o\  =  D-(S)  =  r0  1  M°) 

T_(s)  p  is  L-(s)' 

with  L±(s)  defined  in  (2.96)  and  (2.97)  and  G±(s)  given  by 

GUs)  =  exb  !  —  f1/CsTV  tan-1  fIm[G(^)]\ 

|  7r  Ji/cd^v  \Re[G(ipu;)]/  w  =F  ^ 

To  calculate  the  stress  intensity  factor,  T+(s )  as  s  -»  oo  is  needed.  First,  note  that 

Jl  —  v/cs  i 

L+to  =  T-i/cR7&  38  s^°° 

and  since  G+(0)  =  GL(0)  =  ^(7(0)  =  ^\[g{ 0,p)  then 


(2.100) 

(2.101) 


L+(0)  = 


y/vpCs\fG{^p) 


Therefore, 


/  x  _  -To  1  1-v/cr 

\S)  /•  \a/o  -i/o  /  : 


P 


(*p)3/2  *1/2  ^/l-v/c,  ^Fsy/G{oJ) 
Using  (2.62)  and  the  definition  of  the  stress  intensity  factor  gives 


as  s  — y  oo. 


(2.102) 


K*n{p)  —  fim  V2irx  rl(x,p)  =  lim  e  m^y/2 s  T+(s ) 

X— ►0+  S~>  +  CXD 

_  T0y/2  l-v/cR  p 


P3/2  ^l-v/cs  ^G(0,p)' 


It  can  be  shown  that 


so  that 


Kh(p)  = 


1  1  -  v/cR 


(2.103) 

(2.104) 

(2.105) 

(2.106) 


and  by  Laplace  inversion  the  dynamic  stress  intensity  factor  in  the  time  domain  for  this 
loading  mode  is 

Kn{t ,  v)  =  2 r0  ^/2c577(u)/7t  Vt.  (2.107) 

V1  “  V/°s 
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f 

Graphite 

Epoxy 

E-Glass 

Epoxy 

Boron 

Epoxy 

Epoxy 

cn 

20.77 

8.38 

32.67 

3.91 

C22 

2.18 

2.29 

3.12 

3.91 

C12 

0.52 

0.79 

7.48 

5.5 

6.4 

fTOff 

p  (Kg/m3) 

1600 

2100 

1990 

1260 

Table  2.1:  Mechanical  properties  used  for  the  analysis. 


Denoting  77(0)  the  function  77  evaluated  at  v  =  0,  equation  (2.107)  can  be  rearranged  as 

1  -  v/cR 


Kn(t,  v )  = 


T](v) 


2ro^/2cs77(0)/7r  \fi t 


(2.108) 


yfl  -  vjcs  \  *7(°)  L 

The  expression  in  brackets  corresponds  to  the  stress  intensity  factor  for  a  stationary  semi¬ 
infinite  crack  in  orthotropic  materials  under  impact  shear  loading  derived  in  [12],  and  the 
remaining  factor  will  be  kji(v),  the  universal  function  of  the  crack  tip  speed  that  relates  the 
stress  intensity  factor  for  stationary  and  propagating  crack.  That  is 

1  -  v/cR 


kn(v) 


V(v) 


Kn(t,v)  =  kn(v)Kn(t,  0) 


yjl^vjca  \  V(0) 

The  expression  for  k//(v)  can  be  approximated,  for  0  <  v  <  cR,  by 

1  -  v/cR 


(2.109) 


kn(v) 


71  ~  V/Cs 


(2.110) 


as  in  the  case  of  isotropic  materials. 


2.5  Results  and  Conclusions 


Closed  form  expressions  for  the  dynamic  stress  intensity  factors  Ki(t,v)  and  Kn(t,v )  due 
to  suddenly  applied  uniform  loads  on  the  crack  faces  have  been  determined  for  semi-infinite 
cracks  propagating  at  constant  velocity  v  in  orthotropic  materials.  The  results  presented  here 
are  consistent  with  those  for  stationary  cracks,  i.e.  for  v  =  0,  Ki(t,  0)  and  Ku(t,  0)  developed 
in  a  previous  work  by  Rubio-Gonzalez  and  Mason  [12]  for  semi-infinite  cracks  in  orthotropic 
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materials  subjected  to  spatially  uniform  step  loads  on  the  crack  faces  are  obtained.  Using 
these  results  it  has  been  possible  to  find  expressions  for  the  universal  function  of  the  crack 
tip  speed  kj(v )  and  ku(v)  which  may  be  written  in  approximated  manner  as 


ki(v)  ~ 
kn(v)  » 


1  ~  v/cr 

\fl^vjcd 
1  ~  v/cr 
yjl-v/c. 


Note  that  these  equations  conserve  the  same  form  as  those  for  isotropic  materials  [5],  where 
the  wave  speeds  cs  and  c<*  should  be  computed  and  understood  properly  for  orthotropic 
materials  with  the  planes  of  orthotropy  assumed  here. 

Figures  2.3  and  2.5  illustrate  kj(v )  and  kn(v)  for  different  materials  for  which  the 
properties  are  given  in  Table  2.1.  They  also  show  that  for  orthotropic  materials  the  functions 
ki(v)  and  kn(v)  given  by  equations  (2.70)  and  (2.109)  may  be  approximated  by  the  simplified 
expressions  (2.71)  and  (2.110),  respectively  as  usually  it  is  for  isotropic  materials.  They  are 
very  close  each  other  for  several  materials.  In  fact,  the  factor  y £(w)/£(0)  is  close  to  one  over 
a  large  interval  of  speeds  as  shown  by  figure  2.4.  It  can  be  shown  that  £(w)/£(0)  =  7](v)/r](0) 
and  therefore  Figure  2.4  applies  to  the  in-plane  shear  loading  too. 

It  was  assumed  tacitly  in  the  method  of  solution  of  the  equations  of  motion  that  the  crack 
speed  propagation  is  subsonic,  v  <  cr,  in  order  to  maintain  elliptic  differential  equations. 
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Universal  function  k^v),  Epoxy 


Universal  function  k,(v),  Graphite-Epoxy  x 


Figure  2.3:  Comparison  of  the  universal  function  ki{v)  and  its  approximated  expression 
(2.71)  for  different  materials  in  mode  I  loading. 


Factor  ( 5(v)  /  5(0)  )1/2  for  different  materials 


Figure  2.4:  Factor  y£(v)/t;(0)  for  different  materials.  Graphite-epoxy  x  or  y  means  that  the 
fibers  are  parallel  to  the  x—  or  y— axis,  respectively. 
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k„(v) 


Universal  function  k,(v),  Eglass-epoxy 


Universal  function  k..(v),  Boron-epoxy 


Figure  2.5:  Comparison  of  the  universal  function  ku(v )  and  its  approximated  expression 
(2.110)  for  different  materials  in  mode  II  loading. 
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Chapter  3 


Dynamic  Stress  Intensity  Factor  due 
to  Concentrated  Normal  Loads  on 
Semi-infinite  Cracks  in  Orthotropic 
Materials. 


Co-authored  with  C.  Rubio- Gonzalez  and  submitted  to  the  Journal  of  Composite  Materials 
The  transient  elastodynamic  response  due  to  concentrated  normal  impact  load  on  the 
faces  of  a  a  semi-infinite  crack  in  an  orthotropic  material  is  examined.  In  contrast  to  earlier 
papers  where  numerical  approximations  were  used,  a  closed  form  solution  for  the  stress 
intensity  factor  history  around  the  crack  tip  is  found  here.  Laplace  and  Fourier  transforms 
together  with  the  Wiener-Hopf  technique  are  employed  to  solve  the  equations  of  motion  in 
terms  of  displacements.  Even  though  the  problem  has  characteristic  length,  it  has  been  shown 
in  previous  works  that  the  Wiener-Hopf  technique  can  be  applied.  The  asymptotic  expression 
for  the  stress  near  the  crack  tip  is  analyzed  which  leads  to  the  dynamic  stress  intensity  factor 
in  mode  I.  Similarly  to  the  isotropic  case,  it  is  found  that  the  stress  intensity  factor  has  a 
singularity  and  discontinuity  when  the  Rayleigh  wave  emitted  from  the  load  arrives  at  the 
crack  tip.  Results  are  presented  for  orthotropic  materials  as  well  as  for  the  isotropic  materials. 
The  closed  form  solution  is  given  by  simple  integral  and  algebraic  expressions  and  does  not 
exhibit  the  spurious  oscillations  seen  in  earlier  numerical  solutions. 

Keywords:  stress  intensity  factor,  dynamic  fracture,  orthotropic  materials. 
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3.1  Introduction 


The  growing  use  of  composites  in  many  engineering  applications  demands  the  fundamen¬ 
tal  understanding  of  the  response  of  cracked  orthotropic  bodies  under  impact  loads.  The 
behavior  of  finite  cracks  with  concentrated  loads  on  its  faces  in  orthotropic  materials  has 
been  analyzed  by  Rubio-Gonzalez  and  Mason  [8,  9]  using  integral  transform  methods.  In 
that  work,  the  problem  is  reduced  to  a  Fredholm  integral  equation  in  the  Laplace  transform 
domain  which  is  solved  numerically,  and  the  stress  intensity  factor  is  recovered  in  the  time 
domain  by  numerical  Laplace  inversion.  Although  this  approach  can  be  quite  accurate,  the 
solution  obtained  is  approximated  and  restricted  to  finite  cracks.  It  would  be  desirable  to 
have  an  exact  solution  for  the  stress  intensity  factor  in  cracks  in  orthotropic  materials  sub¬ 
jected  to  impact  concentrated  loads  that  can  be  used  as  a  Green’s  function  in  dynamic  crack 
problems  in  orthotropic  materials. 

In  the  present  work  a  semi-infinite  crack  with  impact  concentrated  loads  on  its  faces 
in  orthotropic  material  is  analyzed,  an  exact,  closed-form  solution  for  the  dynamic  stress 
intensity  factor  is  obtained  for  the  opening  loading  mode.  Laplace  and  Fourier  transforms 
are  used  along  with  the  Wiener-Hopf  technique  to  find  the  stress  ahead  of  the  crack  tip  and 
the  displacements  of  the  crack  faces.  Asymptotic  expression  for  the  stress  near  the  crack 
tip  leads  to  the  stress  intensity  factor  Ki(t).  The  equivalent  problem  for  isotropic  materials 
was  solved  by  Freund  [2]  by  superposition  of  two  different  problems  with  no  characteristic 
fixed  length.  The  superposed  problems  were  solved  employing  the  Wiener-Hopf  technique. 
Kuo  and  Chen  [4,  5]  have  shown  that  the  direct  application  of  the  Wiener-Hopf  technique 
may  be  used  to  solve  that  problem  in  a  more  straightforward  manner.  The  latter  approach 
is  used  in  this  paper  to  extend  the  solution  of  Freund  to  orthotropic  materials. 
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3.2  Governing  Equations 


Consider  the  plane  problem  of  an  infinite  orthotropic  medium  containing  a  semi-infinite 
crack,  Figure  3.1.  Let  Eiy  %  and  i ^  ( i,j  =  1,2,3)  be  the  engineering  elastic  constants  of 
the  material  where  the  indices  1,  2,  and  3  correspond  to  the  directions  (x,y,  z)  of  a  system 
of  Cartesian  coordinates  chosen  to  coincide  with  the  axes  of  material  orthotropy.  The  crack 
faces  are  along  the  negative  x— axis  and  the  origin  of  the  xy  axes  is  the  crack  tip.  The  crack 
faces  are  suddenly  loaded  by  a  pair  of  concentrated  normal  forces  of  magnitude  <r0  located  a 
distance  l  away  from  the  crack  tip,  as  shown  in  Figure  3.1. 


Figure  3.1:  Schematic  of  the  semi-infinite  crack  geometry. 


The  problem  is  restricted  to  two  dimensions  with  wave  propagation  in  the  x  —  y  plane 


only.  By  setting  all  the  derivatives  with  respect  to  z  to  be  zero,  it  is  readily  shown  that  the 


displacement  equations  of  motion  [6]  reduce  to 


d2u  d2u  .  .  d2v 

C"dx>  +  dy*+{1+Ctt)dxdy  - 

1  d2u 

C2  dt2  ’ 

(3.1) 

d2v  d2v  .  .  d2u 

a  2  +  c22  a  2  +  (1  +  c12 )  0  Q  — 

ox1  oyz  oxoy 

1  d2v 
c2  dt 2  ’ 

(3.2) 

where  u  and  v  are  the  x  and  y  components  of  the  displacement  vector  and  cn,  c\2  and  c2 2 
are  non-dimensional  parameters  related  to  the  elastic  constants  by  the  relations: 


Ex 

Cu  yx2{i-(E2/E1yl2y 

C22  =  (E2/Ex)cn,  (3.3) 

C12  =  ^12^22  =  ^2lCn, 
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for  generalized  plane  stress,  and  by 

Cn  =  £k(1  -  Va'/32)’ 

022  =  ^(1  -  "'S"31)- 

E\  (  E2  N, 

Cl2  =  ^("21  +  B^13"32*’ 

A  =  1  —  ^12^21  —  V2ZVZ2  —  —  ^12^23^31  —  ^13  ^21  ^32  > 


(3-4) 


for  plane  strain.  In  the  orthotropic  solid,  cs  =  yfpu/p  represents  the  velocity  of  the  in-plane 
shear  wave  propagating  along  the  the  principal  material  axes  and  p  is  the  mass  density.  The 


stresses  are  related  to  the  displacements  by  the  equations: 


Cf x 

du 

dv 

Hl2 

°ndx 

+  Ci2^~ , 

dy 

(Ty 

du 

dv 

If 

Pl2 

=  Ci2  — 
OX 

+  C22  77*  > 

dy 

Txy 

du 

dv 

P‘12 

dy  + 

dx 

(3.5) 


3.3  Method  of  Solution 


Exploiting  symmetry  and  taking  only  the  upper  half  plane  y  >  0,  the  corresponding  bound¬ 
ary  conditions  are 

oy{x,  0,f)  =  -a0H(t)S(x  +  l)  for  -  oo  <  x  <  0, 

Txy(x,  0,  t)  =  0  for  —  oo  <  x  <  oo,  (3.6) 

v(x,0,t)  =  0  for  x  >  0, 

where  H(t)  is  the  Heaviside  step  function.  In  addition,  the  condition  of  zero  displacements 
at  infinity  and  zero  initial  conditions  are  assumed. 

The  method  of  solution  of  the  governing  equations  presented  here  follows  that  described 
in  Kuo  and  Chen  [4,  5]  for  the  isotropic  case  with  some  significant  differences.  Displacement 
potentials  are  not  used  and  no  assumptions  are  made  about  the  form  of  the  unknown  func¬ 
tions.  In  equations  (3.1)  and  (3.2),  the  time  variable  may  be  removed  by  application  of  the 
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Laplace  transform 


r  oo  1  r 

/»  =  /  m  e~rt  it,  /(<)  =  —  /  f(p )  e*  dt,  (3.7) 

«/ 0  2/7TZ  J  Br 


where  Br  denotes  the  Bromwich  path  of  integration  which  is  a  line  parallel  to  the  imaginary 
axis  in  the  p-plane.  Applying  relations  (3.7)  to  equations  (3.1)  and  (3.2)  and  using  zero  initial 
conditions  for  the  displacements  and  velocities,  the  transformed  field  equations  become 


d2u*  d2u*  d2v* 

Cn a?  +  +  (  +  Cl2)e^ 

d2v*  d2v*  .  .  d2u* 

+  c22~  o'  +  (1  +  c12) 


dx 2 


dy2 


dxdy 


„2 

P  * 

-  — U 
P2 

V 


0, 

0, 


(3.8) 

(3.9) 


where  the  transformed  displacement  components,  u*  and  v*,  are  now  functions  of  the  vari¬ 
ables  x,  y,  and  p.  The  application  of  the  Laplace  transform  to  the  boundary  conditions  (3.6) 
gives 


<7y0 M,P) 

*£*(*»  °»p) 
u*(rr,0,p) 


— a0-S(x  + 1)  for  —  oo  <  x  <  0, 
P 

0  for  —  oo  <  x  <  oo, 

0  for  x>  0. 


(3.10) 


To  obtain  a  solution  of  the  differential  equations  (3.8)  and  (3.9)  subject  to  conditions 
(3.10),  the  Fourier  transform  is  applied, 

/oo  1  roo 

f(x)  elsx  dx,  f(x)  =  —  /  F(s )  e~tsx  ds.  (3.11) 

-oo  z7T  J  —  oo 

It  is  assumed  that  the  displacements  in  the  Laplace  transform  domain  have  the  form 


u*(x,y,p)  = 

roo 

/  A(s,y,p)e~lsx  ds, 

f  — OO 

(3.12) 

v*(x,y,p)  = 

roo 

/  B(s,  y,p)e~lsx  ds, 

f  — OO 

(3.13) 

where  A  and  B  are  the  Fourier  transforms  of  the  Laplace  transform  of  the  displacements, 
u*  and  v*,  respectively,  and  are  yet  to  be  determined.  Substituting  these  transforms  into 
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equations  (3.8)  and  (3.9),  the  functions  A  and  B  are  found  to  satisfy  the  simultaneous 
ordinary  differential  equations 


(cuS2+P2/^-^  +  (1  +  c12)^  =  0, 

(s2+p2/c2s)B-c22^  +  {l  +  cl2)is^  =  0. 

The  solution  of  these  equations  which  vanishes  for  y  -4  oo  is 

A(s,y,p)  =  Ai(s,p)e~'yiy  +  A2{s,p)e~'r2y, 

B(s,  y,p)  =  —A1(s,p)e-™-l^A2(s,p)e-™, 
s  s 

where  Ai  and  A2  are  arbitrary  functions  and  otj  ( s ,  p)  stands  for  the  functions 

,  ,  c„s2+p2/c?- 7j  „ 
aj(s,P)~  (i+Cl2)7.  ’  3  X'2 


(3.14) 

(3.15) 


(3.16) 


(3.17) 


with  7i  and  'y2  being  two  distinct  roots  of  the  quadratic  equation 

C2271  +  [(<4  +  2c12  -  c„C22)s2  -  (1  +  c22)p2/<472  +  (c„S2  +  p2/<4(s2  +  p2/cj)  =  0.  (3.18) 

It  can  be  shown  that  for  many  materials  the  roots  71  and  72  are  real  and  positive  and  the 
expressions  for  the  displacements  in  the  Laplace  transform  domain  become: 


v*  = 


1 

27T  J 

r  00 

/  (. A ie-™  +  A2e-72S,)e"ISI  ds, 

'-OO 

(3.19) 

—i 

2n  . 

roc  g—isx 

/  («!  A^-™  +  a2A2e~l2V) - ds, 

1 — 00  S 

(3.20) 

and  using  (3.5)  the  corresponding  expression  for  r*y  is  given  by 

T*  =  f°°  [(<*!  +  7i)Aie-™  +  (c*2  +  72 )A2e-™]e~isx  ds. 

"y  L1X  J —00 

Applying  the  second  condition  of  (3.10)  to  equation  (3.21)  yields 


(3.21) 


A2(s,  p)  =  -Ai4i(s,p), 
„  <*1  +  7i 


Pi  = 


a2  +  72 


(3.22) 
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Therefore,  the  expressions  for  the  transformed  components  of  displacement  become 

1  f°° 


1 

u’(x,y,p)  =  tt 

Z7T  J— oo 

«•(: r,y,p)  =  ^  f°  (<*i e"™  -  Aa2e-™)'4l(s,p)--i 

Z7T  J —oo 


ds, 


and  the  associated  stress  components  are  given  by 

-*Pl2  roo 


<Ld  = 


°l  = 


27T 

~^12 
27T 
— Ml2 


[  [(cns2  -  ai7iCi2)e  7iy  -  (cns2  -  a272c12)/?ie  72y]  ~  ’ ~ e 
7-00  '  ^ 


(3.23) 

(3.24) 

<$.25) 


-OO 

roo 


/oo 

[(ci2s2  -  ai7iC22)e_7iy  -  (ci2s2  -  a272C22)/?ie~72y] 

-OO 

/OO 

(ai  +7i)[e"7iy  -  e_72y] (s,  p)e~tsxds. 

-OO 


Iy  2tt 

The  following  functions  are  introduced, 


£(s,p)  =  —  (<*i  -/31o;2)Ai(s,p), 
s 


4i(giP)e-fa«^26) 


(3.27) 


F(s,p)  =  _  ^iC^[Ci2S2  -  Q!l7lC22  -  A(Ci2S2  -  <*27^)], 

£  =  ~  7T  ,  ~  wAr  ,  at T { (^12  +  c12  -  CnC22)(Ci2^liV2  -  Cn) 


(3.28) 

(3.29) 

(3.30) 


TV2 

ivl,2 


Cn(l  +  ci2)(Ni  +  N2) 

— c22[ci2Ar2Ar|  +  cn(iV2  +  iVxA^  +  iV| )]}, 

- — {cnc22  —  c22  —  2C12  dr  [(C11C22  —  cf2  —  2ci2)2  —  4cnC22]1^2},  (3.31) 

^C22 


where  the  velocity  c<*  =  ^/c^Cs  represents  the  dilatational  wave  speed  along  the  x— axis.  In 
view  of  the  first  and  third  boundary  conditions  in  (3.10),  equation  (3.24)  and  (3.26)  now  yield 
the  following  pair  of  dual  integral  equations  for  the  determination  of  the  function  E(s,p ) 

E(s,p)e~isx  ds  =  -^6(x  +  l)  —  oo  <  x  (3.3£) 


*/  n  n  *^124  r* 


p 


^s2+p2/^ 

_ j  roc 

v*(x,0,p)  =  —  E{s,p)e~lsx  ds  =  0  0  <  re  <  oo,  (3.33) 

27 r  J— oo 

Let  v*_(x,p)  be  the  unknown  Laplace  transform  of  the  vertical  displacement  on  the 
negative  x— axis,  and  cr+(x,p)  be  the  unknown  Laplace  transform  of  the  normal  stress  on 
the  positive  x— axis,  so  that 


v*(x,0,p) 


0 


for  x  >  0 


v*_ (x,p)  for  x  <  0 


e*Jx 


p)  for  x  >  0 

(a0/p)6(x  + 1 )  for  x  <  0 

(3.34) 
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then  we  can  write  the  Laplace  transform  of  the  normal  stress  and  the  vertical  displacement 

t 

on  the  whole  boundary  y  —  0  as 


ijM2^ 

27T 


/: 


T==  £fap)e-‘“  ds  =  -^S(x  +  l)+a'+(x,p), 
yjs2+p2/c2d  P 


' d 

-i  r°° 


/oo 

E(s,p)e~tsx  ds  =  v*_(x,p), 

-OO 


2?r  j  — oo 

and  by  Fourier  transform  inversion,  these  equations  give 

•Pu  S-r===E(s,p)  =  ——e~isl  +  £+(s), 


^s2+p2/c 


p 


—iE(s,p)  =  V.(s), 


(3.35) 

(3.36) 


(3.37) 

(3.38) 


where 

TOO 

£+(s)  =  /  o*+{x,p)elsxdx ,  (3.39) 

J  o 

V-(s)  =  f°  v*_{x,p)eisxdx,  (3.40) 

J  —  OO 

From  the  physics  of  the  problem  it  is  reasonable  to  assume  that  the  function  a*+{x,p) 
and  v*_(x,p)  are  exponentially  bounded  at  infinity  and  this  ensures  the  existence  of  their 
Fourier  transform  (3.39)  and  (3.40).  In  particular  it  is  shown  by  Noble  [7]  that 

if  |<7+(a;,p)|  <  M\ex~x  as  x  — >  +oo  then  £+(s)  is  analytic  in  Im(s)  =  A  >  A_, 

and  if  |u*  (x,p)|  <  M2ex+X  as  x  — >  — oo  then  VL(s)  is  analytic  in  Im(s)  =  A  <  A+. 


3.3.1  Wiener-Hopf  technique 


Eliminating  E(s,p )  from  (3.37)  and  (3.38)  we  obtain  a  Wiener-Hopf  equation 

2V‘*'  -  E+M  =  put  JM=V.(3)  (3.41) 

P  yj S2+p2/C2d 

which  contains  only  the  two  unknown  functions  E+(s)  and  V_(s),  and  now  the  Wiener-Hopf 
technique  can  be  applied  as  follows.  Suppose  that  the  function  L(s )  is  defined  and  factored 


as 


L(s) 


L-(*)  =  ..  t  F(S>P) 

L+($)  ^Js 2  +  p2/cd 


(3.42) 
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then  equation  (3.41)  becomes 

i 

y e-islL+(s )  -  E  +(s)L+(s)  =  L_(s)U_(s).  (3.43) 

Assume  that  the  function  D(s)  can  be  decomposed  as 

^e~islL+(s)  =  D(s)  =  D+(s)  +  D-{s),  (3.44) 

P 

then  equation  (3.43)  becomes 


D+(s)  -  E  +{s)L+(s)  =  L_(s)V_(s)  -  £L(s)  =  W(s).  (3.45) 


The  first  member  of  this  equation  is  analytic  in  the  upper  half  of  the  s— plane  Im(s)  =  A  >  A_ 
and  the  second  member  is  analytic  in  the  lower  half  of  the  s— plane  Im(s)  =  A  <  A+. 
Therefore,  if  A+  >  A_  the  regions  of  analyticity  overlap.  Using  the  Liouville’s  theorem  to 
determine  W (s),  solutions  for  E+(s)  and  V_(s)  can  be  found. 

After  algebraic  manipulation  it  is  found  that  the  function  F(s,p)  reduces  to 


F(s,p)  = 


{c22Cl27l72  -  7172(1  +  Cx2)c?2s2  + 


£vW^^V5(7i  +  72)  (i  +  C12) 

(cns2  +  p2/c2s)  [022(71  +  7 2)  +  7172(1  +  Ci2)c22  +  S2(ci2  +  c\ 2  -  cuc22)  -  C22p2/c2] } 


The  only  zeros  of  F(s,p)  are  of  the  form  s  =  ±ip/cR  where  cr  is  the  Rayleigh  wave  speed. 
This  can  be  seen  by  substituting  s  =  ip/v  in  F(s,p),  letting  F(s,p)  =  0  and  dividing  by  the 
non-zero  factors,  then  F(ip/v,p)  =  0  reduces  to 


jc 22  | 

(  C11C22  -  c\  2 

1  v 2  v2 

1  v2 

/  Cll 

\  c22 

cV 

N 

h 

oTbO 

1 

ofto  I 

CllC2 

which  is  the  Rayleigh  function  for  orthotropic  materials  [10].  The  roots  of  this  function  are 
v  =  ±cR. 

Consequently,  the  first  step  in  factoring  L(s)  is  to  define 


F(s) 


Fj^p) 

s2+P2/c2r 


(3.46) 
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Im(s) 


Figure  3.2:  Branch  cuts  of  F(s)  in  the  s— plane. 


It  can  be  shown  that  F(s )  -4  1  as  s  -4  oo,  (the  constant  £  in  (3.29)  was  chosen  to  make  this 
possible).  The  function  F(s)  is  regular  and  F(s)  ±  0  in  the  s- plane  cut  as  shown  in  figure 
3.2,  the  only  singularities  are  the  branch  points  shared  with  7!  and  72.  Where  the  branch 
points  of  71  and  72  are 


for  71;  s 

for  72;  s 


ip  _  ±ip 
yj Cll  Cs  Cd 


It  is  well-known  that  factorization  is  accomplished  most  directly  for  functions  that  ap¬ 
proach  unity  as  |s|  — >  00  and  that  have  neither  zeros  nor  poles  in  the  finite  plane;  F(s)  is 
an  example  of  such  a  function.  Therefore,  using  Cauchy’s  integral  formula  it  can  be  shown 


that  [3] 

Ct(s)  =  exP  Jv±  ~~~  ^z} 

where  F(s)  =  F+(s)F-(s)  and  T_  (r+)  is  the  contour  enclosing  the  branch  cut  between 
+ip/cd  and  +ip/cs,  (— ip/cd  and  — ip/cs ).  Using  the  fact  that  F(s)  —  F(s),  one  can  write 


F±(s)=  exp 


—  1  fl/CS 

—  /  tan 
7T  J  1/cd 


/Im[F(ipu;)]\  dw 
\Re[F(zptc)]/  w  =F  ” 


Note  that  by  making  s  —  ip(  in  this  equation,  F±(s  =  ipQ  =  F±( £)  becomes  a  function  only 
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of  £•  That  is 


(3.47) 


F±(Q  =  exp 


f  — 1  rx!c»  dw 

<  —  J  tan  1  [jR(io)] 


where 


R(w)  = 


1  /Cd 

Im[F(ipw)] 


w±C 


R e[F(ipw)] 

Returning  to  the  factorization  of  L(s)  we  have 

rM  _  L-(s)  _  Ms)(s2+P2/cR ) 

L(s)  -  z^r"i2f 

F+(s)F„{s)(s  +  ip/cR)(s  -  %p/cR) 


A*12£- 


y's  +  ip/cayfs  -  ipfcd 


therefore 


(3.48) 

(3.49) 


L-(s) 

(3.50) 

y's  -  ip/cd 

L+(s) 

sjs  +  ip/ cd  1 

(s  +  ip/cR)  F+(s) 

(3.51) 

Using  the  sum  splitting  formula  [7]  for  the  function  D(s)  defined  in  (3.44)  along  with  the 
result  (3.51)  it  is  found  that 


D+(s) 

D-(s) 


1  foo+iXo  D(z)  <Jq  1  foo+iXo  (z  +  ip/Cd)l!2e  tzl 

~CIZ  — 


2m  J—oo+tXo  z 
=  D(s)-D+(s) 


/_ 


p  2m  J-oo+iXo  (z  —  s)(z  +  ip/cR)F+(z) 


dz,  (3.52) 
(3.53) 


where  Ao  is  such  that  A_  <  A0  <  A+. 

Each  side  of  equation  (3.45)  is  analytic  in  one  of  the  overlapping  half  planes,  and  the  sides 
coincide  on  the  strip  of  overlap.  Consequently  each  side  of  (3.45)  is  the  analytic  continuation 
of  the  other  into  its  complementary  half  plane;  so  that  the  two  sides  together  represent  one 
and  the  same  entire  function  TU(s).  The  entire  function  will  be  determined  by  its  behavior 
at  |s|  —>  oo  which  is  related  with  the  behavior  of  physical  quantities  near  x  =  0.  First  note 
that  T+(s)  ~  s-1/2  and  T_(s)  ~  s1/2  as  |s|  — >  oo,  and,  D+(s )  and  (s)  in  (3.52)  and  (3.53) 
are  bounded  in  their  respective  planes  of  analyticity  and  vanish  at  infinity.  Furthermore, 
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a*+(x,p)  is  expected  to  be  square  root  singular  as  x  — )■  0+  and  v*_(x,p)  is  expected  to  vanish 
as  x  -»  0“  to  ensure  continuity  of  displacement.  Consequently  from  the  Abel  theorem  [7] 
relating  asymptotic  properties  of  transforms 


lim  xlt2o*,  (x,p)  ~  lim  s1/,2S+(s) 

x— >0+  +  s~>°° 

lim  \x\-qv*_(x,p)  ~  lim  |s|1+9V_(s) 

rc— >0~  s->-oo 

for  some  q  >  0.  Therefore,  it  is  expected  that  £+(s)  ~  s-1/2  and  VL (s)  ~  s~1-9  as  |s|  ->  oo, 
thus  the  products  £ +(s)L+(s)  and  L_(s)V_(s)  vanish  at  infinity.  Therefore,  each  side  of 
(3.45)  vanishes  as  |s|  — >  oo  in  the  corresponding  half  planes.  According  the  the  Liouville’s 
theorem,  a  bounded  entire  function  is  constant.  In  this  case,  W (s)  is  bounded  in  the  finite 
plane  and  kF(s)  0  as  |s|  -*  oo  so  that  the  constant  must  be  zero;  thus,  W(s)  =  0.  By 
using  (3.45)  and  (3.52),  the  functions  of  interest  are  then  given  by 

S+W  =  (3.54) 

=  £g.  (3.55) 


3.3.2  Stress  Intensity  factor 

To  find  the  stress  intensity  factor,  an  asymptotic  expression  for  the  normal  stress  near  the 
crack  tip  is  sought.  The  Abel’s  theorem  relating  asymptotic  expressions  between  a  function 
and  its  Fourier  transform  is  the  following  [7] 

lim  y/x  al(x,p)  =  lim  £+(s).  (3.56) 

x— >0+  V  +v  s-»+oo  V  7T 


Clearly,  the  behavior  of  £+(s)  as  s  — >  oo  is  needed.  Note  from  (3.51),  (3.52)  and  (3.54) 


that 


<70  1  1  f°°+iX 0  ( z  +  ip/cd)l/2e ‘ 


/_ 


izl 


5]  (s)  —  u  _ 

+  p  S1/2  2ni  J-oo+i\o  (z +  ip/cR)F+(z) 

by  making  z  =  ipC  this  equation  becomes 

-a0Vi  1  1  /•To+io°  (C  +  l/c<i)1/2epCi 


dz  as  s  — »  oo 


(3.57) 


+^~  pV 2  sV2  2-iri  L-ioo  (C  +  l /cR)F+(C)d^  aS  S  °° 


(3.58) 
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where  r0  is  real  and  — 1/q  <  t$  <  0.  Using  this  relation,  the  definition  of  the  stress  intensity 
factor  gives 


*?(P) 


lim  \Z2kx  al(x,p)  =  lim  e 

x^0+  +  s-»+oo 


-in'4V2s  E+(s) 


V2a0 

pl/2 


(  I  rro+ioo  } 

{mL -  /(0eK  dC} 


? 


(3.59) 

(3.60) 


where 


/(O  = 


(C  + 

(C  +  1/cr)F+(C)' 


(3.61) 


To  perform  the  inversion  of  Kj(p)  the  Cagniard-de  Hoop  method  can  be  applied.  The 
central  idea  of  the  Cagniard-de  Hoop  scheme  [1,  3]  is  to  convert  the  integral  in  (3.60)  to  a 
form  which  allow  inversion  of  the  one-sided  Laplace  transform  by  observation.  The  path  of 
integration  is  modified  to  form  a  closed  contour,  figure  3.3,  such  that  the  integrand  in  (3.60) 
is  analytic  inside  of  this  contour.  Applying  Cauchy’s  theorem,  Jordan’s  lemma  and  the  fact 
that  1(C)  =  /(C)  one  can  write  K}  ( p )  as 

*?(p)  =  {;  fZ, Im  [/(c)1  eK'  *}  •  (362) 

Letting  C^  =  —  rj  this  equation  becomes 


K*i(p)  =  L  Im  H ^  -  l/Cd^e  PT>  dr}-  (3-63) 

and  the  inversion  of  the  Laplace  transform  becomes  obvious.  Equation  (3.63)  is  a  product 
of  two  transforms,  so  that  Ki(t)  is  a  convolution  of  the  inverse  of  the  two  transforms,  i.e. 


Ki(t)  =  f/l  Im  {- 

7T  V  Tri  j l/crf  (  I 


(1  /Cd  ~  v)1/2 


drj 


(3.64) 


(1/CR  -  I))f-(l))(t/l  -  T])'/2 

where  it  has  been  used  the  result  F+(~ C)  —  /’-(C)- 

Note  that  for  time  t  >  l/cs,  equation (3.64)  can  be  evaluated  by  complex  integration  as 
shown  by  Freund  [3].  It  equals  the  sum  of  contributions  from  the  pole  at  C  =  1  /cr  and  the 
integral  taken  along  a  closed  contour  of  infinitely  large  radius  by  using  Cauchy’s  theorem. 
The  final  result  for  the  intervals  l/cs  <  t  <  1/cr  and  t  >  1/cr  is 


Ki(t)  =  a0\ 


1 


(1  /cr  —  1  /Cd) 


(i/cR-t/iy/*FZi/cR) 


H(1/cr  -  t/l ) 


(3.65) 
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Figure  3.3:  Contour  of  integration  to  evaluate  the  integral  in  (3.60)  in  the  (—  plane. 

The  second  term  in  equation  (3.65)  becomes  zero  for  t  >  l/cR  and  the  remaining  first  term  is 
the  corresponding  static  solution  for  a  semi-infinite  crack  under  normal  concentrated  forces 
applied  on  its  faces.  For  the  interval  l/cj,  <  t  <  l/cs,  the  stress  intensity  factor  is  obtained 
by  numerical  evaluation  of  the  integral  in  equation  (3.64). 

3.4  Results  and  Conclusions 

The  stress  intensity  factor  Ki(t )  computed  according  to  equations  (3.64)  and  (3.65)  are 
shown  in  figures  3.4  to  3.6  for  different  materials.  The  mechanical  properties  used  for  the 
analysis  are  given  in  table  3.1. 

Figure  3.4  shows  Ki(t)  for  an  isotropic  material,  in  this  case  epoxy.  The  normalization 
factor  is  the  long  time  limit  Kj( oo)  =  aoy2/i:l  and  the  normalized  time  is  cst/l.  The 
behavior  illustrated  in  the  plot  is  the  typical  for  isotropic  materials  as  described  by  [3]. 
After  the  arrival  of  the  dilatational  wave  generated  by  the  load,  the  stress  intensity  factor 
takes  on  a  small  negative  value,  reflecting  the  tendency  of  the  crack  faces  to  move  toward 
each  other.  The  small  negative  value  persists  until  the  shear  wave  front  arrives  at  time 
t  =  l/cs.  Thereafter,  the  stress  intensity  factor  decreases  rapidly  to  a  negative  square  root 
singularity  at  time  t  =  l/cR,  which  is  the  instant  of  the  arrival  of  the  Rayleigh  wave  traveling 
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along  the  crack  faces  from  the  load  to  the  crack  tip.  For  time  t  >  l/cR,  the  stress  intensity 
factor  takes  on  the  constant  value  which  is  the  equilibrium  stress  intensity  factor 

for  the  specified  applied  loading. 

Figures  3.5(a)  and  3.5(b)  show  the  stress  intensity  factor  history  for  a  graphite/epoxy 
composite  with  fibers  parallel  to  the  x— axis  and  y— axis,  respectively.  A  similar  behavior 
is  observed  to  that  shown  by  the  isotropic  material.  After  the  time  t  =  l/cd ,  Kj(t)  takes 
a  small  negative  value  until  t  —  l/cs  and  then  decreases  rapidly  to  a  singularity  at  time 
t  =  1/cr.  After  the  arrival  of  the  Rayleigh  wave  at  the  crack  tip  emitted  by  the  load,  Kj(t) 
takes  its  long  time  limit  AT/( oo)  =  a0^2/7rl.  The  difference  between  Kr(t)  for  isotropic  and 
orthotropic  materials  are  in  the  small  negative  value;  for  the  orthotropic  case  this  value  is 
almost  zero  over  an  ample  range  of  the  interval  l/cd  <  t  <  l/cs.  This  fact  explains  why  figures 
3.5(a)  and  3.5(b)  are  almost  identical.  That  is,  even  though  cd  is  different  for  the  composite 
with  fibers  parallel  to  the  a:— axis  and  to  the  y— axis,  the  arrival  of  the  dilatational  wave 
front  on  the  crack  tip  is  almost  imperceptible  and  K^t)  exhibits  no  appreciable  difference 
in  both  cases. 

Figure  3.6  shows  the  stress  intensity  factor  history  for  E-glass  epoxy  composite.  The 
behavior  of  Ki(t )  is  very  similar  to  that  shown  by  graphite/epoxy.  The  difference  is  the 
instant  when  the  singularity  occurs.  Again  a  very  small  negative  value  is  observed  for  Kd{t) 
in  the  interval  l/cd  <  t  <  l/cs.  Under  closed  inspection,  however,  as  shown  in  figure  3.6(b), 
it  is  observed  that  the  shape  of  the  curve  after  the  arrival  of  the  dilatational  wave  front  is 
similar  to  that  for  isotropic  material  but  with  different  scale. 

Finite  cracks  can  be  considered  as  semi-infinite  for  a  short  period  of  time.  Figure  3.7 
shows  the  stress  intensity  factor  history  for  a  finite  crack  under  impact  concentrated  loads 
on  its  faces  for  orthotropic  material  obtained  by  Rubio-Gonzalez  and  Mason  [8]  using  a 
different  method  of  solution  that  involves  numerical  solution  of  a  Fredholm  integral  equa¬ 
tion  for  the  Laplace  transform  of  the  stress  intensity  factor  followed  by  numerical  Laplace 
transform  inversion.  This  numerical  solution  suffers  from  spurious  oscillations,  due  to  the 
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Laplace  transform  inversion  technique,  at  the  discontinuity  in  Kj(t)  and  at  t  —  0.  The 
closed  form  solution  is  much  simpler  to  evaluate  because  of  its  simple  algebraic  and  integral 
representation,  equations  (3.64)  and  (3.65),  and  does  not  exhibit  spurious  oscillations.  The 
two  methods  show  a  good  agreement,  (Figure  3.7(b)),  in  the  appropriate  time  interval. 
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Isotropic  material 


Figure  3.4:  Stress  intensity  factor  history  for  a  semi-infinite  crack  in  isotropic  material, 
concentrated  loads. 
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(~) '»/(»)  > 


Figure  3.5:  Stress  intensity  factor  history  for  a  semi-infinite  crack  in  graphite/epoxy  com¬ 
posite,  concentrated  loads.  Fibers  parallel  to  the  (a)  a:— axis,  (b)  y— axis 


Graphite 

Epoxy 

E-Glass 

Epoxy 

Boron 

Epoxy 

Epoxy 

Cll 

20.77 

8.38 

KS9 

C22 

2.18 

2.29 

C12 

| 

)Ui2  (GPa) 

6.4 

mum 

p  (Kg/m3) 

1600 

2100 

1990 

1260 

Table  3.1:  Mechanical  properties  used  for  the  analysis. 
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(a)  (b) 


Figure  3.6:  Stress  intensity  factor  history  for  a  semi-infinite  crack  in  E-glass/epoxy  compos¬ 
ite,  concentrated  loads,  (b)  Under  closer  inspection. 


Finite  Crack  with  Concentrated  Loads  at  x^.6,  Orthotropic  Material  Comparison  of  results  for  finite  and  semi-infinite  cracks 


Figure  3.7:  Comparison  between  the  stress  intensity  factor  for  a  finite  crack  of  length  2 a  [8] 
and  a  semi-infinite  crack  under  concentrated  loads  in  orthotropic  materials. 
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Chapter  4 

Dynamic  Stress  Intensity  Factor  for 
Semi-infinite  Cracks  in  Orthotropic 
Materials  due  to  Concentrated  Shear 
Impact  Load 


co-authored  with  C.  Y.  Wang  and  C.  Rubio- Gonzalez  and  accepted  for  publication  in  Inter¬ 
national  Journal  of  Solids  and  Structures 

The  transient  elastodynamic  response  of  an  orthotropic  material  due  to  concentrated 
shear  impact  loads  on  the  faces  of  a  semi-infinite  crack  is  examined,  and  the  solution  for 
the  stress  intensity  factor  history  around  the  crack  tip  is  found.  Laplace  and  Fourier  trans¬ 
forms  together  along  with  the  Wiener-Hopf  technique  are  employed  to  solve  the  equations  of 
motion.  Even  though  the  problem  has  characteristic  length,  it  has  been  shown  in  previous 
works  that  the  Wiener-Hopf  technique  can  be  applied.  The  asymptotic  expression  of  stress 
near  the  crack  tip  is  analyzed,  which  leads  to  the  dynamic  mode  II  stress  intensity  factor. 
The  results  are  presented  for  orthotropic  materials  as  well  as  for  an  isotropic  material. 
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1  Introduction 


The  growing  use  of  composites  in  many  engineering  applications  demands  a 
fundamental  understanding  of  the  response  of  cracked  orthotropic  bodies  to  impact 
loading.  By  using  the  integral  transform  methods,  the  dynamic  behavior  for  a  finite  crack 
under  concentrated  loads  [1,2]  as  well  as  for  a  semi-infinite  crack  under  concentrated  or 
uniform  normal  loads  [3,  4]  on  its  faces  in  orthotropic  materials  have  been  analyzed  by 
Rubio-Gonzalez  and  Mason.  In  the  former,  the  finite  crack  problem  is  reduced  to  a 
Fredholm  integral  equation  in  the  Laplace  transform  domain  which  is  solved  numerically, 
and  in  the  latter  for  the  Weiner-Hopf  technique  was  used  to  find  the  stress  ahead  of  the 
crack  tip  and  the  displacement  on  the  crack  faces.  The  asymptotic  expression  for  the 
stress  near  the  crack  tip  leads  to  the  stress  intensity  factor  K(t) . 

Many  solutions  for  various  loadings  of  cracks  in  orthotropic  materials  have  been 
found  by  applying  transforms  to  the  displacement  formulation  of  the  equations  of  motion. 
Where  previous  researchers  (for  example,  [5])  following  this  approach  had  solved  the 
resulting  dual  integral  equations  using  the  method  of  Sneddon  [6,  7],  Rubio-Gonzalez 
and  Mason  [3]  solve  the  same  equations  by  converting  the  dual  integral  equations  to  a 
Weiner-Hopf  equation.  Solution  of  that  equation  requires  only  a  straight  forward 
application  of  the  method  of  Noble  [8].  In  this  paper,  the  problem  of  concentrated  shear 
impact  loads  on  the  faces  of  a  semi-infinite  crack  in  orthotropic  materials  is  analyzed 
using  the  same  method.  While  it  was  been  thought  that  the  Weiner-Hopf  technique  could 
not  be  applied  to  problems  such  as  this  [9,  10],  this  approach  has  recently  been  shown  to 
be  valid  in  such  cases  [11],  Through  Laplace  and  Fourier  transforms  combined  with  the 
Weiner-Hopf  technique,  a  solution  for  the  stress  ahead  of  the  crack  tip  is  sought.  The 
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asymptotic  expression  for  the  stress  near  the  crack  tip  leads  to  the  mode  II  stress  intensity 
factor,  K„(t). 


2  Description  of  the  problem 

The  plane  problem  of  an  infinite  orthotropic  body  containing  a  semi-infinite  crack 
is  considered  in  Figure  1.  Let  E.,  pi}  and  v..  (i,  j  =  1,  2,  3  )  be  the  engineering  elastic 

constants  of  the  material  where  in  the  indices  (1,  2,  and  3)  respectively  stand  for  the 
directions  (x,  y,  z)  of  a  Cartesian  coordinate  system  chosen  to  coincide  with  the 
orthotropic  axes  of  material.  The  crack  faces  are  suddenly  loaded  by  a  pair  of 
concentrated  shear  forces  with  a  magnitude  (q)  located  a  distance  (a)  away  from  the  crack 
tip,  as  shown  in  Figure  1. 


2.1  Equations  of  motion 

The  following  problem  is  limited  to  a  two-dimensional  case  with  the  inertial  effect  in 
the  x-y  plane  only.  By  setting  all  the  derivatives  with  respect  to  z-direction  to  be  zero,  it  is 
readily  shown  that  the  displacement  equations  of  motion  [12]  reduce  to 


d2u  <?2u  <?2v  1  <?2u 

c"1?+J7h 

d2v  d2' 


-v  „  ,  d2u 

^•  +  CK^r+(l  +  C,2)— =  - 


2  dt2  ’ 


1  d2v 


(1) 


2  dt 2  ’ 


where  u  and  v  are  the  x  and  y  components  of  the  displacement  vector.  In  the  orthotropic 
solids,  cs  =  represents  the  velocity  of  the  shear  wave  propagating  along  the 

principle  material  axes  in  xy-plane  and  p  is  the  mass  density.  Here  also,  we  introduce  the 
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velocity  cd  =  *fc^cs ,  which  will  appear  later,  representing  the  dilational  wave  speed 
along  the  x-axis.  The  non-dimensional  parameters  cu ,  cl2  and  c22  are  related  to  the 
elastic  constants  by  the  following  relations.  For  a  generalized  plane  stress  problem, 


c„  =• 


U„H-(E,/EI )v’„]  ’ 

C22  “(^2/^l)ClI  * 


^12  “  ^12^22  ^2iCll  • 


For  a  plane  strain  problem, 


cn  A  ^23^32)  * 

C22  =  ”"  7"(1  _  ^13^31)  » 

Mn  A 

i?,  Et 

ci2=—Lr(v2  ,+-2-v„vm), 

/^12  A 

A  =  l-vi2v21-v23v32-v31v'13-v,2v23v3,  -v13v2Iv32  , 

=  =  2,  3). 

Ay  A. 

Finally,  the  stresses  are  related  to  the  displacements  by  the  equations, 


O  o\x  ov 

j;-c"~d-x+c"Ty’ 

du  dv 

t  „  du  dv 
v  -  + 


—  =  c„ 


fin  dy  dx 


(2) 


(3) 


(4) 


2.2  Boundary  conditions  and  Initial  conditions 
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For  the  case  shown  Figure  1,  the  pair  of  spatially  concentrated  shear  loads  is  applied 
suddenly  on  the  crack  faces.  Exploiting  symmetry  and  limiting  ourselves  to  upper  half¬ 
plane,  y>0,  the  corresponding  boundary  conditions  become 

O y(x,  0,t)  =  0,  - 00  < x  <  +°° , 

Txy  (x,  0,  t)  =  -qH(t)S(x  +  a) ,  -  °°  <  x  <  0 ,  (5) 

u(x,  0,  t)  =  0  ,  0  <  x  <  +00  , 

where  q  is  the  force  per  unit  length  in  the  z-direction  of  the  opposed  line  loads  acting  on 

the  crack  faces  and  H(t)  is  the  unit  step  function.  The  dirac  delta  function,  <5(x  +  a) ,  has 

the  dimension  of  length'1 .  In  addition,  the  displacement  at  infinity  is  zero,  and  the  body 
is  stress  free  and  at  rest  everywhere  for  t  <,  0 . 

3  Integral  transforms 

3.1  Laplace  transform 

hi  equations  (1),  the  time  variable  may  be  replaced  by  application  of  Laplace 
transform 

f\s)=  f/(0  e-"dt,  f(t)  =  -  V  f  /* (*)  « : mds  =  -5-t  f /* is)  eads ,  (6) 

Jo  l  J Br  2d7C\  Ja-fr* 

where  Br  denotes  the  Bromwich  path  of  integration  which  is  a  line  parallel  to  the 
imaginary  axis  in  the  s-plane.  Thus,  using  the  Laplace  transform  representation,  the 
displacements  have  the  form 

u(x,  y,  t)  -  —V  f  s)  es'ds , 

2  71 1  JBr 

v(x,  y,  t)  =  — V  f  v*(x,  y,  s)  e“ds . 

2ni  JBr 
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Applying  the  relations  (7)  to  equations  (1)  and  assuming  zero  initial  condition,  the 
transformed  domain  equations  become 


<?2u*  d2u 


dx1  dy 


+-r-r-  +  (l  +  C12) 


<?V 

dxdy 


aw 

dx 2 


+  C22 


dW 

dy2 


+  (!  +  %) 


aw 

dxdy 


(8) 


Applying  the  relations  (7)  to  equations  (4),  the  transformed  domain  relations  between 
stresses  and  displacements  become 


du  dv' 

—  _  c„  —  h  cu  — 

/*,2  dx  dy 

<  du  <?v* 

jul2  C 12  +C* 

nn  dy  dx 


(9) 


Finally,  the  application  of  the  Laplace  transform  to  the  boundary  conditions  (5)  gives, 

O  *  (x,  0  ,  s)  =  0  ,  -  00  <  x  <  +00 , 

W (*»  0,  s)  =  —  q/s S(x  +  a) ,  -°°<x<0,  (10) 

u*(x,  0,  s)  =  0 ,  0<x<+°°. 


3.2  Fourier  transform 

To  obtain  a  solution  of  the  differential  equations  (8)  subject  to  conditions  (10),  the 
Fourier  transform  is  applied, 

f(a)  =  j  f(x)  e“"dx,  /U)  =  ^J/(t»)  e-‘°*d6)  (11) 

It  is  noted  that  the  boundary  conditions  (10)  are  explicitly  defined  only  on  half  of  the 
range  of  x.  Consequently,  the  Fourier  transform  can  not  be  applied  to  these  boundary 
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conditions.  To  remedy  this  situation,  the  boundary  conditions  must  be  extended  to  apply 
on  the  full  range  of  x.  Two-unknown  functions  u_'(x,s )  and  z+'(x,s)  are  introduced. 
The  function  u_’(x,s)  is  defined  to  be  the  x-  direction  displacement  of  the  crack  face, 
y=0,  for-°°<jc<0  and  0 <t<°°,  and  to  be  identically  zero  for  0 <*<+«>  and 
0 <t<oo.  Likewise,  the  function  z+\x,s)  is  defined  to  be  the  shear  stress  in  the  x 
direction  on  the  plane,  y=0,  for  0<jc<°°  and  0 < t < °o ,  and  to  be  identically  zero  for 
— oo < jc < 0  and  0<f<°°  .  With  these  definitions,  the  boundary  conditions  (10)  can  be 
rewritten  as 


o  *  (x,  0  ,  s)  =  0  ,  for  -oo  <  x  <  °o , 

tjcy\x,0,s)  =  -q/s6(x  +  a)  +  z+\x,s)  ,  for  -oo < x  < °° ,  (12) 

U  (X,0,  S)  =  u’  (X,  S)  ,  for  -oo  <  X  <  oo  . 

Thus,  the  Laplace  transformed  displacements  have  the  form 


u(x,  y,  s)  =  f  u*(fi),  y,  s)  e'laxdco  , 
2  K  J- 

v*(x,  y,  s)  =  ^-Jv*(ru,  y,  s )e'laxdco  . 


(13) 


Substituting  these  transforms  into  equations  (8),  the  functions  m*  and  v*  are  found  to 
satisfy  the  ordinary  differential  equations 


d2 u*  .  .rfv*  2  s*  A 

— (1  +  c22)ox— - (cuco  +—)u  =0, 

dy  dy  ct 

d2r  /t  ^  .rfu*  .  2  _ 

c^-— — (\  +  cn)ax— — (co  +— )v  =0. 
dy  dy  cs 


(14) 


The  solution  of  equations  (14)  under  the  condition  of  zero  displacement  at  infinity  is 
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(15) 


u '  {a,y,s)  =  Ai{a,s)e~ny  +  A2(a,s)e  ny , 

v *  (a,  y,s)  =  - A.  (a,  s)e ~ny  -  — -A2(a,s)e~ny  , 
a  a 

where  A,(<m,j)  and  A2(a,s)  are  arbitrary  functions,  and  a.(a,s)  stands  for  the 
functions 


,  N  cn<a2  +s2/cs2  -)2 
a,  (a,  s)  =  — - — 1 — — — 

(1  +  c12)r, 


(i  =  l,  2), 


(16) 


with  j  2  and  }  2  being  two  roots  of  the  quadratic  equation, 

+[(c122  +2cI2 -cucn)a2  -(1  +  c22)s2/c2]y2  +(ctia2  +  s2/c2)(a2  +s2/cl2)= 0  . 


(17) 

From  (17),  it  is  readily  found  that  the  functions  }  t(a)  and  }  2(a)  are  multiple 
valued  functions  of  a  in  the  complex  a  -plane  with  branch  points  at  a  =  ±is  /  cs  and 
a  =  ±is/Jc^cs  =  ±i s/Cj  ,  respectively.  It  can  be  shown  that  for  many  materials  that) , 
and  )  2  are  real  and  positive.  Thus,  the  expressions  for  the  displacements  in  the  Laplace 
transform  domain  become 


u  =—  f  (A1e"m  +A2e-r2y)e-imd(0 , 
2  k 


'•=-f 

2n 


(18) 


{a.A.eny  +a2A1e-ny) - da . 

a 


Substituting  these  displacements  into  the  constitutive  law,  (9),  we  obtain 


o 

Mn 


*  •  Mo o  •  MOO  -|GK 

—  =  ~^7cJ(A'e~ny  +A2e~ny)e~>c*d(D  +  l-^)  («,r.Vm  +a2r2A2e-r2y)—dco, 


-n>  +  Ae-ri>}e-i«*(i(0+i£n_  f  (a  y  Ale~ny  +  a  2y  2A2e~riy )- — da, 
jul2  2 n  J-  2k  J—  a 
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♦ 

f~  =  _  2 J (r,A,e  m  +  y2A2e'ny)e'iMd(°-^  J(celAle'r,y  +  a2A2e-ny)e-imdo) . 


t 
Mn 


(19) 

Substitution  of  the  expression  for  o'  the  second  of  equations  (19),  into  the  boundary 


condition  for  o\,  the  first  of  equations  (12),  yields 


A2{,a,s)  =  -^Al{co,s) , 

cJsaj1  -  c,2fl?2  _  cK  (c„fi?2  +  *7c,2  -  yj_ )  -  cn  (1  +  cn  )fl>: 
-cno)2  cn(cnO)2  +s2/c2  -y22)-cn(l  +  cn)6) 


(20) 


Using  this  result  in  the  expression  for  shear  stress,  the  third  of  equations  (19),  and  the 
expression  for  horizontal  displacement,  the  first  of  equations  (18),  we  obtain 


V  =-|jJ[(a,  -0(a2  +y2)e-™]Ale-“*da>, 

u  =—  f  [e'm  - 

2#  J™ 

Applying  the  remaining  boundary  conditions  yields 

-fcji*.  +yl)-JK<x2+y2)]Aie-“*da)  =  -q/sS(x+a)  +  T;(x,s ) , 

—  f  [l-j8]i4I«"'“,</fl>  =  «_*(*,*)  . 

In  J— 

By  Fourier  transform  inversion,  these  equations  become 

-//12[(«,  +j,)->5(a2  +?2)]A,  =-q/se'iwa  +E+(ta) , 

(1  -p)Ax=US(D), 


where 


£+(<y)=  |V(x,  *)«'“<&, 

Jo 

U_{a)  =  ^u-'(x,s)e“*dx. 


(21)  . 


(22) 


(23) 


(24) 
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Eliminating  A,  from  equation,  (23),  we  obtain  a  Wiener-Hopf  equation, 


Mn C  7  ==  j  v.m  =  -«/«""  +  £.(«>)  . 
Ja)2+s2/c,2 


(25) 


where  we  have  introduced  the  function 


Jo)2  +s2/c 2 

R(o),s ) = -  (1_^c*  +  r. )  -  + r2 )] , 


(26) 


c= 


{cncJ2iV12iV22-c122(l  +  c12W1iV2  + 


c22(l  +  c12)NliV2(A(I+Ar2) 

Cn [c22(A i  +  N2  )  +  cn(l  +  cl2)NlN2  +  (cI2  +cl2  — CjjCjj)]}, 


(27) 


Nia1  =T^— {(cuc22  ~cn  -2c,2)±[(cnc22  -cI22  -2c12)2  —  4c, , c22  ] 1/2 } .  (28) 

After  algebraic  manipulation,  it  is  found  that  function  R(co,s)  can  be  written  as 

-c»ii+cM‘  Hc"m‘ +iVc-!> 
[c22(r,2  +r22)  +  cn{\+cl2)yj2  +(c12  +cI22  -cuc22)©2  -c^/c,2]} . 

(29) 


From  the  physics  of  the  problem,  it  is  reasonable  to  assume  that  the  function 
i*(x,  5)  and  u_'(x,  s)  are  exponentially  bounded  at  infinity,  which  ensures  the  existence 
of  their  Fourier  transform  (24).  In  particular,  it  is  shown  by  Noble  [8]  that  if 
|t+*(x,5,)|  <  M{ex~x  as  x->+oo  then  X+(ty)is  analytic  in  Re(tu)  =  X>X_,  and  if 

|«_*  (x,  s)|  <M2ex+x  as  x  then  U _ ( co )  is  analytic  in  Re(ta)  =  A  <  A+ . 

The  function  R(o),s )  is  analytic  everywhere  in  the  complex  plane  except  for  at  the 
branch  cuts  of  ) .  between  Q)  =  ±is/cs  and  co  =  ±is/cd ;  it  is  single  valued  in  the  0)- 
plane  cut  as  shown  in  Figure  2.  It  can  be  shown  that  the  only  zeros  of  R((0,  s)  are  of  the 
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form  ±is/cR  where  cR  is  the  Rayleigh  wave  speed.  This  can  be  seen  by  substituting 

t 

co  =  is/v  in  R(0),s),  letting  R(o),s)  =0  and  dividing  by  the  non-zero  factors,  then 
R(i  s/v,  s )  =  0  reduces  to 


which  is  the  Rayleigh  wave  function  for  orthotropic  materials  [13].  The  roots  of  this 
function  are  v  =  ±cR . 


4  Wiener-Hopf  technique 

To  solve  the  Weiner-Hopf  equation,  (25),  we  use  the  Weiner-Hopf  technique  as 
follows.  The  function  L(co)  is  defined  and  factored  as 


L(co)  =  L+(co)L_(co)  -  MnC 


R(co,s ) 


^a)3+s2/ct2 


(31) 


then  the  Weiner-Hopf  equation  becomes 


L+(co)  L+(a>) 


(32) 


If  the  function  D(a)  can  be  written  and  decomposed  as 


D{co)  =  -  =  D+  (a)  +  D_ {co) , 

L+(co) 


(33) 


then  the  Weiner-Hopf  equation,  (32),  is  further  reduced  to 

D+(co)  +  £+  (<a)/L+  (oo)  =  L_ (co)U_ (co)-D. (oo)  =  W(ta) .  (34) 

The  first  expression  for  W(co)  in  (34)  is  analytic  in  the  right  half  plane,  Re(ftJ)  >  X_ , 
and  the  second  expression  is  analytic  in  the  left  half  plane,  Re(ffl)  <  .  If  A_>  A+,  the 
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regions  of  analyticity  overlap  and  by  invoking  analytic  continuation,  it  is  concluded  that 
W(co)  is  analytic  and  single-valued  in  the  whole  tf)-plane  shown  in  Figure  2. 
Furthermore,  invoking  the  extended  Louisville  theorem,  it  can  be  shown  [8]  that  if  W(o)) 
is  bounded  and  entire  and  W((ti)  ->  0  as  CD  -» °° ,  then  VK(ta)=0.  Hence,  we  can  solve  for 
the  transform  of  stress  ahead  of  the  crack  tip,  Z+(tfJ) ,  and  displacement,  U_{co) ,  behind. 


U_(a»  = 


DM 

lm  ' 


(35) 


Following  this  outline,  the  first  step  is  to  factor  L(a>) ,  (31),  by  defining 


F«o)  = 


R(o),s) 

cd2+s2/cr 2  ' 


Note  that  by  making  co  =  isg  in  equation  (36),  then  F(ffl)  is  changed  to  a  function  only 


(37) 


It  can  be  shown  that  F(f)  — >  1  as  Q  — >  °° ,  ( the  constant  C  in  (27)  was  chosen  to  make 
this  possible).  The  function  F(f)  is  regular  and  F(g)  ^  0  in  the  g  -plane  cut  as  shown  in 
Figure  3.  The  only  singularities  of  F(g)  are  the  branch  points  shared  with  } ;,  (i  =  1,  2) , 
i.e.  g  =  ±l/ca  and  g  =  ±\/cd.  It  is  well-known  that  factorization  of  a  function  is 
accomplished  most  directly  for  functions  which  approach  unity  as  |$|  -»  °°  and  which 


have  neither  zeros  and  nor  poles  in  the  finite  plane.  Indeed,  F(g)  is  an  example  of  such  a 
function.  Therefore,  using  Cauchy’s  integral  formula,  it  can  be  shown  that  [10] 
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FM 


log  F(& 


dt\. 


(38) 


where  T.  is  the  closed  counterclockwise  contour  enclosing  the  branch  points 
+l/Cj  and  +  l/c3 ,  and  also  T+  is  the  closed  counterclockwise  contour  enclosing  the 
branchpoints  —  1 Jcd  and  —  1/c, . 


Using  the  fact  that  F(g)  =  F(g),  =  -H(rj) ,  then 


Vc*  di]  1 

my  11 

v<j 


(39) 


where  H(rj)  =  tan 


Im[  Pm' 
Re[ F(u)b ' 


By  making  g  =  co/  is  in  equation  (39),  then  F±  (g)  is  changed  to  F±  (w) . 
Returning  to  the  factorization  of  L(co) ,  we  now  have 

+ 1  s/cl  «Ja>  -is/cs 


(40) 


Therefore, 


L((0)  =  M12C^j  IS!C.r)  F_{(0) , 
Jo-is/c, 

Jco  +  is/cs 


(41) 


Substituting  equation  (41)  into  equation  (33)  the  function  D(co)  can  be  obtained, 


-q/se-~  >+iVc,  g/se- 

LJ,o>)  (a+is/c,)  FJa) 

Through  Cauchy’s  integral  formula,  we  obtain 


(42) 
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D(ffl>=— f  j_r 

*  2  mJc±z-Ct)  s[2mJc±(z-(D)(z  +  is/cR)F+(z)  J 


,  (43) 


where  C_  is  the  closed  counterclockwise  contour  enclosing  the  branch  points 
+isjcd  and  +  is/ct ,  and  C+  is  the  closed  counterclockwise  contour  enclosing  the  branch 
points  -is/Cj  and  -  is/  cs . 

Thus,  we  can  show  that  W(6))=0,  and  using  equation  (35),  we  can  readily  obtain 
£+  (ffl)  and  U_ (co) .  We  only  give  £+ (co)  here, 


(44) 


5  Stress  Intensity  Factor 


To  find  the  stress  intensity  factor,  an  asymptotic  expression  for  the  shear  stress  near 
the  crack  tip  is  sought.  A  well-known  result  relating  asymptotic  expressions  between  a 
function  and  its  Fourier  transform  (Abelian  theorems)  [8]  is 


lim  -yfxT  ' * (x,  s)  =  lim  e 

JC-+0+ 


-ixi  4 


(45) 


Thus,  the  usual  definition  of  the  stress  intensity  factor  gives 

K’(s)  =  lim  *J2kxt"(x,s)  =  lim  c'"'4 -s/2 6)S+(6))  =  lim  0(6)) , 

0(6))  =  e'"r/4V26)2+(6)) . 


As  6)  +oo ,  0(6))  becomes 


0(6))  =  - 


42g  e~i/4 
^  2m 


*  Jz+Ts/C, 

c+  (z  +  is/cR)  F+(z) 


and,  by  making  z  =  ,  equation  (47)  becomes 


(46) 


(47) 
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<D(<y)  =  - 


_  J2q  If  Vg+Vc,  c* 


LLf 

2m  Jr+ 


(f+v«.)f;(s) 


as  6J  -4  +oo 


Thus, 


K:,{s)  =  -^2q\^y{g)dg\, 


where 


8(0  = 


V^+l/c,  e'Q 
(g+1 /cR)t(0s'n 


Note  that  g(g)  =  g(£) ,  and  g(g)  -»  0  as  |$|  — >  <»  in  the  left  half  plane.  Also,  note  that 
F+  (-g)  =  F_(g) .  Consequently,  equation  (49)  can  be  written  as  the  real  integral 


k;m=  -jm-  \L[g(7f)v^ 

yjt  J-iicj 

it  Jucj  l  (1 /c„-tj)F_(V)  s> 


where  if  is  a  positive  real  number. 

Applying  the  inverse  Laplace  transform  to  fQ(s),  using  the  convolution  formula  and 
letting  rf  a  =  t ,  the  stress  intensity  factor  Kn(t)  can  be  written  as  follows 


K„ (t)  =  -qJZTJSi --  f  Im 

7t  Jucd 


(l/cR-T])Jt/a-riF_(7j) 


drj},  (51) 


The  integrand  of  (51)  has  a  first  order  singularity  at  g  =  l fcR  and  branch  points  at 
g  =  1 Jcd,  g  =  1 jcs  and  g  =  g0  =  t/a  ,  which  is  shown  in  Figure  4.  Here,  pay  attention  to 
those  points  at  g  =  1 jcd,  g  =  1/c,,  g=g0-t/a  and  q  =  1 JcR  ,  which  are  critical  of  the 


integrand  of  (51)  as  shown  in  the  Figure  4. 
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For  t  /  a  <  1  /  cd ,  no  waves  generated  by  the  applied  loads  have  arrived  at  the  crack 
tip,  hence  the  stress  intensity  factor  Kn(t)  is  identically  zero. 

For  l/crf<f/a<l/c,,a  direct  evaluation  procedure  for  the  stress  intensity  factor 
Kn (0  can  not  be  applied  to  the  integral.  It  has  to  be  obtained  by  numerical  evaluation. 

For  tla>\lci,  the  integrand  of  (51)  is  analytic  in  the  entire  £-plane  cut  along 
l/cd  <  Re(£)  <t  /  a,  Im(Q=0,  except  for  the  pole  point  at  g  =  1 jcR  .  In  this  case,  the 

path  of  integration  can  be  closed  around  the  right  side  of  the  branch  cut  shown  in  Figure 
4.  The  integral,  (51),  takes  the  form 

K„it) 

where  r(£0)  denotes  a  closed  counterclockwise  contour  embracing  the  branch  cut  of  the 

integrand.  Cauchy’s  integral  formula  can  now  be  applied  to  show  that  value  of  the  integral 
in  (51)  is  equal  to  the  value  of  the  integral  taken  along  a  closed  counterclockwise  circular 
path  of  indefinitely  large  radius,  (Q),  in  £-plane  shown  in  Figure  4.  Therefore, 


W2  f  yi /c,-g  i  i 

2 mjm  Jr(fo)  (1 /cR  -  g)  Jt/a-g  F_(g) 


(l/c^-l/c,)1'2 

<a/cR-t/a)mF(l/cR) 

0 


if  1/ cs  <t / a<l  / cR 
if  t  !  a>  1  /  cR 


(53) 


The  first  term  on  the  right  side  of  (53)  is  the  contribution  from  the  closed  contour  of  the 
large  radius  and  the  second  term  is  the  residue  contribution.  In  addition,  The  first  term  in 
(53)  is  the  corresponding  static  solution  for  a  semi-infinite  crack  under  shear  concentrated 
forces  applied  on  it  faces,  Ku  (t)  =  -q^Jlfm  . 
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6  Results  and  Discussion 


Using  equations  (52)  and  (53),  the  stress  intensity  factor  K„  (t)  can  be  calculated. 
The  computed  results  are  shown  in  Figure  5-8  for  several  materials,  and  the  mechanical 
properties  [14]  used  for  the  analysis  are  given  in  Table  1.  The  following  tips  should  be 
taken  into  account  when  performing  the  numerical  evaluation  of  the  integral  in  (51). 

(1)  Theoretically,  R(a,s)  in  equation  (26)  can  include  either  the  factor 

tJcq2  +  s2/c2  or  tJco2  +  s2/c2  ,  and  we  should  have  F(g )  =  1  as  .  However, 

experience  shows  that  if  the  factor  Jo2  +s2/c/  is  selected  in  R(a,s) ,  it  is  difficult  to 
numerically  make  F(g)  =  F+(g)F_(g)  at  points  in  the  regions  close  to  [1 Jcd  ,1 jcs  ]  and 
[-1 Jcs ,  -l/cj  ],  where  F(g),  F_(g)  and  F+(g)  are  respectively  calculated  using  equations 

(37)  and  (38).  This  is  true  even  though  F(g)  =  1  and  F(g)  =  F+(g)F(g)  at  points  very 
far  from  [l/cd  ,1/c,  ]  and  [-1/c, ,  -l/cd  ].  For  the  problem  with  concentrated  normal 

loads  [3],  the  factor  ^G)1  +s2/c/  should  be  selected  in  R(a,s) ,  while  the  factor 
tJg)2  +s2/c2  should  be  selected  in  R(co,s)  for  the  problem  in  this  paper.  Otherwise, 
there  may  be  large  errors  in  the  calculations. 

(2)  It  is  best  to  adopt  the  formula  F(f)  =  F(g)/P+(g)  to  indirectly  calculate  the 

values  of  F_(g)  along  [1 Jc4 ,1/c,  ]  because  F_(g)  is  not  analytic  along  [1/c^  ,1/c,  ].  At 
g  =  l/cR  the  value  of  F_  (1 /cR )  can  be  directly  obtained  by  using  the  equation  (39). 
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The  numerical  results  presented  correspond  to  isotropic  and  orthtropic  materials. 
Properties  of  transversely  isotropic  materials  with  fibers  parallel  and  perpendicular  to  the 
x-axis  are  considered  in  the  orthtropic  case. 

For  isotropic  materials,  as  mentioned  in  [3],  note  that  the  mechanical  parameters 
make  } ,  =  y  2  and  further  lead  to/5  =  1  and  R(cq,s)  not  defined.  Consequently,  in  order 
to  obtain  the  results  for  the  isotropic  case  from  the  orthotropic  formulation,  we  let 
El  =  E,  E2  =  (1  -  e)E ,  v12  =  Vjj  =  v  and  fln  =  E2/ 2(1  +  v)  where  E  and  v  correspond 

to  the  isotropic  properties  and  f  is  a  small  quantity  with  £  « 1 . 

Figure  5  shows  K„  (t)  for  an  isotropic  material  in  the  plane-strain  case.  The 

normalization  factor  is  the  long  time  limit  /£H(°o)  =  -q^Tfm  and  the  normalized  time  is 
cd  t/a.  The  behavior  illustrated  in  Figure  5  is  as  expected  from  the  surface  displacements 

in  the  calculated  solution  of  Lamb’s  problem  where  a  concentrated  shear  load  is  applied 
on  a  half-space  [15].  Horizontal  displacements  on  the  surface,  generated  by  the  point 
load,  arrive  at  the  crack  tip  with  the  arrival  of  the  dilatational  wave  causing  a  rapid 
increase  in  the  stress  intensity  factor.  In  addition,  the  results  in  Figure  5  are  very  similar 
to  that  of  the  problem  in  [2],  before  the  arrival  of  the  second  dilatational  wave,  shear 
wave  and  Rayleigh  wave,  where  a  finite  crack  in  infinite  orthotropic  body  is  impacted  by 
two  pairs  of  concentrated  shear  loads.  After  the  arrival  of  the  dilatational  wave,  a  rapid 
increase  in  Kn(t)  is  observed,  and  then  the  stress  intensity  factor  gradually  decreases 
until  the  shear  wave  front  arrives  at  t  =  a/c1 .  Thereafter,  the  stress  intensity  factor 
decreases  rapidly  to  a  negative  square  root  singularity  at  t  =  a/cR  ,  which  is  the  instant  of 
the  arrival  of  the  Rayleigh  wave  traveling  along  the  crack  faces  from  the  point  of 
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application  of  load.  For  tim et>a/cR  ,  the  stress  intensity  factor  takes  on  the  constant 
value  -qjlfnia ,  which  is  the  equilibrium  stress  intensity  factor  for  the  specified  loading. 
Obviously,  the  dilatational  wave  has  an  important  effect  on  Kn  (t)  under  applied 
concentrated  impact  shear  loads,  unlike  applied  concentrated  impact  normal  loads  where 
the  dilatational  wave  has  a  minimum  effect  on  Kf(t)  [1,3].  The  stress  intensity  factor  for 
an  isotropic  material  in  the  plane-stress  case  is  identical  to  the  plane-strain  case  except  for 
differences  in  the  dilatational  wave  speed. 

Figure  6(a)  and  6(b)  show  the  stress  intensity  factor  history  for  a  graphite-epoxy 
composite  in  the  plane-strain  state  with  fibers  parallel  to  the  x-axis  and  y-axis, 
respectively.  Both  figures  are  similar  to  the  results  for  the  isotropic  material  shown  in 
Figure  5,  except  that  there  is  a  more  uniform  plateau  from  t  =  a/cd  to  t  =  a/cs  and  a 

steeper  drop  at  t  =  a/cR.  In  Figure  6(b)  the  plateau  is  lower  than  in  Figure  6(a)  and  the 
duration  of  the  plateau  is  longer  due  to  the  dramatic  difference  in  cd  for  these  two  cases. 

The  Rayleigh  wave  speeds  are  approximately  the  same;  the  shear  wave  speeds  are  exactly 
the  same. 

As  expected,  the  behavior  illustrated  in  Figure  6(b)  is  the  same  as  that  for  the 
corresponding  case  of  Rubio-Gonzalez  and  Mason  [2]  before  the  second  set  of  waves 
arrives  crack  tip.  The  behavior  illustrated  in  Figure  6(a)  can’t  be  directly  compared  with 
that  of  the  corresponding  case  of  Rubio-Gonzalez  and  Mason  [2]  because  in  that  work  the 
second  dilatational  wave  arrived  at  the  crack  tip  before  the  arrival  of  the  first  shear  wave 
and  Rayleigh  wave  leading  to  very  different  loading  conditions.  To  the  knowledge  of  the 
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authors,  this  is  the  first  presentation  of  results  for  this  case  that  are  not  disturbed  by 
reflected  or  secondary  waves. 

Figure  7  and  8(a)  show  the  plane  strain  stress  intensity  factor  history  for  the  E-Glass 
Epoxy  and  Boron  Epoxy  composites,  respectively.  In  both  of  cases  the  fibers  are  oriented 
along  the  x-axis.  Figure  7  and  8(a)  are  very  similar  to  Figure  6(a).  Figure  8(b)  shows  the 
singularity  at  t  =  a/cR  more  clearly.  The  singularity  can  be  quite  difficult  to  capture  using 
the  method  of  Rubio-Gonzalez  and  Mason  [2]  because  the  shear  and  Rayleigh  wave 
speeds  are  nearly  equal.  Here,  equation  (53)  gives  Kn(t )  explicitly. 

7  Conclusions 

(1)  For  impact  concentrated  shear  loads  on  a  semi-infinite  crack  in  an  orthotropic 
material,  the  arrival  of  the  dilational  wave  at  the  crack  tip  makes  an  important 
contribution  to  the  dynamic  stress  intensity  factor  Kn  (f ) ,  unlike  in  mode  I  loading  where 

dilatational  wave  had  a  minimum  effect  on  Kx  (t) . 

(2)  The  long  time  value,  Ku(°°)  =  -qjlfnu ,  is  independent  of  the  mechanical 
properties  of  the  orthotropic  material,  however  the  curve  shape  and  magnitude  of  Kn(t) 
from  t  =  a/cd  to  t  =  a/cR  is  very  dependent  on  the  mechanical  properties  of  materials. 

(3)  From  the  results  given  here  for  an  isotropic  material  and  orthotropic  materials,  we 
can  conclude  that  the  normalized  K„(t)  from  t  =  a/ cd  to  t  =  a/cR  is  less  than  1  except 
at  a  singularity  at  t  =  a/cR .  In  addition,  under  similar  conditions,  when  fibers  are  parallel 
to  the  crack  direction  K„  (t)  is  generally  greater  than  when  the  fibers  are  perpendicular  to 
the  crack  direction. 
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(4)  Finally,  the  method  outlined  here  can  easily  be  extended  to  solve  the  problem  of  a 
semi-infinite  crack  under  more  general  impact  shear  loads  on  the  faces  of  a  semi-infinite 
crack.  Alternatively,  the  solution  presented  here  can  be  used  as  a  Green’s  function  to 
solve  more  general  problems.  Both  methods  should  yield  the  same  results. 
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Table  1  Mechanical  properties  used  for  the  analysis 


Isotropic 

material 

Graphite 

Epoxy 

E-Glass 

Epoxy 

Boron 

Epoxy 

Et  ( Gpa ) 

200 

156.75 

45 

207 

E2  (Gpa) 

200 

10.41 

12 

19 

vl2 

0.3 

0.31 

0.19 

0.21 

^23 

0.3 

0.49 

0.19 

0.21 

Mn  (°Pa) 

76.92 

7.07 

5.5 

6.4 

p(Kg!  m3) 

7840 

1580 

2100 

1990 

Figure  1  Schematic  of  the  semi-infinite  crack  under  concentrated  shear  impact  loads 


Figure  2  Branch  cuts  in  the  co-plane  Figure  3  Branch  cuts  in  the  ^-plane 
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1  -  1/c, ,  branch  point  1  ~  1/cd ,  branch  point 

2  -  1/c, ,  branch  point  2  -  1/c, ,  branch  point 

3  -  1/c„ ,  pole  point  3  -  1/cB ,  pole  point 

4  --  t/a  ,  branch  point  4  --  t/a  ,  branch  point 

(a) For  l/cs  <tla<\lcR  (b)Fort/a>l/ 

Figure  4  Branch  cuts  in  the  ^-plane  and  Integral  path 


dilatational  wave 


Figure  6(a)  Stress  intensity  factor  history  in  Graphite  Epoxy  under  concentrated  shear  loads 


Rayleigh  wave 


(t)/K  ,H 


Figure  6(b)  Stress  intensity  factor  history  in  Graphite  Epoxy  under  concentrated  shear  loads 


dilatational  wave 


Plane  strain 


Figure  7  Stress  intensity  factor  history  in  E-Glass  Epoxy  under  concentrated  shear  loads 


H  » /  (i) 


Figure  8(a)  Stress  intensity  factor  history  in  Boron  Epoxy  under  concentrated  shear  loads 


Figure  8(b)  Stress  intensity  factor  history  in  Boron  Epoxy  under  concentrated  shear  loads 


Chapter  5 


The  Dynamic  Stress  Intensity  Factor 
and  Strain  Energy  Release  Rate  for  a 
Semi-infinite  Crack  in  Rotated 
Transversely  Isotropic  Materials  due 
to  Uniform  Impact  Loading 


co-authored  with  C.  Y.  Wang  and  submitted  to  International  Journal  of  Fracture 

The  transient  elastodynamic  response  of  a  rotated  transversely  isotropic  material  under 
uniform  normal  impact  loading  on  the  faces  of  a  semi-infinite  crack  is  examined.  Three 
loading  modes  are  considered,  i.e.,  opening,  in-plane  shear  and  anti- plane  shear.  Their 
solutions  for  the  stress  intensity  factor  history  around  the  crack  tip  are  respectively  found. 
Laplace  and  Fourier  transforms  together  along  with  the  Wiener-Hopf  technique  are  employed 
to  solve  the  equations  of  motion.  The  analyzed  asymptotic  expression  of  stress  near  the  crack 
tip  leads  to  a  closed-form  solution  of  the  dynamic  stress  intensity  factor  for  each  loading 
mode.  It  is  found  that  the  stress  intensity  factors  are  proportional  to  the  square  root  of 
time  as  in  the  isotropic  case.  Results  for  rotated  transversely  isotropic  materials  converge 
to  known  solutions  for  un-rotated  orthotropic  materials,  un-rotated  transversely  isotropic 
materials  and  isotropic  materials  as  special  cases.  For  shear  loading  on  a  penny  crack  in  a 
transversely  isotropic  material,  the  analysis  is  used  to  find  approximate  strain  energy  release 
rates. 


76 


1.  Introduction 


There  are  many  closed-form  solutions  for  stationary  and  propagating  cracks  in  isotropic 
materials  under  dynamic  loading  [1,  2].  The  stationary  semi-infinite  crack  under  uniform  step 
loading  in  the  crack  faces  was  first  considered  by  Maue  [3].  The  equivalent  propagating  case  was 
studied  by  Baker  [4].  Many  solutions  for  various  loadings  of  cracks  in  orthotropic  materials  have 
been  found  by  applying  transforms  to  the  displacement  formulation  of  the  equations  of  motion. 
Previous  researchers  [5,  6]  following  this  approach  had  solved  the  resulting  dual  integral 
equations  using  the  method  of  Sneddon  [7,  8].  For  a  semi-infinite  crack  under  concentrated  or 
uniform  loads  in  orthotropic  materials,  Rubio-Gonzalez,  Mason  and  Wang  [9,  10  and  11]  solve 
the  same  equations  by  converting  the  dual  integral  equations  to  a  Wiener-Hopf  equation. 
Solution  of  that  equation  requires  only  a  straightforward  application  of  the  method  of  Noble  [12]. 
The  Wiener-Hopf  technique  was  used  to  solve  the  equation,  find  the  stress  ahead  of  the  crack  tip 
and  find  the  displacement  on  the  crack  faces.  The  asymptotic  expression  of  stress  near  the  crack 
tip  leads  to  the  stress  intensity  factor  K(t) . 

The  growing  use  of  composites  in  many  engineering  applications  demands  a  fundamental 
understanding  of  the  response  of  cracked  orthotropic  bodies  to  impact  loading.  In  this  paper,  the 
problem  of  uniform  normal  impact  loading  on  the  faces  of  a  semi-infinite  crack  in  a  rotated 
transversely  isotropic  material  is  analyzed  using  the  same  method  as  [9,  10,  and  11].  While  it  has 
been  thought  that  the  Wiener-Hopf  technique  could  not  be  applied  to  problems  such  as  this  [1, 
13],  the  approach  has  recently  been  shown  to  be  valid  in  such  cases.  Through  Laplace  and 
Fourier  transforms  combined  with  the  Wiener-Hopf  technique,  a  closed-form  solution  for  the 
stress  ahead  of  the  crack  tip  is  sought.  The  asymptotic  expression  of  stress  near  the  crack  tip 
leads  to  the  stress  intensity  factor  K(t) . 
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The  analysis  is  motivated  by  investigations  of  dynamically  loaded  penny  shape  cracks  in 
orthotropic  (or  composite)  materials.  Because  of  the  orientation  of  the  composite  fibers,  or 
principle  axes  of  the  orthotropic  material,  the  penny  shaped  crack  problem  is  not  axisymmetric 
and  can  not  be  solved  using  an  axisymmetric  approach.  As  an  approximation,  in  this  work, 
attention  is  focused  on  regions  very  near  the  crack  front  where  a  plane  strain  approach  becomes 
an  accurate  approximation.  There,  the  problem  can  be  reduced  to  that  of  a  semi-infinite  crack  in 
a  rotated  orthotropic  material  for  short  times.  This  is  only  approximate  because  the  effects  of 
waves  generated  at  other  points  on  the  crack  front  are  neglected.  Thus  the  exact  solutions 
presented  here  for  semi-infinite  cracks  offer  only  an  approximate  solution  for  the  penny  shaped 
crack.  The  solution  is  useful,  however,  in  developing  intuition  about  the  most  probable  location 
of  crack  initiation  for  the  penny  shaped  crack  and,  consequently,  can  give  qualitative  predictions 
of  the  final  shape  of  the  crack  after  finite  duration  loading. 

2  Equations  of  motion 

The  quasi-three-dimensional  problem  of  an  infinite  transversely  isotropic  body  with  a 
semi-infinite  crack  is  considered  in  Figure  1.  The  body  is  initially  stress  free  and  at  rest.  A 
Cartesian  (x,,x2,x3)  coordinate  system  is  chosen  to  coincide  with  the  principle  axes  of  material. 

The  transversely  isotropic  material  has  a  symmetry  plane  on  x3  =  0  and  an  axis  of  symmetry 
along  x3.  Another  right-handed  rectangular  (x,  y,  z )  coordinate  system  is  introduced  in  the  body, 
oriented  so  that  the  z-axis  coincides  with  the  crack  edge.  The  y-axis  is  identical  with  the  x2  -axis, 
consequently,  the  coordinate  system  (x,  y,  z)  is  related  to  the  ( x, ,  x2 ,  x3 )  coordinate  system  by  a 
rotation  with  0  around  the  x2  -axis.  The  crack  faces  are  suddenly  loaded  by  uniform  traction,  as 
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shown  in  Figure  1.  Let  Ejt  G(>.  and  vv  (i,  j  =  1,  2,  3  )  be  the  engineering  elastic  constants  of  the 

* 

material  where  in  the  indices  (1,2,  and  3)  respectively  stands  for  the  directions  (  jc,,  x2,  x3 ). 

The  problem  is  assumed  to  be  two-dimensional  and  to  include  inertial  effects.  However, 
because  of  the  constitutive  law,  the  out  of  plane  displacement  (in  the  z-direction)  can  not  be 
neglected  and  a  coupling  exists  between  in-plane  and  out-of-plane  deformation.  Hence,  three 
components  (u,  v  and  w)  of  the  displacement  vector  exist  in  the  (x,  y,  z)  coordinate  system  and 
are  independent  of  z,  that  is, 
u  =  u(x,y), 

v  =  v(x,y),  (2.1) 

w  =  w(x,y). 


In  the  (x,,jc2,x3)  coordinate  system,  the  non-dimensional 


transversely  isotropic  material  is 


C  = 


'11  M2 


'12  ''ll 


'13  M3 


0 

0 

0 


0 

0 

0 


M3 

Cj3 

c33 

0 

0 

0 


0 

0 

0 

C55 

0 

0 


0 

0 

0 

0 

C55 

0 


0 

0 

0 

0 

0 

C11  ~C12 


(2.2) 


where  C  has  only  five  independent  constants.  The  non-dimensional  parameters  c..  have  the 
following  relations  with  the  engineering  elastic  constants  for  plane  strain  conditions 


'ii 


_Ej(l  Vi3V31  )  ^(1  V12)  _  £/V12  +  V13V31  l  _  fiV13 

GI2(l  +  v12)A’  33  G12A  ’  12  G12(1  +  v12)A  13  G12A 


Gn  Ex 

c55  “  7^  9  A  =  (1  —  vn  —  2v13v31  ) ,  Gn  — 


M2 


2fl  +  v12J  E3  Ex 


Zll_Z3L  £n_£n=1  (23) 
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Next,  by  a  coordinate  transformation  involving  a  rotation  about  x2 ,  the  non-dimensional 
stiffness  matrix  C  in  the  ( xt,x2,x2 )  coordinate  system  is  changed  to  a  new  non-dimensional 
stiffness  matrix  D  in  the  (x,  y,  z )  coordinate  system  [14], 


D  = 


Ml 


dn 

d\ 3 
0 


d 


15 

0 


12 

i22 

^23 

0 

^25 

0 


“l3 

^23 

^33- 

0 

d35 

0 


0 

0 

0 


d 


44 

0 


d 


46 


“15 

^25 

^35 

0 

d55 

0 


0 

0 

0 

^46 

0 

d&6 


where  the  thirteen  different  constants  in  D  are  given  as 
du  =cu  cos4  9  +  2fc13  +  2c55  )cos2  Osin2  9  +  c33  sin4  0 , 


(2.4) 


dn  =c12  cos20  +  c13sin20. 


dn  =ci3  cos*6+(cn  +c33  -4c55)cos20sin20  +  c13sin40  , 
d,5  =  (c„  - c,3  - 2c55 )cos3  0sin0+(c,3  -c33  +2c55  )cos#sin3  9 , 


d  ~  c 

U22  Cll’ 

t?23  =  C|2  sin2  9  +  c13  cos2  9 , 
dis  ~(cn -c13)sin0cos$, 

d33  =  c„  sin4  0  +  2(ci3  +  2cS5  )cos2  0sin2  9  +  c33  cos4  9 , 

d35  =(c13  -c33  +2c55)cos3  ^sin^+fc,,  -c,3  -2c55)cos0sin30, 

du  =—(cn  -c12)sin2  9  +  c55cos2  9 , 

2 

=“(cn  —  ^12  —  2cj5  )cos  9  sin  9 , 
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d55  =  c55cos40  +  (cu  -2c13  +  c33-2c55)cos20sin20  +  cS5 sin40, 

t 

^66  =='T’f^ll  Cj2  ^cos  #  +  c55  sin 


Note  that  for  the  case  0  =  0  the  material  is  isotropic  in  the  plane. 


(2.5) 


(2.6) 


du  —  cu,dl2  —  cn*dn  —  c13,  <i,  5  —  0,  d22  —  Cn,<f23  —  c13,^25 

^33  =  C33>  ^35  =  ^’^44  =  C55»^46  =  ®»  ^55  =  C55  ’  ^66  =  (Cll  —  C12  )/^- 

For  0  =  90°  the  material  is  orthotropic  in  the  plane, 

d\l  =  c33>  d\2  =  c13»  ^13  =  C13>  ^15  =  ^22  =  C\l>  ^23  =  C!2'  ^25  — 

^33  =  C\V  d 35  =  0,  C?44  =  (C,,  —  C|2  )/2,  =  0,  </S5  =  CJ5,  ^gg  —  CJ5. 

Due  to  the  two-dimensional  assumption,  equation  (2.1),  all  the  derivatives  with  respect  to 
are  zero.  Thus,  in  the  (x,  y,  z)  coordinate  system  the  equations  of  motion  [16]  become 
,  <?2u  ,  <?2u  .  ,  d2\  ,  <?2w  ,  <?2w  _  «?2u 

d"^  +  d“^  +  (d“+dn)~&ty+dl5  /5r2  ^  /9v2  _.2^2’ 


d2u 


dy1  c2dt 


d2v  ,  d2\  ,  ,  .  ,  <?2w  <?2- 


(^66  +  ^12  )  -3-5-  +  ^66  3—7  +  ^22  TT  +  (^46  +  ^25  ) 

c*cey  «?  x  ey  i/Aw/  Cj  c/j 

,  <?2u  .  J  <?2u  .  /J  ,  J  N  <?2v  ,  J  <?2w  t  J  <92w  _  <?2w 

^15“^2"  +  rf46  ^y2  +^46  +«25)  ^^.+rf55  ^2+rf44  ^,2  .2^,2 


(2.7) 


<&<?y  55  A:2  44  <7y2  cs2dt7 

where  u,  v  and  w  are  the  x,  y  and  z  components  of  the  displacement  vector.  In  the  rotated 
transversely  isotropic  solids  shown  in  Figure  1,  cs  =  ^Gn/p  represents  the  velocity  of  the  shear 


wave  propagating  in  xxx2 -plane  (cd  =  <Jc^cs  is  the  velocity  of  the  dilatational  wave  propagating 
along  x2  -direction  in  x,x2 -plane  which  will  be  used  later)  and  p  is  the  mass  density.  When 
6  =  0  or  6  =  90° ,  dl5  =d25  =  d^  =  0  and  the  third  equation  becomes  uncoupled  from  the  first 

two  equations,  when  O<0<9(? ,  this  coupled  system  is  required. 

Finally,  the  stresses  are  related  to  the  displacements  by  the  equations, 
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xy  _ 


G12 

Gy 


=  d 


dw 

dy 


,  ,di  dv. 
+  d6 «(-r  +  -T-). 


'12 
G|2 


X  _ 


Tn 


dy  dx 

.  dw  ,  ,dx  dv^ 

d*4  37~  ^46  (37 + 37^  ’ 

ay  dy  ox 

,  da  ,  dv  dw 

:  “12  37  +  "22  3  *"  "25  37"  * 

a*  ay  ox 

,  dx  dv  dw 

d\  \  —  +  «i2  37  +  “15  33 » 
ox  dy  .  dx 


,  dx  , 

—  U]i  _  +  w 


<?v 


+  (l 


dw 


'12 


Vxz  _ 


'12 


M3  a  1  **23  ->  1  “35  > 

dx  dy  dx 

,  di  ,  dv  ,  dw 

d\5  ->  +  “25  “5T  **55  “TT  • 

ox  dy  dx 


(2.8) 


3  Normal  Impact 


For  the  loading  shown  Figure  1,  the  spatially  uniform  normal  traction  with  magnitude  o0  is 


suddenly  applied  on  the  crack  faces.  Exploiting  symmetry  and  limiting  ourselves  to  upper  half¬ 
plane,  y>0,  the  corresponding  boundary  conditions  are 


oy(x,0,t)  =  -o0H(t), 
rxy(x,  0,t)  =  0 
t^Cx,  0,  t)  =  0 
v(x,  0,  t)  =  0 , 


-°°<;c<0, 

-  00  <  X  <  +00  , 

-00  <x<  +00  , 
0  <  x  <  +00 , 


(3.1) 


where  H(t)  is  the  unit  step  function. 

The  displacement  at  infinity  is  zero,  and  the  body  is  stress  free  and  at  rest  everywhere  for 


r<0. 


3.1  Laplace  transform 

In  equations  (8),  the  time  variable  may  be  removed  by  application  of  the  Laplace  transform 
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f\s)  =  ^me-dt,  m=^-lr(s)e"ds=^Tf;y(s)e-ds,  (3.1.1) 

where  Br  denotes  the  Bromwich  path  of  integration  which  is  a  line  parallel  to  the  imaginary  axis 
in  the  s-plane.  Thus,  using  the  Laplace  transform  representation,  the  displacements  have  the  form 


u(*>  y>  t) =  2xi , 

u*(x,  y,  5)  e”ds , 

*  Br 

v’(x,  y,  s) 

*  Br 

(3.1.2) 

,  „  1 
w(x,  y,t)=  . 

2k  i 

I  w*(x,  y,  j)  . 

WBr 

Applying  the  relations  (3.1.2)  to  equations  (2.7)  and  assuming  zero  initial  condition,  the 
transformed  domain  equations  become 


J  <?2U*  J  <?2U*  J  X<?2V*  <?2W*  <?2w*  J 

d"^+d«WH  “  a)^  '1~&r  — 

.  .  <?2u*  .  <?V  .  <?2v*  .  .  ,  .<?2w*  52v* 


/J  J  n^2u*  j  ^2v*  j  ^2y*  /  j  j  n^2W*  52V*  _a 

,2>^  “7^+d22l^+(  *  ^  c/  ’ 

j  <?2U*  J  <?V  ,J  j  x<?2v*  J  <?2w*  J  <?2w*  s2 w*  _A 

d,5“*r+  4677'+C  “  B  A2  “  A2  C.2  ■ 


(3.1.3) 


Applying  the  relations  (3.1.2)  to  equations  (2.8),  the  transformed  relations  between  stresses 


and  displacements  become 
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'O'  _ 


=  d. 


■'12 


<9w* 


+  ^6fi( 


<?y  dx 


T>Z  j  ^W*  j  >*  oV*. 


G  ,  da*  j  dv*  dw* 

—  -  a,  2  — —  +  a 22  —  +  a  25  -z— , 

cx  oy  <zk 


'12 


<7X*  .  ob*  .  oV*  ,  o\v* 

—  1 1  ~  1"  Uo  r  \~  di , 


(3.1.4) 


'12 


Ml  ->  *  “12  n  ^“15  n  * 

ax  ay  ax 


<J  *  _  .  oh*  j  dv*  j  dw* 

"IT  B  55  IT’ 


'12 


G,2  "  15  &  “ 


23  <?y 
<V 


.  av* 
ay +4i  dc  • 


The  application  of  the  Laplace  transform  to  the  boundary  conditions  (3.1)  gives, 


oy  (x,  0  ,t)  =  -a Q/s , 
r^fx,  0,  0=0 , 
x^(x,  0,  0=0 , 
v*(x,  0,  0=0 , 


-°°<jt<0  , 

-°°<JC<+©o  , 
-°°<JC<+<»  , 
0<X<+oo  . 


(3.1.5) 


3.2  Fourier  transform 

To  obtain  a  solution  of  the  differential  equations  (3.1.3)  subject  to  conditions  (3.1.5),  the 
Fourier  transform  is  applied, 

7(®)  =  JV(*)  eia*dx,  f(x)  =  ^-^f(co)e~imdcQ  (3.2.1) 

It  is  noted  that  the  boundary  conditions  (3.1.5)  are  defined  only  on  half  of  the  range  of  x. 
Consequently,  the  Fourier  transform  can  not  be  applied  to  these  boundary  conditions.  To  remedy 
this  situation,  the  boundary  conditions  are  extended  to  apply  on  the  full  range  of  x.  Two- 
unknown  functions  v_*(x,j)  and  o+‘(x,s)  are  introduced.  The  function  v.*(x,j)  is  defined  to 
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be  the  y-direction  displacement  of  the  crack  face,  y=0,  for  <  x  <  0  and  0  <  t  <  °° ,  and  to  be 

l 

identically  zero  for  0  <  x  <  +°°  and  0  <  t  <  °° .  Likewise,  the  function  o+\x,  s )  is  defined  to  be 
the  shear  stress  in  the  x-direction  on  the  plane,  y=0,  for  0<x<°°  and  0<r<°°,  and  to  be 
identically  zero  for  <  x  <  0  and  0  <  t  <  °° .  With  these  definitions,  the  boundary  conditions 
(3.1.5)  can  be  rewritten  as 


<j*(x,  0 ,  s)  =  -—H(  -x)+a+*(x,s) ,  -°o<x<+°°, 

*  r 


T,,.  (x,  0,  j)=0  , 
r)Z*(x,  0,  s)=0 , 
v*(x,  0,  s)=v_*(x,s) , 


-oo<X<+°°  , 
-°°<X<+00  , 


-oo<X<+oo  . 


(3.2.2) 


Thus,  the  Laplace  transformed  displacements  have  the  form 
u(x,  y,  s)  =  — —  f  u*(tf>,  y,  s)  e~laxdo)  , 

2  K  J-o 

v*(x,  y,  s)  =  — —  f  v*(ty,  y,  s)  e~io,xdco  , 

2tc  J~ 

w‘(x,  y,  s) =  —  [ w* (o),  y,  s)  e~'axd(D  . 

2  n  J— 


(3.2.3) 


Substituting  these  transforms  into  equations  (3.1.3),  the  function  n*,  v*  and  w*  are  found 
to  satisfy  the  ordinary  differential  equations 

-(duco2+s/cs  )u  +^66^— 7 -(^66 +^,2)6* -3 — <*15^  w  +d46— ^-  =  0, 

d  y  dy  dy 


-  (<*66  +  <*12  - (<*66^2  +  S2/ Cs  ) V*  +  <*22  - (<*46  +  <*25  )<**' =  0  »  (3'2'4) 


<?y 

<?2s* 


-  dlsO)2U  +  d«  -  (<*46  +  <*25 - (<*55^2  +S2/c2)w*  +d«  =  0  . 

d  y  dy  dy 


This  system  of  equations  can  be  rewritten  as 


(3.2.5) 


(d66  0  f 
0  d22  0  1<?  u 

\d46  0  ^ 

r-du(Dl -s2/cs 


& 


0 

—  (^66  dtt 


’(^66  ~^dl2)ox 


0 


0  (^46  "t”  ^25  )® 

0  ~(rf46+  </„)tfX  0 


0 


-d„(0 


46  '  “25 
2  \ 


<?0* 


d> 


0  ~di6(02  -s2/c2  0 

-dl5o2  0  -dS}6)2  —s2/c2 


P*=0, 


where  U*  =  (u\  v\  w’)T .  This  is  a  second  order  system  of  ordinary  differential  equations  with 
constant  coefficients  having  solutions  of  the  form 

U*  =  A (co,  s)e~n ,  y  >  0 .  (3.2.6) 

where  A(ty,  s)  =  (A, ,  A2 ,  A3  )r 

The  solution  of  equation  (3.2.6)  under  the  condition  of  zero  displacement  at  infinity  is 

3 

u*(<y,  y,s)  =  ^jAi(o),s)e~r,y , 

l=I 
3 

\\co,  y,s)  =  Y,  ioyYiai  (<».  s)A,  (co,  s)e~r,y , 

1=1 

w‘(<y,y^)  =  y^uPi(co,s)Ai(co,s)e~r,y . 


(3.2.7) 


where  A^ra.s)  is  arbitrary  function,  and  a^co^),  yQ,  (<u,5)  stand  for  the  functions 

CC- (co  s )  =  ^46  25X^66 Yi  ~dn®  ~s  /cj  )~(d^ +d12)(c/46/l-  —  d15fi)  ) 

WaTi*  “<*6«®2  -^/c^id^y2  -dl5co2)  +  (d66  +dn)(d46  +d25)co2y2  ’ 

,  _  (<*66  +<f12)2<0V|2  +(<*22ft2  ~^66^2  -^Vc/X^y,-2  -J2/c/) 


(3.2.8) 


Pi(co,s)  =  - 


(d22 Yi2  ~d66C°2  ~^lcs^d^7i  ~dl5(°2)  +  (d66  +  d\2^d46  +d2 ' 


with  y,2  ,  y22  and  y32  being  the  eigenvalues  of  the  sixth  order  equation 


det 


d66Y2~dn(i>2-s2lc1 2  W66+^i2)^y  d«Y2  ~d\%<t>2 

(d(t  +dn)aiy  d22y2  -d^co2 -s2 fc2  (dA6+d25)coiy 

d4eY2 -dlsco2  (dA6  +d25)coiy  dHy2 -dssco2 -s2  /c,: 


=  0. 


(3.2.9) 


This  determinant  is  further  expanded  as 
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(3.2.10) 


6  b  4  c  2  d  _ 
a  a  a  ' 

where 

^22^46  *”^66^44)* 

=  (^22^44  —  ^46  ^^22^66^^44^66^  lCs  "^(”^12  ^44^^11^22^44 

—  Zd^dyid^  +  2<f126?25^46  ^  ^\2^46  “”^25  ^66  ”^12^44^66  “**  ^22^55^66^^  * 

C  =  “[(^22  **“  ^44  ^66 )  ^  /C5  ^”(“^12  J(~^'\\^22  ^25  ~^^\\^44~~ 

^15^46  “^25^46  *^46  “*"^22^55  ““^12^66  +^44^66  d^d^CQ  S  jc$ 

+  (— </j5  ^22  2^12^15^25  ^1 1^25  +  ^\2^\5^46  ~~  ^1 1^25^46  ^11^46 

^12  ^55  ^11^22^55  ^15^25^66  ^11^44^66  “  ^d\2^5$d^)CO  ], 

d  =  (d66G)2+s2/c2)[s4/c4  +(du  +d55)a)2s2/c2  +(dud55 -dl52)0)4]. 

From  (3.2.10),  it  is  readily  found  that  the  functions  y,(to) ,  y2(a>)  and  y3(a>)  are  multiple 
valued  functions  of  tu  in  the  complex  6)  -plane  with  branch  point  at  o  =  ,  <y  =  ±Hs/ks2 

and  <y  =  ±1  s/k^  ,  respectively,  which  are  the  roots  of  the  equation  given  by  d=0,  i.e. 

(^O)2  +  s2/c2)[s4/c ,*  +(dn  +d55)co2s2/cs  +(dnd55 -dx2)co4}  =  0,  (3.2.11) 

so  *c„  =  ^c,, ,  Ki2  =  ,  and  Kj  =  Jd^cs[ , 

where 

_  2(dxldS5  —dl5  ) 

Uj  —  - -  , 

1  +  ^55  +  •y 1  +  4^?,5  —2  dud55  +  d55 
j  _  2(dlld55—dl5  ) 

«2  —  1 - - 

^11  ^"^55  “  V^H  ^"^15  “2^11^55  +^55 

Note  that  for  the  case  0  =  0,  K^,0  =-y[c^cs,  ks2  =  J  - •  - C^-c5  =  cs  and  /r/  =  <{c^cs .  The 
parameters,  r^0 ,  ksZ°  and  Kd° ,  are  the  shear  wave  speed  in  jc2x3 -plane,  the  shear  wave  speed  in 
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xtx2 -plane  and  the  dilatational  wave  speed  along  *2  -direction  in  x,x2 -plane,  respectively;  for 

t 

the  case  6  =  90%  k™  =  f~cs ,  ks2%  =  <Jc^cs  and  a c/°  =  Jc^cs . 

Later  we  will  be  interested  in  the  behavior  of  these  functions  for  CO  ->  °° .  In  this  case, 
y.2  ->  6)2N2  (i  =  1, 2, 3) ,  where  N.1  are  the  roots  of  the  following  equation 

A  A 

^  +  *#<+£**+^.  =  6  (3.2.12) 

a  a  a 


where 


a  ~  d72{dA6  ^66^44)* 

=  ”^12  ^44  ^11^22^44  “  ^^15^22^46  ^12^25^4 6 

2^12^46  ^25  ^66  ^12^44^66  ^  ^22^55^66  » 

c  =  —(—d\*d21  +  2dndl5d25  —  rfn^25  +  ^dndxsd^  —  Id^^d^  —  dud^ 

—  rfj2  ^55  "^^11^22^55  ^"^^15^25^66  ^”^11^44^66  ”^^12^55^66)’ 

a  2 

^  =  d^{dnd55  —  d15  ), 


Similarly,  asta  — » <*>,  af(£a,s)  and  y5,  (<M,s)  become 


a.(fi>,s)->ar~ 
fii(co,s)— »/?, 


(4«  +  4^(4^, 2  -dn)-{d«  +  4I2)(4JY,2  -415)  1 

(4**,* -dJid^N*  -dJ+V*  +412X<*46  +425)Ni-2  <»2  * 
(4«  +412)2//,.2  +(4^, 2  -‘du)(dJ2Ni2  -4«) 

(422N,2  -466)(446Ar,2  -415)  +  (466  +di2)(d46+d25)Ni2 ' 


(3.2.13) 


It  can  be  shown  that  for  many  materials  that  the  roots  of  (3.2.10),  }  ,,  }  2  and  )  3 ,  are  real 


and  positive  and  the  expressions  for  the  displacements  in  the  Laplace  transform  domain  become 
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u  =  —  f  ( V'm .+  +  AJe"w  )e'lwdco , 

2k  J^o 

v*  =  —  f  {yxaxAxeny  +  r2<*  2  Vm  +  y^A^^e^dco , 
2k  J-~ 

w’  =  —  f  (£,  A.*'m  +  p2A2e'™  +  p3A3e1,yy)eimd(0 . 

2n  J— 


(3.2.14) 


Substituting  these  expressions  for  displacements,  (3.2.14),  into  the  constitutive  law,  (3.1.4), 


we  have 


—¥—  =  — —  f” (m,A,e  r,y  +  m2A2e_,'J,’  +  m3A3e  r3>’)e  ,a*d(0, 
Gn  2k  J-»  1  33 

— = — —  r(n,A1e-r,y +/i2A2c-,'j5'  +  n3A3e~r'y )  e~,a*  do) , 
G„  2;r  J-TC 


(3.2.15) 


(kxAxe~r'y  +  k^e-™  +  k3A3e~ny )  ei0*d(d , 


where 


m(.  =(^46  -d^o)2^.  +d4APi)yi, 
ni  =(d66  ~d(Aoi1ai  +d46Pi)yi, 


kj  —  dX2  +  d^i  o-i  +  d^Pi . 


Substitution  of  expression  forT^,T*z,  the  first  and  second  of  equations  (3.2.15),  into  the 
boundary  condition  for  x*  ,  x*  ,  the  second  and  third  of  equations  (3.1.5),  yields 


m,  A,  +  m2A2  +m3A3  =0, 
nx  A,  +  n2A2  +  tt3A3  =  0. 

These  equations  can  be  rewritten  as 

A1=Xx  ( a >,  s) A, ,  A3  =  A2  ( a ),  s) Ax , 


(3.2.16) 


;t,(&u)=  1  3  3i,a2 

m3n2  -m2n3 


.  mxn2  -m2nx 

,X2(0),s)  = - i-2. - =— L 

m3n2  -tn2n3 


(3.2.17) 


As  co  — >  ,  we  have 
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(3.2.18) 


Aj(q},s)  — >  \ 


m~n~  - m~n{’ 
m~n~  - m~n. 


m~n2” 

m~n~  -m~n2 


where  m^oi.s)  — > com"  and  n^m.s)  — » <uw(~  as  ra  — » oo . 

Using  this  result  in  the  expression  for  normal  stress,  the  third  of  equation  (3.2.15),  the 
expression  for  vertical  displacement,  the  second  of  equations  (3.2.14),  and  the  remaining 
boundary  conditions,  we  obtain  " 


(k,  +  A, k2  +A2k3)coAlie  ,OKdco=  -<t0/sH(-x)  +  <t+*(x,5), 

2k  J~°° 

1  f00 

— J  (ai7i  +^2/2^  =  vl(x,5). 

in 

By  Fourier  transform  inversion,  these  equations  become 

-GI2(k,  +  X,,k2  +X,2k3)coA1/  =  -a0/i(os  +  E+(oo,s), 

(a,y,  +  a2y2X,  +  a3y3X,2)a>A,/  =  V_(om). 


(3.2.19) 


(3.2.20) 


where 


X+(<y,j)  = 
V_(cd,s)  = 


f  <J+\x,s)e'mdx , 

Jo 

J  v {x,s)eim dx  . 


(3.2.21) 


Eliminating  A,  from  equation,  (3.2.20),  we  obtain  a  Wiener-Hopf  equation, 

GnXx  =2==  V  (0),  s)  =  a  Jicos  -  E+  (<y,  5) ,  (3.2.22) 

ja)2  +  s2/icd2 


where  the  following  functions  are  introduced 

^  +  5  ! Kd  (^1  +  ^7^3  ) 

+«2M  +a3r34) 

(3.2.23) 

_  *I-+v*2"+w 

^  1  W  ir  OO  ..  ^  oo  00  -  T  /I  OO  5 

OTj  N  2^\  +  <*  3 

(3.2.24) 
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and  X\  ls  dependent  of  6  and  the  material  parameters.  Its  complete  expression  is  given  in  the 

t 

Appendix. 

Here,  note  that  the  case  of  6  =  0 ,  Xi°  =  — ^ — —  (it  is  a  limit  value),  and  6  =  90° , 


90  _  (C|iC33  C|3  )[C33+C„A,^2Ar3(Ar,+7V2+^3)3 

"  cI1c33(A1+iy2)(iV3+^2)(^3+iV1) 


(3.2.25) 


where  A,,  are  shown  in  equation  A.5  of  the  Appendix.  X\  and  X\°  are  the  limit  values  of  X\  at 


6  =  0  and  6  =  90° ,  respectively. 

For  a>  =  0,  after  lengthy  algebraic  manipulation  of  equation  (3.2.23),  we  have 


F(0,s)  = 


■y2  -fin 

i  c. 


*dX\  cd  ' 


(3.2.26) 


From  the  physics  of  the  problem,  it  is  reasonable  to  assume  that  the  function  o+'(x,5)  and 
v_’(x,s)  are  exponentially  bounded  at  infinity,  which  ensures  the  existence  of  their  Fourier 

transform  (3.2.1).  In  particular,  it  is  shown  by  Noble  [12]  that  if  |<t+*(x,j)|  < MxeLx  as 

then  L+(ta) is  analytic  in  Re(ffl)  =£>£_,  and  if  v_*(x,s)|  <M2e^x  as  x  ->  -<»  then  V.(<y)  is 
analytic  in Re(o)  =  £  <£+. 

This  function,  F(a,s) ,  is  analytic  everywhere  in  the  complex  plane  except  on  at  the  branch 
cuts  of  y  i  at  co  =  ±i  s/ks1  ,  co  =  ±i  s/ks2  and  co  =  ±is/Kd  .  It  is  single  valued  in  the  <a  -plane  cut 
as  shown  in  Figure  2.  For  general  anisotropic  materials,  the  Rayleigh  wave  equations  are  given 
by  Nayfeh  [15].  We  found  the  algebraic  equations,  F(tfM)=0  (letting  a  =  is/v)  and  equation 
(3.2.10)  for  }  are  the  same  as  the  Rayleigh  wave  equations  because  F(co,s)=0  is  equivalent  to 


tyz=txy=Oy  =0  which  are  the  necessary  boundary  conditions  used  when  deriving  the 
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Rayleigh  wave  equations.  The  system  of  Rayleigh  wave  equations  is  dependent  on  the  "material 

( 

properties"  referring  to  the  transformed  properties  through  the  rotation  of  the  azimuthal  angle  6 . 
Therefore,  the  Rayleigh  wave  speed  cR  is  varying  with  the  azimuthal  angle  6 .  In  general, 
finding  the  critical  value  of  cR  requires  numerical  calculation  because  of  the  complexity  of  the 
system  of  Rayleigh  wave  equations  for  general  anisotropic  materials  [16].  It  is  true,  however, 
that  the  only  zeros  of  F(co,s)  are  of  the  form  ±is/cR  . 


3.3  Wiener-Hopf  technique 

The  Wiener-Hopf  technique  can  be  outlined  as  follows,  the  function  L(co)  is  defined  and 
factored  as 


Um)^  =  GaX>  F(^> 


M«» 

thus,  the  Wiener-Hopf  equation,  (3.3.1),  becomes 


L  (<y)V  (o))  =  T±L+(a))-?,+((o,s)L+(co) . 

ICQS 


If  the  function  D(a>)  can  be  decomposed  as 


D(o))  =  ^L+(o))  =  ^-\^ 
teas  s  L 


(a)-L+(0)  t  LJO) 


ICO 


1(0 


=  D+(0))  +  D_(co)  , 


with 


D+(o))  =  ^-\L^  L+(Q)1,  D  = 

s  L  io)  J  s  L  ico 

Then  the  Wiener-Hopf  equation,  (47),  is  further  reduced  to 

D+  (a )  -  Z+  (<a)L+  (co)  =  L_  (a) V_  (co )  -  D_  (co)  =  Q(ta) . 


(3.3.1) 


(3.3.2) 


(3.3.3) 


(3.3.4) 
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The  first  member  of  (3.3.4)  is  analytic  in  the  right  half  plane  of  Re(ta)  >  and  the  second 

t 

member  in  the  left  half  plane  Re(ra)  <  £+ .  If  ><*+,  the  regions  of  analyticity  overlap  and  by 
invoking  analytic  continuation,  it  is  concluded  that  Q(co)  is  analytic  and  single-valued  in  the 
whole  ty-plane  shown  in  Figure  2.  Furthermore,  invoking  the  extended  Liouville  theorem,  it  can 
be  shown  [12]  that  if  Q(co)  is  bounded  and  entire  and  £2(ta)  —>  0  as  co  — » <» ,  then  Q(<n)  =0. 
Hence,  we  can  solve  for  the  transform  of  stress  ahead  of  the  crack  tip,  £+(<u) ,  and  displacement, 
V_(ftj) ,  behind, 


V.(fll)  = 


L+(co)  sico 
D_(co)  <x0 


1- 


M  0) 


L_(co)  sico 


L+(co) 

L.(! 0) 


L.m 


(3.3.5) 


Following  this  outline,  the  first  step  is  to  factor  L(co) ,  (3.3.1),  by  defining 


F{a»  = 


F(ta,5) 

co2  +s2/cf 


2  • 


(3.3.6) 


It  can  be  shown  thatF(o)  -^1  as  (the  constant  X\  in  (3.2.24)  was  chosen  to  make 

this  possible).  The  function  F{co)  is  regular  and  F{co)  *  0  in  the  co  -plane  cut  as  shown  in 
Figure  2,  the  only  singularities  are  the  branch  points  shared  with  }  . ,  (i  =  1,  2) ,  i.e.  co=±  is j xs] , 
co  =  ±is/ks2  and  co=±is/Kd  .  It  is  well  known  that  factorization  of  a  function  is  accomplished 
most  directly  for  functions  that  approach  unity  as  co  -4  °°  and  that  have  neither  zeros  and  nor 
poles  in  the  finite  plane.  Indeed,  F(co)  is  an  example  of  such  a  function.  Therefore,  using 
Cauchy  integral  formula,  it  can  be  shown  that  [1] 
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A  A 


F(co)  =  F+(d))F_(co) 


t(<»)  =  expjy-T J 

[2^Jj?±  z-co  J 


(3.3.7) 


where  R_  is  the  closed  counterclockwise  contour  enclosing  the  branch  points 
+is/nsl ,  +  is/ k s2  and  +  is/ Kd  ,  and  also  R+  is  the  closed  counterclockwise  contour  enclosing 
the  branch  points  -is/nsi,  -is/ks2  and  -is/Kd  . 


Returning  to  the  factorization  of  L(co) ,  we  have  now 


F+  (co)F-  (cd)(co + i  s[  cR  ){co  -  i  sj_  cR ) 
■yjco  +  is/jCjijQ)  —  is/fCj 


L(Q})  =  L+(£0)L_(Q))  =  GnX\ 


(3.3.8) 


Therefore, 


L.m = gi2x>  ?  '.^g)F-(<») . 

V a-is/Kt 

M»)- 

(ct)  +  is/cR)F+((0) 


(3.3.9) 


Note  that  as  F+(ry)->l,  F_(tf))->1,  L+  (ry)  =  ,  L_(ty)  =  GnX\J®  and 


L+(0)  = 


s^k/^F+{0)M 


,  and  using  equations  (3.3.5),  we  can  readily  obtain  2+(ty)  and  V_(a)) 


as  Ct)  — i.e. 


^0Cft 


S+(®)->  v - 

(si)A^}J^F+(0) 


(3.3.10) 


V_  (<2>) 


°0CR 


G12^,(5ift))^/r^F+(0) 


(3.3.11) 
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3.4  Stress  Intensity  Factor 

t 

To  find  the  stress  intensity  factor,  an  asymptotic  expression  for  the  stress  near  the  crack  tip 
is  sought.  A  well-known  result  relating  asymptotic  expressions  between  a  function  and  its 


Fourier  transform  (Abelian  theorems)  [12,  17]  is 


lim 4xo* {x, s')  =  lime  ,,r/4 

X->0  6)->  oo 


lim 

jr->0 


V_*(x) 

4~X 


=  lim  e'*3/4 


v_m 


(3.4.1) 


(3.4.2) 


Hence,  the  definition  of  the  stress  intensity  factor  gives 


KUs)=  lim  «J2jdc<j‘ (x,  s )  =  lim  e  '*/4V2tf)E+(G>)  =  ■ 

x — >0+  0) — >-K>° 


V2cr( 


0 


(3.4.3) 


Note  that  F+( 0)  =  F_(0)  =  ^F(Q)  consequently,  equation  (3.4.3)  can  be 


written  as 


*/*(*)  = 


V2<J0 


^2V^“Vf(  0,5) 


(3.4.5) 


52  C. 


Using  equation  (3.2.26),  F(0,$)  = - — ,  it  readily  shown  that 

**X\  c* 


K*,(s)  = 


V20-0 


,3/2 


(3.4.6) 


-11 


Applying  the  inverse  Laplace  transform  to  K*,(s),  the  stress  intensity  factor  K,{t)  can  be 


written  as  follows 


Kj(t)=  2ffQCj  V7  =  2a Q  Vf"  I 

2  Vl 

’ 

7T  C  *  „ 

V 

11 

(3.4.7) 
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where 


r  -  IKK 

’~U 

Thus,  the  stress  oy(x,0,t)  can  be  expressed  as 


<jy{x,0  ,t)  = 


K,{t) 

Jinx 


(3.4.8) 


(3.4.9) 


For  6  =  0,  i.e.  for  the  case  of  transversely  isotropic  material,  C,  reduces  to 


P  o  _  |2crf  2(cn  1) 

'  1  *  C,2 


(3.4.10) 


By  using  c„  =  ^  ^  ,  equation  (3.4.10)  further  reduces  to  the  well-known  result  for  isotropic 
materials  [13],  i.e. 


^  o  _  <Jcd(l-2v)/7t 

vr  - 

(1  — v) 


(3.4.11) 


For  6  =  90° ,  i.e.  for  the  case  of  orthotropic  materials,  C,  reduces  to 


r90_  f2c^190 

7  ~i*  cn  • 


(3.4.12) 


where  X\°  is  shown  in  (3.2.25).  Tables  2  and  3  show  that  the  numerical  values  of  C,90  for 
orthotropic  materials  are  same  as  those  obtained  by  using  formulation  in  [10],  hence,  equation 
(3.4.7)  also  reduces  to  the  result  for  orthotropic  materials  as  expected  for  6  =  90° .  Thus,  the 
value  of  C,°  and  C,90  can  be  accurately  calculated  by  using  (3.4.10)  and  (3.4.12),  or  can  be 

evaluated  by  using  (3.4.8)  and  letting  6  =  e  and  6  =  90°  -  e ,  £  «  0 .  Satisfactory  agreement  of 
the  present  solution  with  known  solutions  is  found. 
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From  (3.4.2),  we  have 


,%  2 


42ttv_  ( x )  ..  /r —  ;«-3/4  2(0'  _  _  .  \ 

lim- - j= - =  lim  V2;r<?  3  — =^VL(<y)  = 

*-*°  V^r  Gi22Ti 


(5) 


(3.4.13) 


thus,  the  displacement  v(x,0,t)  can  be  expressed  as 


v(x,0,t)  = 


^12^1 


(3.4.14) 


For  a  finite  crack  extension  8,  the  following  Irwin’s  crack-closure  integral  (the  definition 
of  the  strain  energy  release  rate  for  single  mode  at  once  )  can  be  evaluated  [18, 19] 


„/x  f*<T>(x,0)v(£-x,0)j  K / 

GAt)  =  lim  — i - - - <£c  = - i — 

' +*Jo  6  2  GnX, 


(3.4.15) 


without  ambiguity  by  using 


Cs 

Jo  V  * 


15  l8-xdx^4K  T{\!Z)  _n 

X~  2  T(2)  2 


(3.4.16) 


Equation  (3.4.15)  further  reduces  to  the  well-known  result  for  isotropic  materials  [20],  i.e. 

(3.4,,) 

E 


4  In-plane  Shear  Impact 

Consider  the  crack  geometry  illustrated  in  Figure  2.  Figure  3  schematically  shows  the  same 
crack  with  a  spatially  uniform  shear  traction  of  magnitude  t0  applied  suddenly  on  the  crack 
faces.  Exploiting  symmetry  and  limiting  ourselves  to  upper  half-plane,  y>0,  the  corresponding 
boundary  conditions  become 
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(4.1) 


txy(x,0,t)  =  -t0H(t), 
<7).(*,0,t)  =  0 

Tyz(x,0,t)  =  0 


-oo<X<0  , 
-00<X<+00  , 
-00<X<+00  , 
0<X<+o»  , 


n(x,  0,  t)  =  0, 

where  H(t)  is  the  unit  step  function.  The  displacement  at  infinity  is  assumed  zero,  and  the  body  is 
stress  free  and  at  rest  everywhere  for  t  <  0 . 

A  method  of  solution  similar  to  that  that  used  for  normal  impact  can  be  used.  Applying  the 
Laplace  transform  to  boundary  conditions  (4.1)  gives 


Tx*  (x,  0,  t)  =  — T0  /  S  , 
cry*  (x,0,t)  =  0, 
Tyz\x,0,t)  =  0, 
u*(x,0,  t)  =  0, 


-o°  <  x  <  0, 

-OO  <  X  <  +°°  , 
-  °o  <  X  <  +°o  , 

0  <  x  <  +°° . 


(4.2) 


The  displacement  field  given  by  (3.2.3)  yields  the  same  system  of  ordinary  differential 
equations,  (3.2.4),  whose  solution  is  given  by  equations  (3.2.7)  where  af(ty,s)  and  y5((CJ,^)  are 

defined  by  (3.2.8)  and  are  obtained  from  the  solution  of  (3.2.10). 

Substitution  of  the  expressions  forcr*,  T*yz ,  the  first  and  third  of  equations  (3.2.15),  into  the 
boundary  condition  for  cr*,  T*z ,  the  second  and  third  of  equations  (4.2),  yields 


m,A,  +  m2A2  +  m3A3  =0, 
klAl  +  k2A2  +k3A3  =0. 


(4.3) 


These  equations  can  be  rewritten  as 

A2  =  Af  (o),  s)Ai ,  A3  =  (o),  s) A, , 


XI  ( a. >,  s)  =  m4L- % 


(4.4) 


m3k2  -  m2k3 


m3k2  —  m2k3 


As  ct)  -»  oo ,  we  have 
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(4.5) 


i //,  x  ,  ,//~  m\  ^3  ~m3  k\  3//^,  .  -)//“  mi  ^“2  m2  K 

4"  (<»,*)-+ A,  =  T— - — r,  /t'2  (<y,s)-+A2  - — - — 

m3  k2  —m2  k3  tn3  k2  -m2  k3 


where  mi  (a),  s)  —>  6)  m°° ,  n;  (<w,  n,~  and  (&),  s)  — »  /c,  ,as  w  . 

Using  this  result  in  the  expression  for  shear  stress,  the  second  of  equation  (3.2.15),  the 
expression  for  horizontal  displacement,  the  first  of  equations  (3.2.14),  and  the  remaining 
boundary  conditions,  we  obtain 


r(n,  +  4"n2  +  A>3)A,  -t0/sH(-x)  +  t+‘(x,s), 

2k  j-oo 

—  [(1 + Af  +  M!  )Ale~ia*da )  =  ul  (x,  s ) . 

2n  J-~ 

By  Fourier  transform  inversion,  these  equations  become 
—  Gj2  (n{  +  Af  n2  +  n3 )  At  =  —  t0  /itus  +  T+  (<y,  5), 

(1  +  4W  +  A?)A,  =U_(6>,s). 


(4.6) 


(4.7) 


where 


T+  {co,  s)  =  P  T+*  (x,  s)ei0*dx , 

U_(ty,j)  =  J°  u_*(x,,s)e““dlx. 

Eliminating  A,  from  equation,  (4.7),  we  obtain  a  Wiener-Hopf  equation, 
G„X ,  .  2=  rU_(a).J)  =  -fr— T »(«>,*), 

\lm1  +  s2/ic/  ■«* 

where  the  following  functions  are  introduced 


(4.8) 


(4.9) 


yl^+sI/icJ{nl+A]n2+X2n3) 

G{co,  s )  = - 7771—  r~r - 

Af2(l  +  ^  +  ^2) 


(4.10) 


Af2 


1+4"“+ AT 


(4.11) 


24 


where  %i  1S  dependent  on  6  and  material  parameters  and  its  expression  is  shown  in  the 
Appendix.  Here,  note  that  for  6  =  0 ,  Xi  =  — ~  >  a°d  for  6  =  90° , 


Y  90  _JL 
Xi  h ' 


(4.12) 


where 


8  —  (C11C33  C13  )c55  [C33  2N ^(N t  +  N2  +  ^3  )]  > 

h  =  clt{c33[cuNlN2N3+c55(Nl+N2+N3)]+ 

A^A^-c,,2  -2  cI3c55  +cnc55(NlN2  +  N3N2  +  N3Nl)]}, 

and  Nj  are  shown  in  equation  A.5  of  the  Appendix.  Xi  and  X™  310  the  limit  values  of  %r  at 

6  =  0  and  6  =  90° ,  respectively. 

Also,  for  £0  =  0,  after  lengthy  algebraic  manipulation  of  equation  (4.10),  we  have 


G(0,s)  =  - 


2  „ 
S  p 


KsiZlCs 


(4.13) 


where 


^VC55 


Equation  (4.9)  contains  only  the  two  unknown  functions,  T+(£0,j)  and  U_(£0,s)  ,  and  can  be 
solved  using  the  Wiener-Hopf  technique  as  in  the  normal  impact  analysis  of  section  3.  The  result 


is 


D+  (£0)  -  X+  (£0)L+  (£0)  =  L_  (£«)U_  (co)  -  D_  (eo )  =  0, 


(4.14) 


where 
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and 


L_  (fl»  =  GnX<i  -y  -  1^Cr  } 
Ja)-is/Ks] 

L+  (0))  -  — i  slKjl 

(a)  +  is/cR)G+  (of) 


G_(co), 


6(a))  = 


G(oj,s) 

a)2  +  s2 / cR2 


D(cd)  =  ^-L+(cq) 

ICOS 


tAL+(0))-L+(Q)  |  L+(Q) 
s  [  ia)  id) 


D+(a))  +  D_(a)), 


with 


0) 

ico 


D_  (ay)  -  — 
s 


(4.15) 


(4.16) 


Hence,  we  can  solve  for  the  transform  of  stress  ahead  of  the  crack  tip,  T+(<sj)  ,  and  displacement, 
U_  (&>) ,  behind, 


L+(co)  sico 

U  (of)  =  — 

L_(co)  sid) 


1- 


L+(0) 


L+m 

MO) 

L_(0)) 


(4.17) 


It  can  be  shown  thatG(6>)  — »  1  as  0)  — »  oo  (the  constant  %2  ln  (4. 1 1)  was  chosen  to  make  this 
possible).  Hence,  G±  (co)  is  given  by 
G(co)  =  G+(cq)G_(o>) 


Jr±  Z—O) 


■dz\. 


(4.18) 
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where  R_  is  the  closed  counterclockwise  contour  enclosing  the  branch  points 
+is/Ksl,  +  is/ns2  and+is/Kd  ,  and  also  R+  is  the  closed  counterclockwise  contour  enclosing 
the  branch  points  -is/ns , ,  -is/ns2  and  -is/Kd  . 


Note  that  as  a)  -+<»,  G+(co)  -+ 1  ,  G_ (a>)  1 


L+(0))  =  -^,  L_(co)  -  Gl2%24a  and 
Vty 


MO) 


/*K/*G+(Q)4i 


,  and  using  equation  (4.17),  we  can  readily  obtain  T+(ry)  and  U_(co) 


as  co  — >  <» 


'I0CR 


T+(ft>)  -» - p 

(siY^J^G+(0) 


(4.19) 


U_(<y)-> 


'I0CR 


Gnz2(sico/2rcs/iG+(0) 


(4.20) 


From  equations  (4.7)  and  (4.4),  the  coefficients.  A, ,  A2  and  A3 ,  as  6)  — >  °°  become 


A,  =  /"U_(<y)  A2  =  f"G_(co)  and  A3  =  /3"U_(<0) 


(4.21) 


where  /”  =  J  ....  .  /,” - ^ 


—  and  fJ  =- 

//“  ,  o//“  '  '  ‘  i  .  i//~  .  o//~  3  i  ,  i//~  ,  2ir 


1  +  4"  +  ^' 


i+4"  +4' 


1+4"  +4 


Substituting  (4.21)  into  (3.2.7),  the  dual  transformed  displacement,  w"  (a>,0,  s) ,  as  6)  — >  ©° 
becomes 


w"*(<y,0, j)  -+  W_" (a»  =  /ffU_(a>) 


(4.22) 


****  f*  =  PC  fx*  +  Pi~  fi  +  P*~  fj 


From  Abelian  theorem  [12, 17], 


limV*T  *(jc,s)  =  lim  e~'*/4J—T.(co). 

*— >o  y  ^ 


(4.23) 
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% 


limiLW  =  lime,„„2® 
*-*°  sx  ^ 


■u_m 


(4.24) 


w''  (a;) 


2(0 


X 


lim  "V"'  =  lim  einV4  ^=-  W_"  (©) 


x->0 


4~x 


4t? 


(4.25) 


Using  (4.21)  and  the  definition  of  the  stress  intensity  factor  gives 


42?, 


0 


K'u (5)  =  lim 42m?' (jc, s )  =  lim  e^'4 42aT+(co)  =  ■ 

saJ44g+(  0) 


(4.26) 


Note 


that  G+  (0)  =  G_  (0)  =  ^ G(0)  =  (  consequently,  equation  (4.26)  can  be 


rewntten  as 


K;,(s)  = 


42a Q 


(4.27) 


Using  equation  (4.13),  G(0,  s) 


2  „ 
s  p 


Ks\XlCs 


-,  it  readily  shown  that 


v*  _  V^o  [fsZ 2 

V  p  ' 


(4.28) 


Applying  the  inverse  Laplace  transform  to  K?„(s),  the  stress  intensity  factor  K„(t)  can  be 
written  as  follows 

lo  o  v~  I 

(4.29) 


where 


r  -  2c>Xi 

'■'ii  ~  \  — 

\  P 


(4.30) 


Thus,  the  stress  *^(*,0,0  can  be  expressed  as 
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(4.31) 


Txy(x,0,t)  = 


41 


'.7DC 


For  6  =  0,  i.e.  for  the  case  of  transversely  isotropic  material,  C„  reduces  to 


^  o  _  c~j  2(c„  -1) 

—  i  _ 

V  n  cn 


(4.32) 


By  using  cM  =  ^  ^  ,  equation  (4.32)  further  reduces  to  the  well-known  result  for  isotropic 


(l-2v) 


materials  [13],  i.e. 

C 


r  0  _  2Cf 

"  "V*(l- 


v) 


(4.33) 


For  6  =  90° ,  i.e.  for  the  case  of  orthotropic  materials,  C„  reduces  to 


V  K 


90 


(4.34) 


where  1S  shown  in  (4.12).  Tables  2  and  3  show  that  the  numerical  values  of  C,,90  for 
orthotropic  materials  are  same  as  those  obtained  by  using  formulation  in  [10],  hence,  equation 
(4.26)  also  reduces  to  the  result  for  orthotropic  materials  as  expected  for  6  =  90° . 

From  (4.24),  we  have 


V2/ru_  ( x )  r- —  .  3/4  2 (d^1  2  * 

hm - t= - =  hm  sjlne ,Ki'  — =-[/_(&)  = - K„ 


lim—  -  —  =  lim  ■sllKe"13"*  =z:=-U  (co)  =  K .,  (s)  (4.35) 

*■*  ~~  GnXl 


thus,  the  displacement  un(x,0,t)  can  be  expressed  as 


uV'°’')=^J£x"(,) 


(4.36) 


Similarly,  from  equation  (4.25)  the  displacement  w“(x,0,r)  can  be  expressed  as 
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(4.37) 


wn(x,0,0'= 


5  Anti-plane  Shear  Impact 

Consider  the  crack  geometry  illustrated  in  Figure  1.  Figure  4  schematically  shows  the  same 
crack  with  a  spatially  uniform  shear  traction  of  magnitude  i0  applied  suddenly  on  the  crack 
faces.  Exploiting  symmetry  and  limiting  ourselves  to  upper  half-plane,  y>0,  the  corresponding 
boundary  conditions  become 


lyzix,  0,t)  =  -z0H(t), 
cry(x,0,t)  =  0 

Txy  (*»  0,  t)  =  0 

w(x,  0,  t)  =  0 , 


-oo<  *<0, 
-OO  <  X  <  +°°  , 
-  oo  <  X  <  +oo  , 

0  <  X  <  +oo  , 


(5.1) 


where  H(t)  is  the  unit  step  function. 

The  displacement  at  infinity  is  zero,  and  the  body  is  stress  free  and  at  rest  everywhere  for 


t<  0. 


The  method  of  solution  is  similar  to  that  that  used  for  normal  impact.  Applying  the  Laplace 
transform  to  boundary  conditions  (5.1)  gives 


Tyz‘(x,  0,  t)  =  -T0/s , 
<Ty\x,0,t)  =  0, 
Txy*(x,0,t)  =  0, 


-oo  <  X  <  0  , 

-  oo  <  X  <  +oo  , 

-  oo  <  X  <  +oo  , 

0  <  X  <  +oo  . 


(5.2) 


w  (x,  0,  t)  =  0 , 

Assuming  the  displacement  field  given  by  (3.2.3)  yields  the  same  system  of  ordinary 
differential  equations,  (3.2.4),  whose  solution  is  given  by  equations  (3.2.7)  where  a;(co,s)  and 


y0;(ta,^)  are  defined  by  (3.2.8)  and  are  obtained  from  the  solution  of  (3.2.10). 
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Substitution  of  expression  for<7*,T^,  the  second  and  third  of  equations  (3.2.15),  into  the 


boundary  condition  for  tr*,  t*  ,  the  second  and  third  of  equations  (5.2),  yields 

n,A,  +n2A2  +  n3A3  =0 
ktA{  +k2A2  +  k3A3  0 


These  equations  can  be  rewritten  as 

A2  =  (CD,  s)A, ,  A3  =  X2  (a,  s)A, , 

V  (m,  .  Jlf  (<*»)  =  - 

n3  k2  —  n2k3 
As  co  -) <*> ,  we  have 


nxk2  —  n2kl 
n3k2  -  n2k3 


(5.3) 


(5.4) 


©o  •  oo  00,00 

1  III  / „  _N  v  -!///~  _  W1  *3  —  n3  *1  0/W/„  v  0///° 


(©,*)->  4  = 


oo  .  oo  oo  •  o 

n3  k2  -n2  k3 


-,  X"'(co,s)-*W  =- 


OO  J  OO  OO  7  O 

”l  ^2  ~n2  *1 


oo.oo  oo  * 

no  —  n- 


(5.5) 


where  mi  (CD,  s)~)  CD  m°° ,  ni  (CD,  s)-)  CD  n~  as  CD -)<*> . 

Using  this  result  on  the  expression  for  normal  stress,  the  secondary  of  equation  (3.2.15),  the 
expression  for  horizontal  displacement,  the  first  of  equations  (3.2.14),  and  the  remaining 
boundary  conditions,  we  obtain 


j7m,  +^"m2  ro3)A,  e  i(axdcD  =  -r0/sH(-x)+r+‘(i,i), 

2k  J^o 

—  f(A  +  M”  +  )\e~i(*dcD  =  w!(jc,5). 

In  J— 

By  Fourier  transform  inversion,  these  equations  become 

-G12(m,  +X!i"m2  +X l"m3)Al  =  -r0/icas  +  r+(CD,s), 

(fl  +M®  )4  =W_((D,s), 


(5.6) 


(5.7) 


where 


T+  (cd,  s)  =  T*  (x,  s)eie*dx , 

W_((D,s)  =  Jf>  w  *  (x,s)e,0Kdx. 
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Eliminating  A,  from  equation,  (99),  we  obtain  a  Wiener-Hopf  equation, 

t 

E(<0,s) 


- - . - W_  (<0,  s)  =  — — r+  (<y,  s ) , 


(5.8) 


where  the  following  functions  are  introduced 


+  s^/kJ  (m,  +  ^m2  +  ^m3 ) 

~*2(A+M+M> 


*3  l+^f+Af" 


(5.9) 


(5.10) 


where  is  dependent  on  0  and  material  parameters  and  its  expression  is  shown  in  the 
Appendix  .  Here,  note  that  for  6  =  0 ,  ,  and  for  6=90°,  X-T  =  when  using  N, 


shown  in  equation  A.5  of  the  Appendix.  Xi  and  Xi  are  the  limit  values  of  x2  at  6  =  0  and 
6  =  90° ,  respectively. 

For  <»  =  0,  after  lengthy  algebraic  manipulation  of  equation  (5.9),  we  have 


E(0,s)  =  - 


q 


Ks2%3  Cs 


(5.11) 


where 


(^C55  C55  ) 

q  =  —j= — 1 - r - - — • 

yc55  +  c55  sin  6  +  cos  6 

Equation  (5.8)  contains  only  the  two  unknown  functions,  T+  (<0,  s)  and  W_  (0J,s) ,  and  can  be 
solved  using  the  Wiener-Hopf  technique  as  in  the  normal  impact  analysis.  The  result  is 
D+(<0)-E+(<0)L+(<0)  =  L_(0))W_  (<0)  — D_(<0)  =  0.  (5.12) 

where 
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L- m  =  Gl2z3  E_ (cd) , 

V co-is/ks1 


K  (co) 


^J(0  +  is/fcs2 
(a)  +  is/cR)E+((D) 


E((D)  = 


E(fi^s) 
co2  +s2/cR 


2  ’ 


and 


D((o)  =  ^L+(co)  =  ^ 
ICOS  s 


'L+(n?)-L+ 

ico 


~M 0)  ,  L+(0) 
ico 


=  D+(co)  +  D_(o)), 


(5.13) 


with 


D+(co)=  — 


L+(co)-L,(  0) 


1(0 


,  D-  (co)  =  — 


MO) 

ico 


Hence,  we  can  solve  for  the  transform  of  stress  ahead  of  the  crack  tip,  r+(<M),  and 
displacement,  W_(co) ,  behind, 


L+(co)  sico 

D  (co)  Tn 
W_  (co)  =  — — — -  =  — — 
L_(co)  sico 


1- 


M 0) 


M«) 
M  Q) 

L_(co) 


(5.14) 


It  can  be  shown  that  E(<y)  — »  1  as  — » <»  (the  constant  in  (5. 10)  was  chosen  to  make  this 


possible).  Hence,  E±  (co)  can  be  given  by 
E(co)  =  E+(co)E_(co) 

A 

E±(co)  =  exp 


_Lris«e*  . 

2  mh±  z-co 


(5.15) 
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where  R_  is  the  closed  counterclockwise  contour  enclosing  the  branch  points 

t 

+is/ns] ,  +Is/ks2  md+is/Kd  ,  and  also  R+  is  the  closed  counterclockwise  contour  enclosing 
the  branch  points  -is/tcsl,  -is/ns2  and  -is/Kd  . 

Note  that  as  -»<*>,  E+(o))  -4  1  ,  E_ (ty)  -4 1  ,  L+(co)  =  -^=,  L_(co)  =  Gl2x3'Jco  and 

V<» 


L+(0)  =  - 


/2k/*E+(0)41 


,  and  using  equation  (108),  we  can  readily  obtain  T+  (co)  as  0)  -4  °° 


'I0CR 


r+(o})  -4 - y/ 

(siV2^^E+(0) 


(5.16) 


W _(&)-> 


*0  CR 


Gaisin)3/2  ks/zE+(0) 


(5.17) 


From  equations  (5.7)  and  (5.4),  the  coefficients,  A, ,  A2  and  A3 ,  as  (0  -4  «>  become 


A,  =  gf'W_(<y)  A2  =  g'"W _m  and  A3  =  g"'W _(o) 


All 


Illy 


(5.18) 


where 


tr = 


_  9  O'"  — _ 

-  1  ///**  62  n  - 


r+«"“+A-Af* 


and  _/«  _ 

*3 


Substituting  (5.18)  into  (3.2.7),  the  dual  transformed  displacement,  u®  (6),0, s) ,  as  <2)  -4 <» 


becomes 


um  (co,0,s)  -4  U‘“  (©)  =  gm W_(<y) 


///, 


(5.19) 


where  g,„  =  g("  +g'2  +  g"1 


From  Abelian  theorem  [12,  17], 


lim-y/xr  *(jc, 51)  =  lime  "r/4n|— r  (ty). 
x->o  y 


(5.20) 
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lim  !— 
x^°  4x 


=dime'*J,4^l 

<u->~  4k 


._  w_  (x)  _  __ijr3/4 


W  (©) 


(5.21) 


lim  U-"^  =  lim  e,Vr3/4  U"f  (fl>) 
x_>0  vx  m~*°°  4k 


(5.22) 


Using  (5.18)  and  the  definition  of  the  stress  intensity  factor  gives 

K*m (5) =  a/2 ~kx.t'{x,s)  =  lim  e~,n'^ =  _^2£^± - 

x_>0+  *"”*•  s/2Jk^E+(  0) 


(5.23) 


Note  that  E+( 0)  =  £_( 0)  =  -^£(0)  =  Cr^^(Q’s)  ,  consequently,  equation  (5.20)  can  be 


rewritten  as 


V2t0 


(5.24) 


Using  equation  (5.1 1),  E(0,  s)  = - — ,  it  readily  shown  that 

KslXz  Cs 


Cw  = 


V2t0  fcjfr 


,3/2 


(5.25) 


Applying  the  inverse  Laplace  transform  to  t£ufs) ,  the  stress  intensity  factor  A^/f)  can  be 
written  as  follows 

_  _  lo  ,  v  I 

(5.26) 


where 


Ml 

9 


Thus,  the  stress  t)>z(x,0,r)  can  be  expressed  as 


(5.27) 
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(5.28) 


For  6  =  0,  i.e.  for  the  case  of  transversely  isotropic  material,  Cw  reduces  to 


C  0  =• 
'-/// 


1 2-\/C55  c 


(5.29) 


By  using  c55  =  1  ,  equation  (5.29)  further  reduces  to  the  well-known  result  for  isotropic  materials 
[13],  i.e. 


(5.30) 


For  6  =  90° ,  i.e.  for  the  case  of  orthotropic  materials,  Cin  reduces  to 


r  90- 

'-■///  — i 


2c  *  -\/C55 


(5.31) 


Tables  2  and  3  show  that  the  numerical  values  of  C,,,90  for  orthotropic  materials  are  same  as 
those  obtained  by  using  formulation  in  [10],  hence,  equation  (5.26)  also  reduces  to  the  result  for 
orthotropic  materials  as  expected  for  6  =90°.  From  (5.29)  and  (5.31),  Cm  is  same  as  Cw  , 
which  is  because  d55  (0° )  =  d55  (90°) . 


From  (5.21),  we  have 


lim  (X)-  =  lim  W  (0))  =  — Kin\s) 

V*  GnX-i 


(5.32) 


thus,  the  displacement  wm(x,0,r)  can  be  expressed  as 


(5.33) 


Similarly,  from  equation  (5.22)  the  displacement  u 111  (x, 0,0  can  be  expressed  as 
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um(x,0,t) 


(5.34) 


6  Results  and  Conclusion 

Closed  form  expressions  of  the  dynamic  stress  intensity  factors  Ky(t) ,  K„(t)  and  K„,(t) 
have  been  determined  for  semi-infinite  cracks  in  rotated  transversely  isotropic  materials  under 
the  three  loading  modes,  i.e.  opening,  in-plane  shear  and  anti-plane  shear  loadings,  respectively. 
The  stress  intensity  factors  for  the  three  modes  are  proportional  to  the  square  root  of  time;  there 
is  no  equivalent  quasi-static  problem  for  semi-infinite  crack  under  uniform  loadings  and  there  is 
no  long-time  equilibrium  value.  The  method  of  solution  in  this  paper  differs  from  that  typically 
used  in  the  isotropic  case.  The  displacement  formulation  of  the  equations  of  motion  is  solved 
without  the  use  of  Helmholtz  potentials.  This  allows  solutions  to  be  found.  As  0  =  0 °  and 
0  =  90°,  we  can  use  the  formulae  in  [10]  for  the  coefficient  of  stress  intensity  factor  for  the  three 
modes  to  check  our  result.  The  results  of  the  comparison  are  shown  in  Tables  2  and  3.  The 
results  found  here  are  completely  the  same  as  those  in  Table  2,  as  expected.  Thus,  it  has  been 
shown  that  the  rotated  transversely  isotropic  formulation  includes  the  case  of  the  crack  located  in 
the  transversely  isotropic  plane,  the  case  of  the  fiber  along  the  crack,  and  the  case  of  isotropic 
materials  as  special  cases.  That  is,  for  each  loading  mode,  the  dynamic  stress  intensity  factor  for 
these  special  materials  is  recovered  from  the  rotated  transversely  isotropic  expressions  with  the 
proper  substitution  of  elastic  constants. 

The  calculated  results  for  C,  (0) ,  C„  (0)  and  Cai(0),  where  Ki{t,0)  =  2ciCi{0)4t ,  are 

shown  in  Figure  5-8  for  several  materials,  and  the  mechanical  properties  [16]  used  for  the 
analysis  are  given  in  Table  1.  In  these  figures,  it  can  be  seen  that  the  fiber  reinforcement  of 
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epoxy  leads  to  an  increase  of  the  stress  intensity  factor  for  all  three  inodes,  i.e.,  the  stress 
intensity  factors  K(t)  for  Graphite-Epoxy,  Eglass-Epoxy,  and  Boron-Epoxy  composites  are 
greater  than  the  corresponding  K(t)  for  epoxy  at  all  rotation  angles.  For  mode-I  and  mode-II, 
the  coefficient  of  stress  intensity  factor  increases  as  the  rotation  angle  9  increases,  and  at  9  =  0° 
and  9  =  90°  the  coefficient  of  stress  intensity  factor  reaches  its  minimum  value  and  maximum 
value,  respectively. 

Figure  7  shows  the  coefficient  of  stress  intensity  factor  for  mode-M  with  rotation  angle  9 . 
The  shape  of  the  curve  is  different  from  those  for  mode-I  and  mode-II.  Surprisingly,  at  6  =  0° 
and  9  =  90°  the  coefficients  of  stress  intensity  factor  are  same,  which  can  be  obtained  by  proper 
substitution  of  elastic  constants.  This  is  because  dS5(0°)  =  dss(90°),  i.e.  shear  on  the  principle 

axis  is  independent  of  rigid  body  rotation  through  90° .  Also  at  9  =  0°  and  9  =  90”  the 
coefficient  of  the  stress  intensity  factor  has  its  minimum  value.  The  point  of  the  maximum  value 
for  coefficient  of  stress  intensity  factor  is  between  9  =  0°  and  9  =  90° .  The  maximum  value  and 
its  corresponding  rotation  angle  are  clearly  dependent  on  the  material  properties. 

Figure  8  shows  the  coefficient  of  the  stress  intensity  factors  K,(t) ,  K„(t)  and  Km(t)  for 
isotropic  (epoxy)  and  orthotropic  (graphite-epoxy)  materials,  note  that  with  same  load  amplitude 
the  coefficient  of  stress  intensity  factor  K„  (t)  and  Km(t)  are  greater  than  K,(t) .  By 
comparison  of  both  cases,  it  can  be  seen  that  the  introduction  of  fibers  results  in  an  increase  in 
the  stress  intensity  factor.  This  effect  is  more  dramatic  in  the  shear  modes  than  in  the  normal 
mode. 

In  terms  of  the  penny  shaped  crack  problem,  these  results  offer  some  intuitive  insights.  If  the 
crack  is  loaded  in  mode  I,  Figure  5  leads  us  to  expect  the  crack  front  to  advance  in  the  direction 


38 


of  the  fibers,  provided  the  toughness  is  constant  in  all  directions.  If,  however,  the  crack  is  loaded 
in  shear,  the  problem  is  more  complex  because  the  crack  front  sees  mixed  shear  mode  loading. 
There  are  components  of  both  mode  II  and  mode  III  loading  along  the  crack  front.  Letting  the 
shear  loading  with  magnitude,  ?0,  be  applied  on  the  surface  of  the  penny  crack,  and  letting  (p  be 
the  angle  of  that  shear  loading  relative  to  the  principle  axis  of  the  material  (fiber-orientation)  as 
shown  in  Figure  9.  The  shear  loading  at  a  point  6  on  the  crack  front  is  given  by  a  mode  II  and 
mode  III  component, 

^//=^o  cos(<p  —  6)  (6  1) 

^  in  =  t0  sm((p — 6) 


In  this  case,  due  to  the  presence  of  both  modes  at  once,  hence  there  are  cross-influences. 

For  a  finite  collinear  crack  extension,  6,  the  following  Irwin’s  crack-closure  integral  (the 
definition  of  the  strain  energy  release  rate  for  both  modes  at  once)  can  be  evaluated  [18, 19] 


,•  f5?  (*,0)[un(£-x,0)  +  uni(£-x,0)] 

(6.2) 

„  ,,  f*  r„  <1, 0)[w m (S  -  x,0) +»“<«-  x,0)] 

G„, «)  =  l.m  _[  -2 - - - * 

(6.3) 

Using  (3.4.16),  and  the  ?^(;c,0),  i yz (x,0) ,  un(x,0),  um(x,0),  w 11  (*,0)  and 

wni(x,0)  derived 

above,  G„  ( t )  and  G,„  ( t )  are  expressed  as 

G„(t)=  1  K„2+  8,11  KnKm 

2Gn%2  2GuZ3 

(6.4) 

G,„«)=  1  K,„2  +  J"  K„K„, 

2GnZ3  2GnXi 

(6.5) 

where  equations  (6.4)  and  (6.5)  further  reduce  to  the  known  result  for  isotropic  materials  [20], 

i.e. 
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2 


(6.6) 


G„( 0  = 


d-v2)^ 


G,„(0  = 


(l+v)*,,,2 


(6.7) 


Substituting  the  expression  for  K„  (/)  and  K (r)  derived  here,  (4.29)  and  (5.26),  the  strain 


energy  release  rates,  Gu  (t)  and  G,„  (t) ,  become 


G„(t)  = 


4c  t  r„  4c J 


_  »// 


Gyjt  P  G^Jt 


8111*11*111. 


X  2 

XiPq 


(6.8) 


4c«r  t,„  Acj  r  I  Xi 
J  //  *  //  *  w  ^ 


G,„(t)  =  *-  +  -'5 


The  total  energy  release  rate  is  given  by 

G total  (0  =  Gt,  (t)  +  G,„  (t) , 


ZzPV 


thus, 


.  2  ^  total 

4csr0  r 


(0 


cos  2{(p-8)  sin  2{(p-0) 

= - 1 - b 


(6.9) 


(6.10) 


cos(9>-0)sin(9>-0) 


8ni. 


1  Xi  ,  f 
XiPq  "\X2Pq 


(6.11) 


For  a  given  shear  loading  angle,  <p , 


Gnn 


4csf0  t 


maximum  energy  release  rate,  max 


Gn*  r  a\ 

2  ^ total  CO 


4cst0  t 


G toiai  (0  is  a  function  of  only  6 .  Thus,  the 


,  occurs  at  6  =  6  ( 0  <6  <n ) )  and  6  is  a 


function  of  <p .  The  functions  are  periodic  in  (p  and  6  with  a  period  of  n  . 

Figure  10,  1 1  and  12  show  the  normalized  strain  energy  release  rates  for  0 <6  <  180°  and 
0  <  <  180°  in  Graphite-Epoxy,  Eglass-Epoxy  and  Boron-Epoxy,  respectively.  There  are  two 
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different  peaks  for  given  (p  in  each  of  these  figures,  where  the  loading  is  most  severe. 
Meanwhile  we  also  find  valleys  in  these  figures  where  the  loading  is  least  severe.  The  most 
probable  direction  of  propagation  for  each  case  and  the  corresponding  value  of  the  normalized 

71 

strain  energy  release  rate  are  shown  in  Figures  13  and  Figure  14.  It  is  seen  that  6C  ~  —  +  (p  for 

each  material,  that  is,  the  crack  will  most  likely  propagate  in  the  direction  perpendicular  to  shear 
loading. 
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Appendix 


(1)  Roots  of  the  characteristic  equation,  (3.2.10) 

The  solutions  of  the  cubic  equation  (3.2.10)  of  y2,  y,2 ,  y  2  and  y32,  are  related  to  a,  b,  c  and  d 
by  the  following  relations, 


2  2  2  0  222  22  22  2  2  C 

r,  +r t  +Yi  = — .  TiTiY,  = — .  n  r2  +r2  r3  +r.  r3  =-• 

a  a  a 


(A.l) 


where  these  expressions,  y,-  ,  can  be  given  as 


and, 


2  _  2 s2/cs2  +(1  +  css)co2  +  (1  - c55)co2 cos20 
*  2 
2  (2c55  -2cu)s2/c2  +(Bt  +  B2cos26)co2  -<Jb^ 

y2  = 

^°11C55 

2  (2c55  -2 cu)s2/c2  +(Bl  +B2 cos26)co2  +J~B ^ 

r3  = - - - 

4cucS5 


—  C,J  2c13C55  +CuC33  +2cuc55, 


J?2  —  c13  c,,c33  +  2cnc55  +  2c13c55  , 

fi3  =  [2c55  $2/c,2  -c13  ft>2  -  2c13c55  +  c,,(2$2/cj2  +c33G)2  +2c55<»2)  +  B2®2cos2  2^]2  _2cucj 
[8^  /c,  +4(C|,  +c33  +  2c 55 ) co  s  jcs  +  (*■* cJ3  +c,|C33  +  3C|,cSJ  —  2c]3c55  +3c33c5J  + 
4(c,,  — c3i)co2(s2/c 2  +  c550)2)cos29  +  (c132  —  cnc33  +  cnc55  +2c13c55  +c33c55)fi>4  cos4# 


As  co  =  0,  equation  (3.2.10)  becomes 


a6  br\  ()4  Crt  a2  dn  _ 

y  +—y  +—7  +— =  0 
a  a  a 


(A.2) 


where 


a  =  d  (d  2  -d  d  )  =  -C|l^n — £12! 

U  “22VU4«  U44U66  '  2  ’ 
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.  _ft  1  J  2  ,  ,  .  .  x  •y2  _  (c|l  CI1C12  3Ci  |  Css  £12^55) 

^0  —  (^22^44  **46  ‘  “22^66  **■  ^44^66 )  2  —  o 

C  2  C. 


--(A  ±A  4- A  \S*  -  (~3C>I  +C|2  ~2Css)  ^ 

-  (a22  +«44  +«6«)  4-  o  4’ 


4>=-t> 


Three  roots  of  equation  (A.2)  are  given  as. 


Y, 


0  2  s  s 

Yi  =  ,7 - r— =  — 

Ucn-cl2)  cs  cs 


0  s  s 
Yi  ~  I —  — 

Vcii  cs  cd 


We  note  that  the  three  roots  are  independent  of  6 . 


(2)  Solution  of  the  characteristic  equation  for  co  -»  °°  (3.2.12) 
For  any  6 ,  we  have 

K1 2  _  l  +  css  +(1-cS5)cos20 

N\  “  o  ’ 


n22  = 


Bx  +  B2  cos  29  -  2^  sin2  6 


4C1IC55 


(A.3) 


N32  = 


Bx  +  B2  cos  26  +  2  *[bI  sin2  6 


4cnc55 


where  Bx  and  B2  are  given  above 

54=(c,3  —  cnc3'i  )(cn  ~cwci'i  +  4CJ3C55  +  4c55  ). 


Note  that  for  the  case  6  —  0 , 

n2  =  n2  =  n32  =  1 , 

and  for  the  case  6  =  90° , 


2C  llc55 


which  are  same  as  those  in  [10]. 


(A.4) 


(A.5) 


(3)  Expression  for  %x 

After  much  algebraic  manipulation  of  equation  (3.2.24),  it  is  found  that  %\  can  be  written  as 

X,  =—  (A.6) 

Xu 

where 

Xu  =(K~K~ -K-K-K-KSN^N,  +K~(Nl  +  N2+N3)]  + 

a,  {-K4~  [— AT6~  2V,  AT22V3  (//, 2  +  N2  +  N32  +  NxN2  +NxN3  +N2N3)  + 

^5“(N,3  +N2  +N3  +  N1N2N3+NxN2(N1+N2)  +  NxN3(Nx  +N3)  +  N3N2(N3  +N2))] 

+  K~[K~Nl2N22N32(Nl  +N2  +N3)  +  K~(NlN2N3(Nl  +N2  +  N3)  +  2NxN2N3(NxN2  + 
NXN3  +  N2N3)  +  N2N22(Nx  +  N2)  +  N2N2(Nx  +  N3)  +  N2N2(N3  +  N2))]}, 

Xu  =a,(«2^4~  +K~K~\nxN2N3{Nx  +  N2)(N2  +N3)(N1  +N3)\ 
and, 

al=d22{.d46  ~dl^d66), 

a2  =  ^12^46  ~  ^25^66  > 

/sT,  —  ~dx2diA+dxxd22d^—dx3d22d^()+dl2d2SdA()+2dX2di(>  —  2dx2d44d66, 
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2 

K2  =  ^12^15^46  —  ^11^25^46  —  ^11^46  "**^11^44^66’ 

^3  =~d\2  ^46  *^^11^22^46  “^15^22^66  ^"^12^25^66’ 

^4  =  ^12^15^66  —^11^25^66  ’ 

^5  =  ^66^11^25  —  ^12^15  )  > 

/STfi  =  — rf|2  ^4$  ■^'^11^22^46  “"^15^22^66  ^  ^12^25^66  » 

^7  =  “^12^15  +^11^25  *^11^46  ”^15^66  ’ 

the  expressions  for  JVf  are  shown  in  equation  (A.3)  in  the  Appendix. 


(4)  Expression  of 

After  tediously  long  algebraic  manipulation  of  equation  (4.11),  it  is  found  that  £2  can  be 
written  as 


(A.7) 


where 


*2„  =  (L,~L3~  +L2-L4~)[-L4~AWV3  +  L5“(AT,  +  W2  +AT3)]  + 

bx  i(L“)2Nl2N22N:}2  (iV,  +1V2  +/73)  +  (L3“)2[iV13  +  W23  +  W33  + 
W,Ar2lV3  +JViJV2(Ar,  +  A72)  + +  N3)  +  N3N2(N3  +  N2)]  + 
L~L~[NlN2N3(NlN2  +  NtN3  +  N2N3 )  + 
tf,2JV22(Ar,  +  W2)  +  +  W3)  +  N32N2(N3  +  N2)]], 
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XU  =-L1“[L4~[^^.2^22^32+V(^^2+^.^3+^2^3)]  + 

L3~[L~(NlN2+NiN3+N2N3)  +  b2[NlN2N3(Nl  +  N2+N3)+ 

N2N22  +N*N32  +  N32N2]]]-bl[L4m’[-a2N2N22N2(NlN2  +  NtN3  +  N2N3)+ 
L~[NlN2N3{Nl  +N2+N3)  +  2NlN2N3(NiN2+NlN3  +  N2N3)  + 

N2N2  (AT,  +  N2)  +  N2N32(Nl+N3)  +  N32N22  (N3  +N2))  +  L~[  L5~ 

[N2N22  +  N2N2  +  N2N2  +  (N2  +  N22)NtN2  +  (N2  +  N^N.N,  + 
(N32+N2)N3N2+2NlN2N3(Ni+N2+N3)]-b2[Nl3N2+Nl3N33+N2N33 
+  N2N2N3(Ni  +  N2)  +  N.N^N^N,  +  N2)  +  ^12AT2^32(iV1  +  N3)]]  + 
L2~[-L4~[L6~  +b2NlN2N3{Nl+N2+N3)]  +  L3~[L~  +b2(N2+N2+N32  + 
NlN2  +  NlN3  +  N2N3)]], 


~  ^22  (^46  ^44^66  )  > 

^2  =  ^22^46  » 

oo  2  2 

Lj  =  (^12  ^44  "”^11^22^44  “*"^15^22^46  ““^12^25^46  “”^12^46  +  2rfj2^44^66 )  » 

oo  2 

^2  =  ^12^15^46  ”^11^25^46  —<^1 1^46  +^11^44^66’ 

•^3  =  ^66^11^25  —  ^12^15)  ’ 

L4  =~dn  dtf  +dud22d46  *"^15^22^66  ‘*'^12^25^66  » 

=  d\3d22  —dftdft  ~dl2dA6  ~d2sd66 , 

^6  =  d\sd(fi , 

the  expressions  for  A/,,  are  shown  in  equation  (A.3)  in  the  Appendix. 


(5)  Expression  for  ^3 

Finally,  the  expression  for  can  be  written  as 
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z3n  =-(Mn2NlN2N3[-M~  +  cxNxN2N3(Nx+N2  +N3)]  +  M~M~[-M~NlN2N3  + 

M  +N2  +  N3)  +  cx[N2N2(Nx  +  N2)  +  N2N2(Nx  +  N3)  +  N2N2(N3  +  N2) 

+  NxN2N3(NxN2+N1N3+N2N3)]]-(M~)2[M3~(Nx  +N2  +  N3)  +  cx(Nx3  +N2  +  N3 
+  Nx  N2N3  +  NXN2(NX  +N2)  +  NxN3(Nx  +  N3)  +  N3N2(N3  +  N2 ))], 

Xu  =-M~M2[M~  -M^°{NxN2  +nxn3+n2n3)-c2(nx2n22 +n2n2  +n2n2)]  + 
(M~)2[c2N2N22N2  +M~(NxN2  +NxN3  +  Ar2yv3)]  +  (M2“)2[-M6“  + 
c2(N2  +  N2  +  N2  +  NxN2  +  NxN3  +  N2N3] 

+  NxN2N3  +  NxN2  (Nx  +N2)  +  NxN3  (Nx  +N3)  +  N3N2  (N3  +  N2 ))], 

and, 

C\  ^  ^22^46  ^44^66  )  ’ 

C2  “  ^22^66  ’ 

M\  =““^12  ^46  +  ^11^22^4 €  ~  ^15^22*^6 6  +  ^12^25^66’ 

M 2  =  ^66(^12^15  “^11^25)  » 

Af3  =—  (rf12  ^44  “"^11^22^44  +  ^15*^22^46  “^12^25^46  “^12^46  +2^12^44^66)  , 

^4  =  ^12^15^46  ”  ^11^25^46  “  ^11^46  +  ^11^44^66  ’ 

^5  =^11^66» 

^6  =  “(^12  "”^11^22+^^12^66)’ 

the  expressions  for  N,  are  shown  in  equation  (A.3)  in  the  Appendix. 


Coefficient  of  stress  intensity  factor  ^  (m/s)  -  Coefficient  of  stress  intensity 


0  10  20  30  40  50  60  70  80  90 


Rotation  angle  0° 

Figure  7  Coefficient  of  stress  intensity  factor  for  Mode-Ill 
for  a  semi-infinite  crack  in  different  materials 


0  10  20  30  40  50  60  70  80  90 

Rotation  angle  0° 

Figure  8  Coefficient  of  stress  intensity  factor  for  three  modes 
for  a  semi-infinite  crack  in  different  materials 
The  isotropic  corresponds  to  Epoxy  and  the  orthotropic  to  Graphite-epoxy 
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Table  1  Material  parameters 


^EEESESSBSBI 

£,  (Gpa) 

10.41 

12. 

19. 

5.21 

E1  (Gpa) 

156.75 

45. 

207. 

5.21 

0.49 

0.19 

0.21 

0.328 

^31 

0.31 

0.19 

0.21 

0.328 

Gn  (Gpa) 

3.49 

5.04 

7.85 

1.96 

G, j  (Gpa) 

7.07 

5.5 

6.4 

1.96 

p(Kg/m3) 

1580 

2100 

1990 

1260 

Table  2  Coefficient  of  stress  intensity  factor  for  three  modes 
at  6 =0°  and  6  =  90"  in  different  materials  by  using  of  the  formulation  in  [10] 


(with  proper  substitution  of  elastic  constants) 


Graphite-Epoxy 

Eglass-Epoxy 

Boron-Epoxy 

Epoxy 

o 

o 

11 

i 

26.65 

27.37 

31.02 

24.45 

n 

37.68 

34.43 

39.20 

34.38 

m 

36.70 

32.10 

33.79 

28.18 

o 

8 

ii 

i 

34.56 

29.36 

32.35 

24.45 

n 

75.24 

49.96 

77.33 

34.38 

m 

36.70 

32.10 

33.79 

28.18 

Table  3  Coefficient  of  stress  intensity  factor  for  three  modes  at  6 =0°  and  6  =  90"  for  several 


materials  by  using  of  the  formulation  in  the  paper  with  6  =  0°  and  6  =  90" 


Graphite-Epoxy 

Eglass-Epoxy 

Boron-Epoxy 

Epoxy 

II 

o 

© 

I 

26.65 

27.37 

31.02 

24.45 

II 

37.68 

34.43 

39.20 

34.38 

m 

36.70 

32.10 

33.79 

28.18 

(9  =  90" 

i 

34.56 

29.36 

32.35 

24.45 

ii 

75.24 

49.96 

77.33 

34.38 

m 

36.70 

32.10 

33.79 

28.18 
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Figure  10  Normalized  strain  energy  release  rate  for  Graphite-Epoxy 


57 


Figure  1 1  Normalized  strain  energy  release  rate  for  Eglass-  Epoxy 
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Figure  12  Normalized  strain  energy  release  rate  for  Boron-  Epoxy 


-90  -60  -30  0  30  60  90  120  150  180  210  240  270 


Angle  of  the  shear  loading  direction  <p 

Figure  13  Maximum  strain  energy  release  rate  as  a  function  of  shear  loading  angle,  <p 


Angle  of  the  shear  loading  direction  cp 


Figure  14  Direction  of  maximum  strain  energy  release  rate,  6C ,  measured  relative  to  the  X,  axis 


